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I. INTRODUCTION 

The observation of the top quark by the CDF and D0 
collaborations in the spring of 1995 [1,2] was the culmi- 
nation of a long and intensive search that began following 
the discovery of the r lepton in 1976 [3] and the bottom 
(b) quark in 1977 [4]. The discovery of these two particles 
gave a firm foundation to the existence of a third family, 
originally proposed by Kobayashi and Maskawa in 1973 
to account for the occurrence of CP violation within the 
standard model [5]. The b quark was shown to possess a 
charge of Qb = — \e. [6-8] and a weak isospin of I3 = — ^ 
[9-11]. Within the standard model (SM), this demanded 
the existence of a partner to the b quark with a charge 
of +|e and a weak isospin of +5. This partner is called 
the "top" quark. 

Initial searches for the top quark were carried out at 
e + e~ colliders. These searches looked for a narrow reso- 
nance (if a bound tt state was produced), an increase in 
the rate of e + e~ — ► hadrons (if a bound tt state was 
not produced), or events with more spherical angular 
distributions which differentiate top quark events from 
the more planar angular distributions expected from the 
lighter quarks. As shown in Fig. 1(a), experiments at the 
PETRA [12,13], TRISTAN [14], and SLC/LEP [15,16] 
colliders raised the lower limit on the top quark mass 
{ra t ) from 15 GeV/c 2 in 1979 to 45.8 GcV/c 2 in 1990. 
In the late 1980's, in the absence of a signal, the focus 
of the top quark search shifted from e + e~ colliders to pp 
colliders and higher center-of-mass energies. Unlike e + e~ 
colliders, pp colliders cannot provide direct limits on the 
mass of the top quark, but rather upper limits on the 
tt production cross section. By assuming a relationship 
between mass and cross section (as provided by SM the- 
ory), these cross section upper limits can be turned into 
lower limits on the mass. The UAl collaboration pro- 
vided the first such limit in 1988, setting a lower bound 
on the top quark mass of 45 GeV/c 2 [17]. This limit 
was followed in 1990 by an updated limit from UAl (60 
GeV/c 2 ) [18] and new limits from UA2 and CDF (69 [19] 
and 77 [20] GeV/c 2 respectively). In 1992, CDF raised 
the lower limit on the top quark mass to 91 GeV/c 2 [21], 
and in 1994, D0 set a lower bound of 128 GeV/c 2 [22]. 

The first evidence for tt production was claimed by the 
CDF collaboration in April of 1994 [23]. With an inte- 
grated luminosity of 19.3 pb -1 , CDF observed twelve 
candidate events with an expected background of about 
six events and estimated a 0.26% probability for the back- 
ground to fluctuate to at least twelve events. The ex- 
cess was assumed to be due to tt production and the 



cross section was determined to be ov 



13.9 



+6.1 

-4.8 



pb for 



m t = 174 GeV/c 2 . The D0 analysis in mid-1994 [24] 
based on 13.5 pb -1 yielded 7 events with an expected 
background of 3.2 ± 1.1 events. The D0 and CDF sen- 
sitivities (expected number of events for a given cross 
section) and expected significance (signal to background 
ratio) were the same. The small excess seen in D0, if 
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FIG. 1. (a) Lower limit on the top quark mass from 1978 to 
1994 [12-22]. (b) Published quark cross section results from 
1994 to 2001 [23,1,2,25-27]. The solid triangle marker with 
the dashed line uncertainty corresponds to the unpublished 
D0 tt cross section in mid-1994 [24]. 



interpreted as being due to tt production, gave a cross- 
section of 6.5 ±4.9 pb for m t = 180 GeV/c 2 . At the time 
of the top quark discovery the following year, the CDF 
and D0 collaborations reported tt production cross sec- 
tions of a tl = 6.8±|;^ pb for m t = 176 GeV/c 2 [1] and 
a tl = 6.4 ± 2.2 pb for m t = 199 GeV/c 2 [2], respec- 
tively. These results were updated by D0 (1997) and 
CDF (1998) to o tl = 5.5 ± 1.8 pb [25] for m t = 173.3 
GeV/c 2 and a ti = 7.6±i;| pb [26] for m t = 175 GeV/c 2 , 
respectively. In 2001, the CDF collaboration reported 
aa = 6-5ti;4 pb for m t = 175 GcV/c 2 [27] as their final 
tt production cross section based on the 1992-1996 run 
of the Tevatron. The corresponding result from the D0 
collaboration, reported in this article, is a t f = 5.7±1.6 pb 
for m t = 172.1 GeV/c 2 . 

At the Tevatron center-of-mass energy of 1.8 TeV, top 
quarks can be produced singly or in pairs. The two cross 
sections are of similar magnitude [28] but single top quark 
events are much more difficult to distinguish from back- 
ground and have not yet been observed [29,30]. This 
paper is thus concerned only with tt pair production. 

The pp — > tt production cross section can be factorized 
in terms of the parton-parton cross section and the par- 
ton distribution functions for the proton and anti-proton, 
and is written [31] 



a(pp -> tt) 



*,3 



dx l dx j ff{x ll ii 2 )f I '{x j ,n' 1 )a ij {s,^ 2 ,mt), (1.1) 
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FIG. 2. Lowest order Feynman diagrams for production of 
tt pairs at the Tevatron. At Tevatron energies, the diagram 
involving quark-antiquark fusion dominates over those involv- 
ing gluon-gluon fusion. 



where the summation indices i and j run over the light 
quarks and gluons, Xi and Xj are the momentum frac- 
tions of the partons involved in the pp collision, ff (x^, p 2 ) 
and fj(xj,p 2 ) arc the parton distribution functions, 
and (Tij{s, p 2 ,m t ) is the parton-parton cross section at 
s = XiXjS. The renormalization and factorization scales, 
typically chosen to be the same value p, are arbitrary 
parameters with dimensions of energy. The former is in- 
troduced by the renormalization procedure and the latter 
by the splitting of the cross section into perturbative (it) 
and nonperturbative (f p , f?) parts. An exact calculation 
of the cross section would be independent of the choice of 
p, but current calculations are performed to hnite order 
in perturbative QCD and are thus dependent on jjb, which 
is usually taken to be of the order of m t . Theorists typ- 
ically estimate the uncertainty introduced by truncating 
the perturbation expansion by varying p over some arbi- 
trary range, usually mt/2 < p < 2mt (the range used for 
all theoretical cross sections referred to in this paper). 

In leading-order QCD (LO), 0{a 2 ), tt production pro- 
ceeds through qq — > and gg — > tt processes (see Fig. 2). 
At y/s = 1.8 TeV, the qq — > tt process dominates, con- 
tributing 90% of the cross section with the gg — > tt pro- 
cess contributing only 10%. The first calculations of the 
LO cross section a were performed in the late 1970's 
[32-37] . Calculations of the tt production cross section at 
next-to-leading order (NLO), O{o? s ), began to appear in 
the late 1980's [38-44]. The 1990's saw the introduction 
of calculations which attempt to estimate the contribu- 
tion of the higher order terms through a technique known 
as resummation, in which the sums of the dominant loga- 
rithms from soft gluon emission to all orders in perturba- 
tion theory are calculated, thus reducing the dependence 
of the cross section on the value of p. The first such cal- 
culations [45,46] summed only leading- log (LL) contribu- 
tions. Increased precision was soon achieved through cal- 
culations [47,48] which incorporated summations through 
next-to-leading- log (NLL) contributions. The most re- 
cent calculations [49,50] sum contributions through next- 
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FIG. 3. Expected distributions for ep dilepton events 
of (a) electron Et or muon pr, (b) I$t, and (c) lep- 
ton g = tanh _1 (cos#) (two entries per event). The 
solid histograms are tt — ► ep + X signal events (gener- 
ated with herwig [54] with m t = 175 GeV/c 2 for pp col- 
lisions at ^J~s — 1.8 TeV). The dashed histograms are 
Z + jets — > tt + jets — > ep + jets events (also generated 
with herwig). All histograms are normalized to unity and 
all events are required to have p e T > 10 GeV/c, Ifr > 10 
GeV, and at least two jets with Et > 15 GeV and \g\ < 2.0. 

to-next-to- leading-log (NNLL). Although the NLL and 
NNLL calculations have reduced the scale dependence, 
kinematic-induced ambiguities lead to estimated unccr- 
taintied of about 7% (these latter uncertainties arc not 
included in the theoretical cross section predictions given 
in this paper). 

In the SM, the top quark is expected to decay pre- 
dominantly into a W boson and a b quark. Decay mecha- 
nisms whereby the top quark decays into a charged Higgs 
boson are not considered here, but are investigated in 
Rcfs. [51-53]. The channels in which the top quark is 
sought are thus determined by the decay modes of the 
two W bosons in the tt event. The W boson can decay 
lcptonically into an electron, muon, or a r lepton (and 
associated neutrino), and hadronically into ud, us, lib, 
cd, cs, or cb pairs. 

The channels can be classified as follows: the dilepton 
channel where both W bosons decay leptonically into an 
electron or a muon (ee, p/j,, ep), the lepton + jets chan- 
nel where one of the W bosons decays leptonically and 
the other hadronically (e+jets, /z+jets), and the all-jets 
channel where both W bosons decay hadronically. This 
paper will focus primarily on the dilepton and lepton + 
jets channels. The all-jets channel is discussed in detail 
in Ref. [56] and is only summarized here. The tt channels 
containing a tau lepton are not explicitly considered, al- 
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FIG. 4. Expected distributions for lepton+jets events of 
(a) electron Et and muon pr (two entries per event), (b) 
Ifr, and (c) lepton r\. The solid histograms are ti signal 
events (generated with herwig with m t = 175 GeV/c 2 for 
pp collisions at y/s = 1.8 TeV). The dashed histograms are 
W+ > 4 jet events (generated with VECBOS [55]). All his- 
tograms are normalized to unity and all events are required 
to have p e T > 15 GeV/c, Ifr > 15 GeV, and at least four jets 
with E T > 15 GeV and \r/\ < 2.0. 



though events containing r — > evv and r — » jivv decays 
do contribute to the efficiency of all channels containing 
an electron or a muon. Similarly, the inability to distin- 
guish between a hadronic tau decay and a hadronic jet, 
contributes to the efficiency of the lepton+jets channels. 
As is indicated in Figs. 3 6, the leptonic channels are 
characterized by high transverse- momentum (pr) leptons 
and jets as well as missing transverse momentum ($t) 
due to highp-r neutrinos (see Sec. IV D). The plots show 
the distributions of several kinematic quantities expected 
from ti decay compared with those expected from the 
leading background for the e\x (Figs. 3 and 5) and lep- 
ton+jets (Figs. 4 and 6) channels. Initial search strate- 
gies are based on previous studies and analyses [57,23,58]. 

The paper is structured as follows: Sec. II gives a 
brief overview of the D0 detector and indicates those 
aspects which were employed in the dilepton and lep- 
ton+jets analyses. Section III describes the triggers used 
in the first stage of the event selection. Event recon- 
struction and particle identification are the subjects of 
Sec. IV. Section V discusses the simulation of the ti sig- 
nal and background. The dilepton channels are described 
in Sec. VI and the lepton+jets channels are described 
in Sec. VII. The all-jets channel is described briefly in 
Sec. VIII. Section IX discusses the systematic uncertain- 
ties. The ti cross section results are summarized and tab- 
ulated in Sec. X and the conclusions to be drawn from the 
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FIG. 5. Expected distributions for e/i dilepton events of 
(a and b) the transverse energies of the two leading jets and 
(c) the jet r\ (two entries per event). The solid histograms 
are ti — > e/i + X signal events (generated with herwig with 
m t = 175 GcV/c 2 for pp collisions at y/s = 1.8 TeV). The 
dashed histograms are Z + jets — > rr + jets — > e/i + jets 
events (also generated with herwig). All histograms are 
normalized to unity and all events are required to have 
p e T > 10 GeV/c, If T > 10 GeV, and at least two jets with 
E T > 15 GeV and \r]\ < 2.0. 
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FIG. 6. Expected distributions for lepton+jets events of 
(a-d) the transverse energies of the four leading jets and 
(e) the jet r\ (four entries per event). The solid histograms 
are tt signal events (generated with herwig with m t = 175 
GeV/c for pp collisions at y/s = 1.8 TeV). The dashed his- 
tograms are W+ > 4 jet events (generated with vecbos). All 
histograms are normalized to unity and all events are required 
to have p e T > 15 GeV/c, tfr > 15 GeV, and at least four jets 
with E T > 15 GeV and \q\ < 2.0. 

combined analyses are presented in Sec. XI. Appendix A 
describes the corrections applied to the jet energy scale; 
Appendices B and C describe the main-ring veto and re- 
covery; Appendix D presents an independent Neural Net- 
work based analysis of the e/z channel; and Appendix E 
describes in detail the handling of the uncertainties and 
the correlations between them. 



II. THE D0 DETECTOR 

D0 is a multipurpose detector designed to study pp 
collisions at high energies. The detector was commis- 
sioned at the Fermilab Tevatron Collider during the sum- 
mer of 1992. The work presented here is based on approx- 
imately 125 pb _1 of data recorded between August 1992 
and February 1996. A full description of the detector 




FIG. 7. Cutaway view of the D0 detector, showing the 
tracking chambers, calorimetry, and muon system. 

may be found in Ref. [59]. This section describes briefly 
those properties of the detector that are relevant for the 
tt production cross section measurements. 

Spatial coordinates are specified in a system with the 
origin at the center of the detector and the positive z- 
axis pointing in the direction of the proton beam. The 
x-axis points radially out of the Tevatron ring and the y- 
axis points upward. Due to the approximate cylindrical 
symmetry of the detector, it is also convenient to use the 
variables r (the perpendicular distance from the beam- 
line), <j) (the azimuthal angle with respect to the x-axis), 
and 9 (the polar angle with respect to the z-axis). The 
polar direction is usually described by the pseudorapid- 
ity, defined as r\ = tanh - (cos#). 

In the previous section it was noted that the final state 
from tt decay may contain electrons, muons, jets, and 
neutrinos. The D0 detector was designed to identify 
and measure the energy/momentum of all of these ob- 
jects. As shown in Fig. 7, the detector has three major 
subsystems: the central tracking chambers, a uranium 
liquid-argon calorimeter, and a muon spectrometer. The 
detector design was optimized for high-resolution, nearly 
hermetic calorimetry that provides the sole measurement 
of the energies of electrons and jets. Due to the compact 
design of the calorimeter, the inner tracking volume is 
relatively small, and there is no central magnetic field. 

The central tracking detectors measure the trajectories 
of charged particles and aid in the identification of elec- 
trons. The former function is performed using three wire- 
chamber systems, and the latter by a transition-radiation 
detector (TRD). The three wire-chamber systems consist 
of two concentric cylindrical chambers centered on the 
interaction point and a set of two forward drift cham- 
bers that are situated at the ends of the cylinder. These 
chambers provide charged-particle tracking over the re- 
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gion \rj\ < 3.2, measuring the trajectories of charged par- 
ticles with a resolution of 2.5 mrad in (f> and 28 mrad 
in 6. The position of the interaction vertex along the 
beam direction (z) can be determined with a resolution 
of 8 mm. These chambers also measure the track ioniza- 
tion for distinguishing singly charged particles and e + e~ 
pairs from photon conversions. Concentric with, and ra- 
dially between, the two central chambers is the TRD. By 
measuring the amount of radiation emitted by single iso- 
lated particles as they pass through many thin sheets of 
polypropylene, this detector aids in the separation of elec- 
trons from charged pions and 7r ± /7 overlaps (since the 
amount of emitted transition radiation is proportional to 
the value of E/m for the particle). This device provides 
a factor of ten rejection of pions while retaining 90% of 
isolated electrons. 

Surrounding the central tracking system is the 
calorimeter, which is composed of plates of uranium and 
stainless steel/copper absorber surrounded by liquid ar- 
gon as the sensitive ionization medium. The calorime- 
ter is divided into three parts, the central calorimeter 
(CC), \r]\ < 1.2, and two end calorimeters (EC), which 
together cover the pseudorapidity range \rj\ < 4.2. Each 
consists of an inner electromagnetic (EM) section, a fine 
hadronic (FH) section, and a coarse hadronic (CH) sec- 
tion, housed in a steel cryostat. Each EM section is 21 ra- 
diation lengths deep and is divided into four longitudinal 
segments (layers) . The hadronic sections are 7-9 nuclear 
interaction lengths deep and are divided into four (CC) 
or five (EC) layers. The outer layer of each hadronic 
calorimeter is known as the "outer hadronic layer" . 
The calorimeter is transversely segmented into pseudo- 
projective towers with Ar] x Acj) = 0.1 x 0.1. The third 
layer of the EM calorimeter, in which the maximum of 
EM showers is expected, is segmented twice as finely into 
cells with Ar/xAcj) = 0.05x0.05. With this fine segmenta- 
tion, the azimuthal position resolution for electrons with 
energy above 50 GeV is ab out 2.5 m m. The energy reso- 
lution is a(E)/E = 15%/v/£(GeV)e0.4% for electrons. 

For charged pions the resolution is about 50%/ y/E(GcV) 
and for jets it is about 80%/y/E(GeV) [59]. For min- 
imum bias data, the resolution for each component of 
I/t, Ifix and Ifiy, has been measured to be 1.08 GeV + 
0.019 (Y i Et) 1 where EEr is the scalar sum of the trans- 
verse energies in all calorimeter cells. In order to improve 
the energy resolution for jets that straddle two cryostats, 
an inter-cryostat detector (ICD) made of scintillator tiles 
is situated in the space between the EC and CC cryostats. 
In addition, separate single-cell structures called "mass- 
less gaps" (MG) are installed in the intercryostat region 
in both the CC and EC calorimeters. 

The D0 muon detection systems cover \rj\ < 3.3. Since 
muons from top quark decays predominantly populate 
the central region, this work uses only the wide-angle 
muon spectrometer (WAMUS) which consists of four 
planes of proportional drift tubes (PDT) in front of mag- 
netized iron toroids with a magnetic field of 1.9 T and two 



groups of three planes each of proportional drift tubes be- 
hind the toroids. The magnetic field lines and the wires 
in the drift tubes are oriented transversely to the beam 
direction. The WAMUS covers the region \r]\ < 1.7 over 
the entire azimuth, with the exception of the central re- 
gion below the calorimeter (|rj| < 1, 225° < (j> < 315°), 
where the inner layer is missing to make room for the 
calorimeter support-structure. The WAMUS system is 
divided into the central iron (CF), |?7| < 1.0, and end 
iron (EF), 1.0 < \r/\ < 1.7, regions. As will be discussed 
in Sec. IV B, the EF region was used for only part of 
the Run 1 data set. The total thickness of the mate- 
rial in the calorimeter and iron toroids varies between 
13 and 19 interaction lengths, making background from 
hadronic punchthrough negligible. The tracking volume 
is small, thereby reducing backgrounds to prompt muons 
from in-flight decays of it and K mesons. The muon mo- 
mentum p is measured from its deflection angle in the 
magnetic field of the toroid. The momentum resolution 
is limited by multiple Coulomb scattering in the material 
traversed, the position resolution in the muon chambers, 
and uncertainty in the magnetic field integral. The typi- 
cal resolution in 1/p is approximately Gaussian and given 

by 

5{l/p) = 0.18(p- 2)/p 2 ®0.003 (2.1) 

(with p in GeV/c). 

As shown in Fig. 7, a separate synchrotron, the Main 
Ring, sits above the Tevatron and passes through the for- 
ward muon system and the outer hadronic section of the 
calorimeters. During data taking, it was used to acceler- 
ate protons for antiproton production. Losses from the 
Main Ring can deposit energy in the calorimeters and 
muon system, increasing the instrumental background. 
As discussed below (Sees. Ill, VI, and VII), these "Main- 
Ring events" are removed during the initial selection of 
all channels. Nevertheless, as discussed in Appendix C, 
and Sees. VI A and VII A, several analyses have been able 
to recover some, or all, of these events. 

III. TRIGGERS 

During normal operation, the Tevatron maintains two 
counter-rotating beams, one consisting of six bunches of 
protons and the other consisting of six bunches of an- 
tiprotons. Proton and antiproton bunches collide at the 
D0 interaction region every 3.5 /is (286 kHz). The D0 
trigger system is used to select the interesting events and 
reduce this to a rate of approximately 3-4 Hz, suitable 
for recording on tape. 

The D0 trigger system is composed of three hardware 
stages (level 0, level 1, and level 1.5) and one software 
stage (level 2) [59,58]. The first stage (level 0) consists 
of hodoscopes of scintillation counters mounted close to 
the beam on the inner surfaces of the end-calorimeter 



9 



cryostats and registers hits consistent with a pp interac- 
tion. This stage is typically used as an input to level 
1, but level is not required to fire before an event can 
proceed to the next stage. In addition, level is used to 
measure the luminosity. The next stage (level 1) forms 
fast analog sums of the transverse energies in calorimeter 
towers. These towers have a size of Arj x A<f> = 0.2 x 0.2 
and are segmented longitudinally into electromagnetic 
and hadronic sections. Based on these sums and pat- 
terns of hits in the muon spectrometer, the level 1 trigger 
decision takes place within the space of a single beam 
crossing, unless a level 1.5 decision is required (see be- 
low). Events accepted at level 1 are digitized and passed 
on to the level 2 trigger which consists of a farm of 48 
general-purpose processors. Software filters running on 
these processors make the final trigger decision. 

At both level 1 and level 2, the triggers are defined in 
terms of specific objects: electron/photon, muon, jet, $ T . 
Tables 1-4 show the triggers used for tt event selection. 
Table 5 shows the triggers used for the muon tag-rate 
studies discussed in Sec. VII B. As noted above, level 
is treated as an input term to level 1. Level 1 triggers 
that do not demand a level pass are denoted "NoLO." 

At level 1, the triggers for electrons (and photons) 
require the transverse energy in the EM section of the 
calorimeter to be above programmed thresholds: Et = 
Esin9 > T, where E is the energy deposited in the tower, 
9 its angle with the beam as viewed from the center of 
the detector (z = 0), and T a programmable threshold. 
The level 2 electron triggers exploit the full segmenta- 
tion of the EM calorimeter to identify electron showers. 
Using the trigger towers above threshold at level 1 as 
seeds, the algorithm forms clusters that include all cells 
in the four EM layers and the first FH layer in a region of 
At] x A<p = 0.3 x 0.3, centered on the highest Et tower. 
It checks the shower shape against criteria on the fraction 
of the energy found in the different EM layers. The Et 
of the electron is computed based on its energy and the 
z position of the interaction vertex as determined from 
the timing of hits in the level hodoscopes. The level 
2 algorithm can also apply an isolation requirement or 
demand an associated track in the central detector. 

During the later portion of the run, the level 1.5 trig- 
ger processor became available for selecting electrons and 
photons. For this purpose, the Et of each EM trigger 
tower passing the level 1 threshold is summed with the 
neighboring tower that has the most energy and a cut is 
made on this sum. The hadronic portions of the two tow- 
ers are also summed and the ratio of the EM transverse 
energy to the total transverse energy of the two towers 
is required to be greater than 0.85. The use of a level 
1.5 electron trigger is indicated in Tables 1-5 as an "EX" 
tower in the level 1 column. 

Muon triggers make use of hit patterns in the muon 
chambers at level 1 and provide the number of muon 
candidates in different regions of the muon spectrometer. 
The algorithm searches for hit patterns consistent with 
a muon originating from the nominal vertex (z = 0). 



A level 1.5 processor is also available and can be used to 
place a pt requirement on the candidates (at the expense 
of a slightly increased dead time). The use of a level 1.5 
muon trigger is indicated in Tables 1-5 as an "MX" muon 
in the level 1 column. 

At level 2, muon tracks are reconstructed using the 
muon PDT hits and the z position of the interaction ver- 
tex from level 0. Valid muon track selection is based on 
the muon pt and quality requirements (similar to those 
of Sec. IV B 1). The level 2 muon trigger can also require 
the presence of a minimum ionizing particle trace in the 
calorimeter cells along the track. This requirement is in- 
dicated in Tables 1-5 by "cal confirm." In addition, in 
between Run la and Run lb, layers of scintillator were 
added to the exterior of the central muon system to veto 
cosmic rays. The muon triggers indicated by "scint" re- 
quired the scintillator timing to be in a window of 30 ns 
before to 70 ns after the beam crossing. 

Jet triggers use projective towers of energy deposi- 
tion in the calorimeter similar to the EM trigger towers 
but including energy from the hadronic portion of the 
calorimeter. Level 1 jet triggers require the sum of the 
transverse energy in the EM and FH sections of a trig- 
ger tower (jet tower) to be above programmed thresh- 
olds: Esin9 > T, where E is the energy deposit in the 
tower, 9 its angle with the beam as seen from the center 
of the detector (z = 0), and T a programmable thresh- 
old. Alternatively, level 1 can sum the transverse energies 
within "large tiles" of size 0.8 x 1.6 in r\ x <p an d cut on 
these sums. The level 2 jet algorithm begins with an Et- 
ordered list of towers that are above threshold at level 
1. A level 2 jet is formed by placing a cone of radius 
AR = \J At] 2 + A<p 2 around the seed tower from level 1. 
If another seed tower lies within the jet cone, then it is 
passed over and not allowed to seed a new jet. Using 
the vertex position measured by the level hodoscopes, 
the summed Et in all of the towers included in the jet 
defines the jet Et- If any two jet cones overlap, then 
the towers in the overlap region are added into the jet 
candidate that was formed first. 

Jf^ , the missing transverse energy as measured in the 
calorimeter (see Sec. IV D for definition), can be com- 
puted at both level 1 and level 2. At level 1, the z position 
is assumed to be z — 0. At level 2, the vertex position 
from level is used. In the offline reconstruction, the 
determination of uses the z position as determined 
by the tracking system. Therefore, the resolution of If^ 1 
at the trigger level is significantly poorer than that in the 
offline reconstruction. 

As noted in Sec. II, the Main Ring passes directly 
through a portion of the outer hadronic calorimeter and 
muon system. Particles lost from the Main Ring can 
affect the measurements in these subsystems. Several 
schemes were employed at the trigger level to reduce 
or eliminate these effects; these are described in Ap- 
pendix B. 

In addition to the complications introduced by the 
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TABLE 1. Electron triggers used in collection of the tt signal sample. Column 1 gives the trigger name, column 2 gives the 
run period for which it was applied, column 3 gives the exposure in pb _1 (see text for definition), columns 4 and 5 give the 
level 1 and level 2 definitions, and column 6 lists the channels that used each trigger. See Appendix C for definitions of the 
MR veto terms: GB, mrbs, ml, and GC. Channel names are defined in Sees. VI and VII. 



Name 


Run 


Expsr. 

(Pb- 1 ) 


Level 1 


Level 2 


Used by 


ELE-HIGH 


la 


11.0 


1 EM tower, E T > 10 GeV 
GB 


1 isolated e, E T > 20 GeV 


e + jcts/topo 


ELE-.IET 


la 


14.4 


1 EM tower, K T > 10 GeV, |rj| < 2.6 
2 jet towers, E T > 5 GeV 
MRBS 


1 e, E T > 15 GeV, |r/| < 2.5 
2 jets (Ail = 0.3), E T > 10 GeV, \r)\ < 2.5 
> 10 GeV 


ee, e/^, ei^ 
e + jets 
e + jets/ fj. 


ELE- JET-HIGH 


lb 


98.0 


1 EM tower, E T > 12 GeV, \-q\ < 2.6 

2 jet towers, E T > 5 GeV, \r)\ < 2.0 

ML 


1 e, E T > 15 GeV, |r/| < 2.5 
2 jets (Ail = 0.3), E T > 10 GeV, \r)\ < 2.5 
#p' > 14 GeV 


ee, e/^, ei/ 
e + jcts/topo 
e + jets/ju 


ELE- JET-HIGH 


lc 


1.9 


1 EM tower, E T > 12 GeV, \r)\ < 2.6 

2 jet towers, E T > 5 GeV, \r)\ < 2.0 

ML 


1 e, E T > 15 GeV, \ V \ < 2.5 
2 jets (Ail = 0.3), E T > 10 GeV, |r)| < 2.5 
$5?' > 14 GeV 


ee, e^i, ei/ 
e + jets/// 


ELE- JET- HIGH A 


lc 


11.0 


1 EM tower, E T > 12 GeV, \r)\ < 2.6 

2 jet towers, E T > 5 GeV, \r)\ < 2.0 

1 EX tower, E T > 15 GeV 
ML 


1 e, E T > 17 GeV, |r/| < 2.5 
2 jets (Ail = 0.3), E T > 10 GeV, |r;| < 2.5 
^jf 1 > 14 GeV 


ee, e/j-, ei/ 
e + jets//x 


EMl-EISTRKCC-MS 


lb 


93.4 


1 EM tower, E T > 10 GeV 
1 EX tower, E T > 15 GeV 
GC, NoLO 


1 isolated e w/track, E T > 20 GeV 
^r al > 15 GeV 


e^ 

e + jcts/topo 


MU-ELE 


la 


13.7 


1 EM tower, E T > 7 GeV 
1 P, \V\ < 2-4 
MRBS 


1 e, E T > 7 GeV 
1 p, p T > 5 GeV/c, |r;| < 2.4 


e/i 




lb 


93.9 


1 EM tower, E T > 7 GeV 
1 MX p, \ti\ < 2.4 
GC 


1 e, E T > 7 GeV, < 2.5 
1 p, p T > 8 GeV/c, < 2.4 


e/i 


MU-ELE-HIGH 


lc 


10.6 


1 EM tower, E T > 10 GeV, |rj| < 2.5 
1 MX p, |r/| < 2.4 
GC 


1 e, E T > 10 GeV, |r/| < 2.5 
1 n, p T > 8 GeV/c, |t7| < 1.7 





TABLE 2. Muon+jet triggers used in collection of the tt signal sample. Column I gives the trigger name, column 2 gives 
the run period for which it was applied, column 3 gives the exposure in pb -1 (see text for definition), columns 4 and 5 give 
the level 1 and level 2 definitions, and column 6 lists the channels that used each trigger. See Appendix C for definitions of the 
MR veto terms: GB and GC. Channel names are defined in Sees. VI and VII. 



Name 


Run 


Expsr. 

(Pb" 1 ) 


Level 1 


Level 2 


Used by 


MU- JET-HIGH 


la 


10.2 


1 M, M < 2.4 
1 jet tower, E T > 5 GeV 
GB 


1 ju, p T > 8 GeV/c, |7/| < 1.7 
1 jet (Ail = 0.7), E T > 15 GeV 


ep, pfi 

p + jcts/topo 
p + jets/ju 


lb 


66.4 


1 ft, pt > 7 GeV/c, |r/| < 1.7 
1 jet tower, E T > 5 GeV, |r)| < 2.0 
GC 


1 fi, pt > 10 GeV/c, \r)\ < 1.7, scint 
1 jet (Ail = 0.7), E T > 15 GeV, |r,| < 2.5 


ep, pfi 

p + jets/topo 
p + jets/p 


MU-JET-CAL 


lb 


88.0 


1 fi, p T > 7 GeV/c, |77| < 1.7 
1 jet tower, E T > 5 GeV, |r)| < 2.0 
GC 


1 p, p T > 10 GeV/c, |tj| < 1.7 
cal confirm, scint 
1 jet (Ail = 0.7), E T > 15 GeV, |r,| < 2.5 


pp 

p + jets/topo 
p + jets/ju 


MU- JET-CENT 


lb 


48.5 


1 M, M < i.o 

1 jet tower, E T > 5 GeV, |r)| < 2.0 
GC 


1 n, p T > 10 GeV/c, \r)\ < 1.0, scint 
1 jet (Ail = 0.7), E T > 15 GeV, |r;| < 2.5 


e/i, ju/i 

/i + jets/topo 
p + jets/ju 


lc 


8.9 


1 H, |rj| < 1.0 
1 jet tower, E T > 5 GeV, r; < 2.0 
2 jet towers, E T > 3 GeV 
GC 


1 n, pt > 12 GeV/c, |7j| < 1.0, scint 
1 jet (Ail = 0.7), E T > 15 GeV, < 2.5 


e/i, /i/t 


MU-JET-CENCAL 


lb 


51.2 


i n, \n\ < i o 

1 jet tower, E T > 5 GeV, |r)| < 2.0 
GC 


1 n, p T > 10 GeV/c, |>j| < 1.0 
cal confirm, scint 
1 jet (Ail = 0.7), E T > 15 GeV, < 2.5 


pp 

p + jets/topo 
/i + jets/n 


lc 


11.4 


1 H, M < 1.0 
1 jet tower, E T > 5 GeV, \r)\ < 2.0 
2 jet towers, E T > 3 GeV 
GC 


1 fj,, p T > 12 GeV/c, < 1.0 
cal confirm, scint 
1 jet (Ail = 0.7), E T > 15 GeV, < 2.5 


e/i, /i/i 
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TABLE 3. Jet triggers used in collection the ti signal sample. Column 1 gives the trigger name, column 2 gives the run 
period for which it was applied, column 3 gives the exposure in ph" 1 (see text for definition), columns 4 and 5 give the level 1 
and level 2 definitions, and column 6 lists the channels that used each trigger. See Appendix C for definitions of the MR veto 



terms: ml, MB, and mrbs. The lepton+jets channels are defined in Sec. VII. 



Name 


Run 


Expsr. 

(Pb- 1 ) 


Level 1 


Level 2 


Used by 


JET-3-MU 


1 K 

lb 


11.9 


o jet towers, C/T -> o Ogv 
#5?' > 20 GeV 
ML 


o jets [i\K — v.i), Et > lo UcV, \rj\ < Z.o 
^r al > 17 GeV 


+ jcts/topo 
H + jets/ /i 


JET- 3- MISS-LOW 


ib 


57.8 


3 large tiles, E T > 15, \r)\ < 2.4 
3 jet towers, E T > 7 GeV, |r)| < 2.6 
MB 


3 jets (AH = 0.5), E T > 15 GeV, \r)\ < 2.5 
> 17 GeV 


+ jcts/topo 
+ jcts//j 


jet-3-l2mu 


ib 


25.8 


3 large tiles, E T > 15, \r)\ < 2.4 
3 jet towers, E T > 7 GeV, |r)| < 2.6 
MB 


1 fi, p T > 6 GeV/c, < 1.7 
cal confirm, scint 
3 jets (AB, = 0.5), E T > 15 GeV, |r/| < 2.5 

^cal > 1? GcV 


+ jcts/topo 
fj. + jcts//j 


JET-MULTI 


la 


14.6 


4 jet towers, E T > 5 GeV 
MRBS 


5 jets (AR = 0.3), E T > 10 GeV, \r)\ < 2.0 


all-jets 


lb 


96.6 


3 large tiles, E T > 15 GeV, |>j| < 2.4 
3 jet towers, E T > 7 GeV, |r;| < 2.6 
and 1 jet tower, E T > 3 GeV 
ML 


5 jets (AB. = 0.3), £ T > 10 GeV, \ v \ < 2.5 
T,E T > 100 GeV for jets with \rj\ < 2.5 


all-jets 


lc 


11.3 


3 large tiles, E T > 15 GeV, \r/\ < 2.4 
3 jet towers, E T > 7 GeV, |r)| < 2.6 
and 1 jet tower, E T > 3 GeV 
ML 


5 jets (AH = 0.3), E T > 10 GeV, |r)| < 2.5 
T,E T > 120 GeV for jets with |r;| < 2.5 


all-jets 



TABLE 4. $ T triggers used in collection of the ti signal sample. Column 1 gives the trigger name, column 2 gives the run 
period for which it was applied, column 3 gives the exposure in pb" 1 (see text for definition), columns 4 and 5 give the level 1 
and level 2 definitions, and column 6 notes that these triggers were used only by the ev channel. See Appendix C for definitions 



of the MR veto terms: mrbs and GB. The ev channel is defined in Sec. VI. 



Name 


Run 


Expsr. 

(Pb" 1 ) 


Level 1 


Level 2 


Used by 


MISSING-ET 


la 


13.7 


t^' > 30 GeV 
1 jet tower, E T > 5 GeV, \-q\ < 2.6 
MRBS 


#5?' > 35 GeV 


ev 




lb 


83.6 


#5? J > 40 GeV 
1 jet tower, E T > 5 GeV, |r;| < 2.6 
GB 


#5?' > 40 GeV 


ev 


MISSING-ET-HIGH 


lc 


0.7 


#5?' > 50 GeV 
1 jet tower, E T > 5 GeV, |r;| < 2.6 
GB 


#5?' > 50 GeV 


ev 
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TABLE 5. Triggers used to study the £+jets//n backgrounds and tag rate function (see Sec. VII B). Column 1 gives the 
trigger name, column 2 gives the run period for which it was applied, column 3 gives the exposure in pb _1 (see text for 
definition), columns 4 and 5 give the level 1 and level 2 definitions, and column 6 notes that these triggers were used only for 
£+jets background studies. See Appendix C for definitions of the MR veto terms: GB, MRBS, ml, and GC. The lepton+jets 
channels are defined in Sec. VII. 



Name 


Run 


Expsr. 
(Pb- 1 ) 


Level 1 


Level 2 


Used by 


JET-MIN 


lb 


0.007 


1 jet tower, E T > 3 GeV 
GB 


1 jet (AR - 0.3), E T > 20 GeV 
prescale — 20 


£ + jets/ fi 
bkg 


JET-3-MON 


lb 


0.92 


2 jet towers, E T > 5 GeV 
and 1 jet tower, Et > 3 GeV 
GB 


3 jets (AR = 0.3), E T > 10 GeV 
prescale — 5 


£ + jets//j 
bkg 


JET-4-MON 


lb 


4.6 


2 jet towers, E T > 5 GeV 
and 1 jet tower, Ex > 3 GeV 
GB 


4 jets (AR = 0.3), E T > 10 GeV 


£ + jets/ /i 
bkg 


JET-MULTI 


la 


14.6 


4 jet towers, E T > 5 GeV 
MRBS 


5 jets (AR = 0.3), E T > 10 GeV, |tj| < 2.0 


£ + jets//j 
bkg 


lb 


96.6 


3 large tiles, E T > 15 GeV, \r)\ < 2.4 
3 jet towers, E T > 7 GeV, |r;| < 2.6 
and 1 jet tower, Et > 3 GeV 
ML 


5 jets (AR. = 0.3), E T > 10 GeV, |r)| < 2.5 
SE T > 100 GeV for jets with 1 77 j < 2.5 


£ + jets//j 
bkg 


ELE-1-MON 


lb 


3.1 


1 EM tower, E T > 7 GeV, |r;| < 2.5 
1 jet tower, E T > 3 GeV 
GC 


1 e, E T > 16 GeV/c 


£ + jets/ /j 
bkg 


GIS-DIJET 


lb 


93.5 


1 EM tower, E T > 10 GeV, < 2.5 
1 jet tower, E T > 3 GeV 
GC 


1 isolated e/7, E T > 15 GeV/c, \-q\ < 2.0 
3 jets (AR = 0.7), E T > 15 GeV, |rj| < 2.0 
T,E T > 70 GeV for jets with |r;| < 2.0 


£ + jets/^ 
bkg 


EMl-EISTRKCC-ESC 


lb 


91.9 


1 EM tower, E T > 10 GeV, |r f | < 2.5 
1 jet tower, E T > 3 GeV 
GC 


1 e (no shape cuts), Et > 16 GeV 
and 1 isolated e w/track, Et > 20 GeV 


£ + jets//j 
bkg 



Main Ring, there are also effects due to multiple in- 
teractions. At the mean luminosity (7.5 x 10 30 /cm 2 /s), 
there are on average 1.3 interactions per bunch crossing. 
Since the cross section for the production of high-py in- 
teractions is small compared to that for minimum bias, 
it is very unlikely that more than one high-p^ interac- 
tion will be present in any given bunch crossing. These 
additional minimum-bias interactions are usually not in- 
cluded in the Monte Carlo models, but do contribute to 
mismeasurement of the primary interaction vertex, and 
therefore to mismeasurement of lepton and jet transverse 
energies /moment a. The systematic uncertainty due to 
multiple interactions is discussed in Sec. IX A 7. 

The Run 1 data were acquired in three separate run pe- 
riods: Run la from 1992-1993, Run lb from 1994-1995, 
and Run lc from 1995-1996. The period appropriate to 
each trigger is given in the second column of Tables 1-5. 

The integrated luminosity L was determined from the 
counting rate in the level hodoscopes (Rlo) as 



C 



-ln(l - tR L0 ) 
ra L0 



(3.1) 



where r = 3.5 fis is the time interval between beam cross- 
ings and ctlo is the effective pp cross section subtended by 
the level counters. As described in detail in Ref. [60], 
clo = 43.1 ± 1.9 mb is obtained from the level trigger 
efficiency and geometrical acceptance, and from a "world 
average" pp total inelastic cross section of 57.39 ± 1.56 
mb based on results from the CDF [61], E710 [62], and 
E811 [63] collaborations at Fermilab. The level trigger 
efficiency is determined using samples of data collected 
from triggers on random beam crossings and the geomet- 
rical acceptance from Monte Carlo studies. It should be 



noted that the CDF luminosity determinations are based 
solely on its own measurement of the pp inelastic cross 
section. As a result, luminosities reported by CDF are 
6.2% lower than those currently reported by D0, and 
consequently, all CDF cross sections are ab initio 6.2% 
larger than all D0 cross sections. Earlier D0 cross sec- 
tions (and all previous D0 tt cross sections) were based 
on a pp inelastic cross section determined only from the 
CDF and E710 measurements and are 3.2% lower than 
current D0 cross sections. 

The integrated luminosity (exposure) seen by each of 
the triggers is given in the third column, labelled "Ex- 
psr.," of Tables 1-5. These values include luminosity 
losses due to Main-Ring vetos and prescale factors (if 
appropriate), but do not include the loss to the offline 
GOOD-BEAM requirement or losses from runs rejected at 
later stages of the analysis (see Appendix B for a discus- 
sion of the Main- Ring veto schemes). 



IV. EVENT RECONSTRUCTION 

A. Electron identification 

Electrons and positrons arc identified by the distinc- 
tive pattern of energy that electromagnetic showers de- 
posit in the calorimeter and by the presence of a track 
from the interaction vertex to the cluster of hit calorime- 
ter cells. The algorithm for clustering calorimeter energy 
and quantities used to distinguish electrons from back- 
grounds are described in Ref. [58]. The present analysis 
includes two additional features: the separation between 
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electrons and backgrounds has been improved by the in- 
troduction of a multivariate discriminant, and, for the 
dilepton channels, use is made of information from the 
TRD. 

The electromagnetic energy scale was calibrated using 
Z — > ee, J/ip — + ee, and 7r° — > 77 decays to a precision of 
0.08% at E = M z /2 and to 0.6% at E = 20 GeV [64,65]. 

The complete set of identification variables, efficien- 
cies, and misidentification rates is discussed below. Un- 
less otherwise indicated, electrons specified to be in the 
CC region of the detector span the range < \r]\ < 1.2 
and electrons specified to be in the EC region of the de- 
tector span the range 1.2 < \rj\ < 2.0 (with the region 
between the cryostats, 1.2 < \r]\ < 1.5, having only a 
minimal acceptance). Since the central tracking system 
does not measure the charge of particles, it is not possible 
to distinguish between electrons and positrons. There- 
fore, for the remainder of this paper, "electron" shall be 
used to indicate both electrons and positrons. 



layer lying in a six by six array around the tower contain- 
ing the highest energy cell; the logarithm of the cluster 
energy; and the z position of the interaction vertex. The 
elements of the covariance matrix depend on rj and were 
determined using the D0GEANT [66] model of the detec- 
tor (see Sec. V). 

4- Cluster-track match significance (auk) 

Calorimeter clusters are required to lie along the tra- 
jectories of charged particle tracks reconstructed in one 
of the inner tracking chambers. The cluster-track match 
significance 0- tr k is a measure of the distance between the 
cluster centroid and the intersection of the extrapolated 
track to the third layer of the EM calorimeter. 

5. Track ionization (dE/dx) 



1. Electromagnetic energy fraction 

Electromagnetic energy clusters are formed by com- 
bining calorimeter towers using a nearest-neighbor algo- 
rithm with EM tower seeds. The electromagnetic energy 
fraction /em of a cluster is the ratio of its energy found 
in EM calorimeter cells to its total energy. All electron 
candidates are required to have /em > 0.9. 



2. Isolation fraction (I) 

Electron showers are compact and mostly contained in 
the core of EM cells within a radius R = 0.2 in {r},(j)) 
around the shower center. The isolation fraction I is 
defined as the ratio of energy in non-core EM and FH 
cells (-Etot) within a cone of 0.4 around the center to 
energy in the EM cluster core (-Eem) 



ff tot (0.4)-g E M(0.2) 
£em(0.2) 



(4.1) 



This quantity tends to be substantially smaller for elec- 
trons from the decay of W and Z bosons than for the 
background, most of which originates from hadronic jets 
where the electron candidate is usually accompanied by 
nearby energetic particles. 



Photons that convert to e + e~ pairs before the 
calorimeter produce pairs of tracks that match an EM 
cluster well and are too close together to be resolved. 
Such double tracks can be identified by the amount of ion- 
ization produced along the track (dE/dx); photon con- 
versions typically deposit twice the charge expected from 
one minimum ionizing particle. 



6. TRD efficiency (e t ) 

The response of the TRD is characterized by the vari- 
able e t : 



et (AE) 



ir E m(E)dE 

J °°™(E)dE> 



(4.2) 



where AE is the difference between the total energy 
recorded in the TRD (E) and that recorded in the layer 
with the largest signal (this is done to reduce sensitivity 
to (5-rays) and ^ is the electron energy spectrum from 
a sample of W — > ev events [67,68]. Hadrons generally 
deposit energy mainly in a single layer (giving a small 
value for AE) and electrons deposit energy more evenly 
(giving a larger value for AE). Therefore, hadrons tend 
to have values of et near unity whereas the distribution 
from electrons is roughly uniform over the allowed range 
from to 1. 



3. Covariance matrix (xt) 

A covariance matrix is used to compute a x 2 variable 
(Xe) representing the consistency of the cluster shape 
with that of an electron shower. The covariance matrix 
uses 41 variables: the fractions of energy deposited in 
the first, second, and fourth layers of the EM calorime- 
ter; the fractions of energy in each cell of the third EM 



7. Likelihood ratio (L 4 ,L 5 ) 

In order to attain the maximum background rejection 
while keeping a high efficiency for real electrons, the 
variables /em, Xe, CT trk, and dE/dx are combined into 
an approximate four-variable likelihood ratio L 4 for the 
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FIG. 8. (a) - (e) Electron identification variables used in 
the Z/4 and L5 likelihood ratios, (f) Isolation, and (g) and 
(h) 4-variable and 5-variable likelihood ratios. The open his- 
tograms are from electron candidates from Z — > e + e~ events 
and the shaded histograms are from electron candidates from 
EM clusters in inclusive jet data (mainly background). Ar- 
rows indicate the position of the cuts on isolation, L4, and 
is . All quantities are for the CC region of the detector only. 



The distributions associated with all the above vari- 
ables for electrons in the CC region of the detector are 
shown in Fig. 8. 



8. Selection 

Based on these quantities, four classes of electron can- 
didates are denned: 

1. extra-loose electrons are denned as objects satisfy- 
ing Jem > 0.9, X< 0.3, and xl < 300. 

2. minimal electrons are defined as objects satisfying 
Jem > 0.9 and I< 0.1. 

3. loose electrons are defined as the subset of the 
extra-loose sample that satisfies the additional re- 
quirements X< 0.1 and L$ < 0.5 for CC and EC 
clusters. 

4. tight electrons are denned as the subset of the 
extra-loose sample that satisfies the additional re- 
quirements X< 0.1 and L 4 < 0.25(0.3) for CC (EC) 
clusters. 

The loose definition is used for the final selection in 
the dilepton channels (ee, e/j,, ev). The tight definition 
is used for the final selection in the e+jets channels. 



9. Efficiency 



hypotheses that a candidate electron is signal or back- 
ground. Similarly, the variables /em, Xe, 0trk) dE/dx, 
and Ct are combined into an approximate five-variable 
likelihood ratio L5. These likelihood ratios are defined 
using the Neyman-Pearson test for signal (e) and back- 
ground (b) hypotheses, where an EM cluster is considered 
to be an electron if it satisfies 



Pn(x\b) 
p n (x\e) 



< k, 



(4.3) 



where x is the vector of observables, p n (x\H) is the prob- 
ability density for x if the hypothesis H is true, and k 
is the cutoff value. The probability densities arc com- 
puted by forming the joint likelihood of the four or five 
variables: 

p 4 (x\H)=p(f EM \H)-p( X 2 e \H)-p(a tlk \H) 

■p(dE/dx\H), (4.4) 

p 5 (x\H)=p(f EM \H)-p( X 2 e \H)-p(a tTk \H) 

■p{dE/dx\H)-p{e t \H), (4.5) 

where p(y\H) is the probability density for a single vari- 
able y if the hypothesis H is true. These signal and 
background hypotheses are constructed respectively from 
inclusive Z — ► e+e~ data and inclusive jet production. 



The efficiencies for electron identification are obtained 
by using the Z — > ee mass peak. The procedure is based 
on a sample of events from the emI-eistrkcc-esc trig- 
ger (see Table 5) that has two reconstructed electromag- 
netic clusters, each with Et > 20 GeV. From this sample, 
one of the electron candidates, denoted as the "tag," is re- 
quired to be a good electron (xl — 100, X < 0.15). If the 
other electromagnetic cluster, denoted as the "probe," 
satisfies X < 0.1, then the invariant mass of the pair, 
m(tag, probe), is recorded. This is done separately for 
probes in the CC and EC regions of the calorimeter. 
The number of entries in the Z boson mass window, 
80 GeV/c 2 < m(tag, probe) < 100 GeV/c 2 , with back- 
ground subtracted, and in the instrumented region of the 
central tracking system, defines the number of true elec- 
tron probes [69]. The track finding efficiency e t rk is de- 
fined as the ratio of the number of true electron probes 
with a track to the total number of true electron probes. 
This efficiency varies with the number of interactions 
per event (see Sees. Ill and IX A 7). Typical values are 
82.7 ± 1.1 % for electrons in the CC and 85.2 ± 1.0 % 
in the EC. The electron identification efficiencies, given 
in Table 6, are defined by the ratio of the number of 
true electron probes with a reconstructed track that pass 
the given identification requirements to the total num- 
ber of true electron probes with a reconstructed track. 
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These efficiencies do not include geometric factors due 
to uninstrumented fiducial regions of detector. The ge- 
ometrical acceptance for electrons in the D0 detector is 
(87.6 ± 0.5)% in the CC and (79.2 ± 1.4)% in the EC. 

TABLE 6. Definition of loose and tight electron identifi- 
cation criteria and the corresponding efficiencies (Eff) and 
misidentification rates (i? m i s ). 





Loose 


Tight 


Region 


cc 


EC 


CC 


EC 


Def 


L 5 < .5 


L 5 < .5 


L 4 < .25 


L 4 < .3 


Eff(%) 


88.0 ± 1.6 


63.8 ±2.3 


81.1 ± 1.0 


51.4 ± 1.8 


Rmis(%) 


4.6 ±0.1 


8.0 ±0.1 


2.2 ±0.1 


4.0 ±0.3 



10. Misidentification rate (R m i s ) 

The electron misidentification rates (R m is) given in Ta- 
ble 6 are measured from a sample of QCD multijet events 
that contained one electromagnetic cluster passing the 
extra-loose electron identification requirements defined 
above. From this sample of extra-loose electron candi- 
dates, the fraction passing the loose/tight electron iden- 
tification is obtained separately for the CC and EC re- 
gions of the calorimeter and defined to be the rate for an 
extra-loose electron candidate to pass the loose/tight cri- 
teria. Note that the multijet backgrounds due to electron 
misidentification arc handled differently in the e±jets 
analyses and are discussed in Sees. VII A and VII B. 

B. Muon identification 

Muon tracks are reconstructed using the muon system 
PDTs. Additional information about the interaction ver- 
tex, matching tracks in the central tracker, and minimum 
ionizing traces left in the calorimeter is also available. 

As noted in Sec. II, the decay products from the ti 
pair are emitted at central rapidities and the muon iden- 
tification is therefore restricted to the central (WAMUS) 
portion of the D0 muon system, |ry| < 1.7. Due to ineffi- 
ciencies caused by radiation damage, the forward muon 
region (EF) with 1.0 < |?7| < 1.7 was not used in these 
analyses for Run la (w 10 pb _1 ) or the early part of 
Run lb (s=y 49 pb" 1 ). The chambers were subsequently 
cleaned and returned to full efficiency for the remainder 
of Run lb and Run lc. In the discussion below, the pre- 
cleaning period of Run lb is denoted as "predean" and 
the post-cleaning period as "postclean" . 

Several categories of muons are used in the analyses. 
The primary categories correspond to the selection of iso- 
lated muons arising dominantly from W — * fiv decay 
and non-isolated (tag) muons from b — > fi + X decays. 
Isolation implies a separation of the muon track from 



nearby jet activity. Isolated muons fall into two cate- 
gories, tight and loose. The selection requirements for the 
three types vary slightly over time and arc summarized 
in Tables 7 - 9 for Run la, Run lb(preclean), and Run 
lb±c(postclean) respectively. Tight and loose muons 
share most requirements except that tight muons have 
the additional requirements of an impact parameter cut 
and a minimum magnetic field path length (see below). 
The pt and Ai?(/i,jet) requirements for isolated muons 
are characteristic of what is expected from W — > fiv de- 
cay. The selection requirements for tag muons are very 
similar to those for loose-isolated muons except for the 
lower momentum threshold of pr > 4 GeV/c and the 
non-isolation requirement of Ai?(/x,jet) < 0.5. These pr 
and AR requirements select muons characteristic of those 
expected from heavy-flavor decays. 

The momentum of the muon is computed from the de- 
flection of the muon trajectory in the magnetized toroid. 
The momentum calculation uses a least squares method 
that considers seven parameters: four describing the po- 
sition and angle of the track before the calorimeter (in 
both the bend and non-bend views), two describing the 
effects due to multiple scattering, and the inverse of the 
muon momentum I /p. This seven-parameter fit is ap- 
plied to sixteen data points: vertex position measure- 
ments along the x and y directions, the angles and posi- 
tions of track segments before and after the calorimeter 
and outside of the iron, and two angles (one in the bend 
view, one in the non-bend view) representing the deflec- 
tion due to multiple Coulomb scattering of the muon in 
the calorimeter. Energy loss corrections are applied using 
the restricted energy loss formula parametrized in GEANT 

The muon momentum resolution depends on the 
amount of material traversed, the magnetic field inte- 
gral, and the precision of the measurement of the muon 
bend angle. As noted in Sec. II, the resolution function 
shown in Eq. 2.1, was based on studies of Z — > /ifi data. 
The first term in the resolution function reflects mul- 
tiple Coulomb scattering in the iron toroids and is the 
dominant effect for low momentum muons. The second 
term reflects the resolution of the muon position mea- 
surements. The muon momentum scale was calibrated 
using J/ip — > fin and Z — > [in. candidates and has an 
uncertainty of 2.5 %. 

The complete set of identification variables and 
misidentification rates is discussed below. 



1. Muon quality (Q) 

For each found track in the muon system, the recon- 
struction makes a set of cuts on the number of modules 
hit, the impact parameters, and the hit residuals. The 
number of cuts which a track fails is defined as the muon 
quality, Q (for a perfect track Q — 0). A similar param- 
eter is also produced by the level 2 trigger. If a track 
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fails more than one (CF) or any (EF) of the cuts on the 
above quantities, then it is of insufficient quality and is 
rejected. Tracks that have hits only in the inner layer 
of the muon system (inside the toroid) are also rejected. 
This eliminates almost all hadronic punchthrough from 
the calorimeter into the muon system. If a muon track 
is not bent by the toroid, muon momentum cannot be 
measured (as is the case for tracks which only have hits 
in the inner layers). 

2. Calmip/MTC requirement 

As a muon passes through the calorimeter it deposits 
energy through ionization along its path. This minimum 
ionizing trace should match to the track found in the 
muon and central tracking systems and can serve as a 
very powerful tool for the rejection of backgrounds. Dur- 
ing the course of the run this was used in two ways. For 
Run la, it is accomplished by checking the energy in the 
calorimeter towers along the expected path of the muon: 
For events in which a track match is found in the cen- 
tral tracking system within Arj < 0.45 and Acfr < 0.45 
of the muon track, an energy deposit of at least 0.5 GeV 
is required in the calorimeter towers along the track plus 
its two nearest neighbor towers; for muons without a cen- 
tral detector track match, at least 1.5 GeV is required (to 
allow for tracking inefficiencies in the region near \r/\ s=a 
1 where the coverage of the central tracking system is 
incomplete). This requirement is denoted by "calmip" 
in Tables 7-9. For data from Runs lb and lc, a more 
sophisticated procedure is employed. This procedure, de- 
noted "MTC," is based on muon tracking in the calorime- 
ter. The track from the muon system is used to define 
a path through the calorimeter to the position of the 
interaction vertex. A 5 x 5 road of calorimeter cells is 
defined along this path. Any cell with an energy two 
standard deviations above the noise level is counted as 
hit. The longest contiguous set of hit cells constitutes 
the calorimeter track. Muon candidates are required to 
have tracks with hits in at least 70% of the possible layers 
in the hadronic calorimeter. If a track does not have hits 
in all the layers, then it is also required that at least one 
of the nine central cells in the outermost layer of the 5x5 
road be hit [69]. These requirements reject both combi- 
natoric background and cosmic rays. The MTC criteria 
cannot be used on the Run la data because the required 
information is not supplied by the la reconstruction. For 
the /z+jets channels (which uses the tight muon identi- 
fication criteria) the Run la raw data was reprocessed, 
incorporating the needed information. Thus, in Table 7, 
MTC refers to the tight identification and the tag iden- 
tification for the fx + jets channels and calmip refers to 
the loose identification and the tag identification for the 
e+jets/^i channel. 



TABLE 7. Definitions of and identification efficiencies for 
loose, tight, and tag CF (\n\ < 1.0) muons for Run la. For 
calmip/MTC: e/j,, fi/i (loose) and e + jets//i (tag) use calmip; 
fi + jets/topo (tight) and /i + jets//! (tight and tag) repro- 
cessed the la data and therefore use MTC. The two efficien- 
cies given for tag muons reflect inclusion of calmip or MTC 
requirements respectively. 







fi id Run la (CF) 


definition: 


Loose 


Tight 


Tag 




15 


20 


4 


Q < 


1 


1 


1 


calmip/MTC 


yes 


yes 


yes 


IP < 




20 cm 




J Bdl > 




1.83 Tm 




AR(n,jet) 


> 0.5 


> 0.5 


< 0.5 


Eff (%) 


64 ±6 


46 ± 7 


80 ± 6/77 ± 6 



3. Impact parameter (IP) 

An impact parameter requirement for the muon tra- 
jectory relative to the interaction vertex provides fur- 
ther rejection against cosmic rays and misreconstructed 
muons. Here IPbv and IPnb are the two-dimensional 
distances-of-closest approach between the muon and its 
associated vertex in the bend and non-bend projections 
re spectively. Th ese are combined in quadrature, IP = 
\fPP~Bv + IPnb-> an d IP i s required to be less than 20 
cm. 

4- Minimum magnetic path length (J Bdl) 

Muons that pass through the thinner part of the iron 
toroid near \rj\ w 0.9 have poorer momentum resolution 
and may also be contaminated by a small background 
from punchthrough. Excluding these thin regions, the 
punchthrough fraction is < 2 % and is essentially negligi- 
ble for muons with p T > 5 GeV/c. The J Bxdl require- 
ment ensures that muons traverse enough field (> 1.83 
Tm) to provide an acceptable pt measurement. 

5. Isolation 

A muon is considered isolated if it is well separated 
from any reconstructed jet. Isolation, or Aii(/i,jet), is 
the distance in (77, <j>) space between a muon and the near- 
est 0.5 cone jet with E T > 8 GeV. 

6. Efficiency 

The total muon-finding efficiency is the product of 
the muon geometrical acceptance and the muon identi- 
fication efficiency. The muon geometrical acceptance is 
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TABLE 8. Definitions of and identification efficiencies 
for loose, tight, and tag CF (\rj\ < 1.0) muons for Run 
fb(preclean). 









\i id Run lb preclean 


(CF) 


definition: 




Loose 


Tight 




Tag 


P T > 




15 


20 




4 


Q < 




1 


1 




1 


MTC 




yes 


yes 




yes 


IP < 






20 cm 






J Bdl > 






1.83 Tm 






AR(jj,, jet) 




> 0.5 


> 0.5 




< 0.5 


Eff (%) 


65 ±5 


46 ±7 


76 ±6 



TABLE 9. Definitions of and identification efficiencies for 
loose, tight, and tag muons for CF {\r)\ < 1.0) and EF 
(1.0 < < 1.7) regions for run lb+c (postclean). 





ji id Run lb+c postclean 




Loose 


Tight 


Tag 


definition: 


CF EF 


CF EF 


CF EF 


P T > 


15 


20 


4 


Q < 


1 


1 


1 


MTC 


yes 


yes 


yes 


IP< 




20 cm 




J Bdl > 




1.83 Tm 




AR(n,jet) 


> 0.5 


> 0.5 


< 0.5 


Eff (%) 


73 ±3 68 ±5 


49 ± 7 52 ± 16 


84 ±4 62 ±15 



(73.7±0.4)% for the CF and (64.1±1.1)% for the EF. The 
total muon-finding efficiency is well-modeled by a mod- 
ified version of D0GEANT. These modifications include 
input from measured muon resolutions and efficiencies of 
the PDTs. The muon identification efficiency is obtained 
from this modified version of D0GEANT, but is further 
corrected to account for time dependent detector ineffi- 
ciencies and incorrect modeling of the muon track finding 
efficiency. As can be seen in Tables 7-9, the muon iden- 
tification efficiency varies with muon category and run 
period. 

C. Jets 

Jets are reconstructed using a cone algorithm [58,71,72] 
with cone sizes, AR(= \J Arf + A</> 2 ), of 0.3 and 0.5. 
The cone size of AR = 0.3 is used only in the level- 
2 trigger and for certain aspects of the all-jets analy- 
sis (see Sec. VIII); all other analyses use a cone size of 
AR = 0.5 to maximize the efficiency for reconstructing 
ti events. The algorithm is executed as follows. Starting 
from an £/r-ordered list of calorimeter towers, the tow- 
ers within AR w 0.3 and with Et > 1 GcV are grouped 
into preclusters. The energy within a given cone (0.3 
and 0.5 for the analyses presented here) centered on the 
precluster is summed, and a new "i?T-wcightcd" center 
is obtained. Starting with this new center, the process 
is repeated until the center stabilizes. A jet is required 
to have Et > 8 GeV. If two jets share energy, they are 
combined or split, depending on the fraction of Et shared 
relative to the Et of the lower Et jet. If the shared frac- 
tion exceeds 50%, the jets are combined. 

The jet energy resolution is obtained from studies of 
Et balance in dijet and photon+jet data in different rj 
regions [58]. As shown in Fig. 9, the fractional resolution 
(<t(Et) I Et) in the central region varies from 20% at a jet 
E T of 30 GeV to 8% at a jet E T of 100 GeV. Resolutions 
in the other detector regions are similar. The absolute 
jet energy scale is discussed in the following section. 

D. Missing E T (rfr) 

Neutrinos escape the detector without interacting. 
Similarly, muons pass through the calorimeter deposit- 
ing very little energy. The presence of a high-energy 
neutrino can be inferred from an imbalance in transverse 
energy/momentum as measured in the calorimeter and 
muon systems. 

Missing transverse energy in the calorimeter, ^f? 1 , is 
defined as 

Ef = ^E^ 2 + Ef\ (4.6) 



where 



18 



0.2 



W 0.15 H 



0.1 



0.05 - 




ti I < 0.5 



50 100 150 200 250 300 350 400 450 



Scale Corrected Jet E, 



(GeV) 



FIG. 9. Jet E T fractional resolution for j?7| < 0.5. The 
circles and solid line correspond to the nominal resolution; the 
dotted lines are the systematic uncertainty on the resolution 
measurement. The stars correspond to resolutions obtained 
from HERWIG+D0GEANT Monte Carlo and are used at high 
Et where dijet data are not available. 
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(4.7) 



(4.8) 



The first sum is over all cells in the calorimeter and ICD, 
and the second sum is over the corrections in Et ap- 
plied to all electrons and jets in the event (see Sec. A). 
The missing transverse energy ($r) resolution of the 
calorimeter is parameterized as [58] 



cr($r al ) = 1.08 GeV + 0.019 ^E T , 



(4.9) 



where Et is the scalar Et, which is defined to be the 
scalar sum of all calorimeter cell Et values. 

For events that contain muons, the transverse momen- 
tum of the muon is subtracted from IfZp to compute the 
total missing Et'- 



e x = e^-J2p£i 



/cal 



(4.10) 



(4.11) 



V. EVENT SIMULATION 

The ti signal efficiencies and several rare background 
rates are computed via Monte Carlo methods. The pri- 
mary event generator for the signal is herwig [54] with 
CTEQ3M [73] parton distribution functions (pdf). Tests 
were also performed with three values of Aqcd, and us- 
ing the MRSA' pdfs [74], but no significant variation in 
ti acceptance was seen, herwig chooses the momenta 



out of the initial hard scattering according to matrix ele- 
ment calculations and models initial and final state gluon 
emission using leading- log QCD evolution [75] . Each top 
quark is then made to decay into a W boson and a b 
quark, and the final state partons are hadronized into 
jets. Underlying spectator interactions are also included 
in the model. As a cross-check, acceptances were also 
computed using the isajet [76] event generator (also us- 
ing the CTEQ3M pdfs), and the difference between the 
two is incorporated into the systematic uncertainties on 
a per channel basis (see Sees. IX A 8 - IX A 9 for details). 
ISAJET also chooses the momenta out of the hard scat- 
tering based on matrix element calculations, but models 
the initial and final state gluon emission using Feynman- 
Field fragmentation [77]. 

herwig was chosen as the primary generator because 
it provides good agreement with data in D0's color co- 
herence [78] and jet-shape [79] analyses. As discussed in 
Sec. X, within available statistics, the leptonic top candi- 
dates found in the current analysis are in good agreement 
with expectations from herwig. However, it should be 
noted that version 5.7 of herwig (the version used for 
the present analyses) is based on leading-log matrix ele- 
ments, and is therefore not in complete agreement with 
higher-order predictions [80,81]. 

herwig and isajet samples were generated with top 
quark masses between 90 and 230 GeV/c 2 . To increase 
event-processing efficiency, two samples were made for 
each mass and generator: (1) one in which both of the 
W bosons were required to decay leptonically (e,/i, r), 
from which only those that resulted in a final state of 
ee, /Xyit, or e/j, were kept, and (2) one in which one of the 
W bosons was forced to decay leptonically (e, fi, r), from 
which those with no final state electrons or muons were 
rejected, as were one-half of the dilepton events (in order 
to preserve the proper branching ratios). 

For the dilepton channels, backgrounds from Z — > rr, 
Z — > jU/i, WW, WZ, and Drell-Yan production are de- 
termined with PYTHIA [82] and with isajet, and the dif- 
ference used as a measure of systematic uncertainty. 

Background levels from VF+jets production are deter- 
mined from data. However, as discussed in Sec. VII A, 
shape information from the vecbos [55] Monte Carlo 
program is used to determine the survival probability 
for the latter stages of the £+jet/topo analyses. For the 
H+jets/n analysis (see Sec. VII B), vecbos is used to 
determine the Z — > fifi background. In both cases, the 
CTEQ3M [73] pdfs are used, vecbos incorporates the 
exact tree-level matrix elements for W and Z boson pro- 
duction, with up to four additional partons, and supplies 
the final state partons. In order to include the effects of 
additional radiation and underlying events, and to model 
the hadronization of the final state particles, the VECBOS 
output is passed through herwig's QCD evolution and 
fragmentation stages. Since herwig requires information 
about the color labels of its input partons, both programs 
were modified to assign color and flavor to the gener- 
ated partons. The flavors are assigned probabilistically 
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by keeping track of the relative weights of each diagram 
contributing to the process. The color labels are assigned 
randomly. To estimate the systematic uncertainty, sam- 
ples were also generated using isajet, instead of herwig, 
to fragment the VECBOS partons. 

The output of an event generator is typically pro- 
cessed through a GEANT [70] simulation of the detector 
(D0GEANT). However, such a detailed simulation is ex- 
tremely computationally intensive and does not allow for 
generation of the necessary high-statistics samples. As a 
compromise, the full D0GEANT simulation was run on a 
large sample of electrons and hadrons, and the resultant 
calorimeter showers were stored in a library [67]. These 
showers are binned in five quantities representing the in- 
put particle: 

• z vertex position (6 bins) 

• T] (37 bins matching calorimeter segmentation) 

• momentum (7 bins) 

• <p region: The calorimeter is largely symmetric in 
4>, the exceptions being the cracks between mod- 
ules in the central electromagnetic calorimeter and 
the region through which the Main Ring passes in 
the hadronic calorimeter. Hence, there are only 
two bins in </>, representing the "good" and "bad" 
regions. 

• particle type: Energy depositions in the calorime- 
ter for electrons/photons and hadrons are stored 
separately in the shower library. 

A total of 1.2 million events was used to populate the 
library. As events are sent through the library version of 
the simulation, a shower from the library is selected to 
model the calorimeter response of each individual parti- 
cle. The total energy of the shower is scaled by the ratio 
of the energy of the particle to be simulated to that of the 
library particle which created the shower. Since the full 
D0GEANT simulation for muons is not as time-consuming 
(owing to their minimum- ionizing nature), muons are not 
included in the shower library but are instead tracked 
through the detector just as in non-library version of the 
simulation. 

For the muon system, the efficiency is overestimated 
and the resolution is underestimated by D0GEANT. The 
next step in the simulation procedure therefore smears 
the muon hit timing information so that the Monte Carlo 
hit position resolution matches that in Z — > fifi data, and 
randomly discards hits to model the chamber inefficiency 
more accurately. In addition, the muon-system geometry 
in the Monte Carlo is misaligned in order to reproduce 
the correct overall momentum resolution. 

For several of the analyses, a final step in the sim- 
ulation models the level 1 and level 2 triggers (trigger 
simulator). As discussed in Sec. Ill, the level 1 trigger is 
a collection of hardware elements interfaced to an AND- 
OR network. The level 1 simulation therefore consists 



of simulated trigger elements and a simulated AND-OR 
network. Level 2 is a software trigger that runs in the 
online data acquisition environment. The level 2 sim- 
ulation consists of exactly the same code but has been 
ported to the offline environment. The level 1 and level 2 
simulations are typically used as a single entity, referred 
to simply as the trigger simulator. 

VI. ANALYSIS OF DILEPTON EVENTS 

As discussed in Sec. I, the ee, e/i, and /i/i dilepton sig- 
natures are characterized by two isolated high-py charged 
leptons, I$t, and two or more jets (from the b quarks 
and initial and final state radiation). Figures 3 and 5 
show Monte Carlo distributions for the lepton and jet 
Et/pt and \r]\, and the $ T expected in ti — > e/i events 
with rat = 170 GeV/c 2 . Background events with a sim- 
ilar topology are relatively rare and arise primarily from 
Drell-Yan production of (Z/7)+jets, TVIV+jets, and lep- 
tonic IV+jets events in which the second lepton arises 
from the misidcntification of one of the jets. Therefore, 
requirements based on the above characteristics form the 
initial selection for all three channels (see Tables 10, 12, 
and 14). Additionally, for the ee and /x/t channels there 
are cuts designed to reject Z — > ee, /i/i events. 

To attack the remaining background, variables were 
selected based on a series of cut optimization studies. 
These are designed to maximize the significance, defined 
as S = signal/ ^background, and result in the introduc- 
tion of a new transverse energy variable, defined as 

H^ = J2 E t + (leading electron E T ) (6.1) 

jets 

for the ee and e/i channels and as 

H T = E t (6.2) 

jets 

for the mi channel. The sums are over all jets with 
E T > 15 GeV and |r?| < 2.5. The optimized H T and 
Ht cut values are given in the event selection tables in 
Sees. VIA, VI B, and VIC. An additional result of the 
optimization studies was the requirement that, for the 
ee, e/i, and /t/i channels, there should be at least two jets 
with Et > 20 GeV. As discussed below, both of these re- 
quirements are very effective in distinguishing ti events 
from background. 

In addition to the ee, e/i, and /i/i channels, a new 
channel was introduced that is designed to catch dilep- 
ton events in which one of the leptons cither fails the 
Pt requirement or escapes detection (perhaps by pass- 
ing through an uninstrumcntcd region of the detector). 
This "ez/" channel selects events that contain one high-py 
electron, significant missing transverse energy, and two 
or more jets. The analogous \iv channel has not been 
considered. 
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Acceptances for all four dilepton channels were com- 
puted from Monte Carlo events generated by the her- 
WIG program for 24 top quark mass values (raj = 90 - 
230 GcV/c 2 ) and then passed through the full D0 detec- 
tor simulation (see Sec. V). The expected number of ti 
events passing the selection for a given channel is 

N = a t - t {m t ) /2 A{ i> ^ mt ^ ' £ « ^ 6 - 3 ) 

i— runs j—dct 

where a t i is the theoretical ti cross section at a 
top quark mass of m t [45], dj is the integrated 
luminosity for run i and pair of lcpton detector 
regions j (for ee j=CC+CC,CC+EC,EC+EC, for 
e/x j=CC+CF,CC+EF,EC+CF,EC+EF, and for nfj, 
j=CF+CF,CF+EF,EF+EF), and the acceptance, A, is 

A = e trig • e pi d • £soi G-B, (6.4) 

where £trig(*, j, m t) is the trigger efficiency, e p a(i,j) is 
the efficiency for identifying the two leptons, s se i(i,j, m t ) 
is the efficiency of the selection criteria, G(i,j) is the 
geometric acceptance, and B is the branching fraction 
for the sample being studied. Trigger efficiencies are ob- 
tained from data or Monte Carlo, depending on the chan- 
nel, and are discussed in greater detail below. Particle 
identification efficiencies are obtained from data in the 
case of electrons (as discussed in Sec. IV A), and from 
a combination of data and Monte Carlo in the case of 
muons (as discussed in Sec. IV B). The selection efficien- 
cies e sc \ and geometrical acceptances G are calculated 
from Monte Carlo. As will be discussed in Sec. X, it 
is the acceptance, rather than the expected number of 
ti events, that is used to calculate the ti cross section. 
Typical values for acceptance, often denoted as the "ef- 
ficiency times branching fraction" (e x B), for all eight 
leptonic channels, are tabulated in Sec. X for seven values 
of top quark mass. The numbers of ti events expected 
in the four dilepton channels are tabulated in Sees. VIA, 
VI B, VIC, and VI D, for the same set of top quark 
masses. Systematic uncertainties on the acceptances are 
discussed in Sec. IX. 

Whenever possible, backgrounds are measured directly 
from data. If not, then the backgrounds are determined 
from Monte Carlo events in which the initial cross sec- 
tions are normalized either to measured or theoretical 
values: 

B = cr b kg (etrig • £ P id • £sci ■ G^j ■ Ci^ (6.5) 

i— runs j— det ^ 

where Cbkg is the measured or theoretical cross section 
for the background under consideration. 

A. The ee channel 

The signature for an event in the ee channel consists of 
two isolated high-i^T electrons, two or more jets (from the 



b quarks and initial and final state radiation) , and signif- 
icant iftrp (from the neutrinos) . The trigger for this chan- 
nel was (depending on run period) ELE-JET(la), ele- 
jET-HiGH(lb), or ele-jet-higha(Ic), requiring an elec- 
tron, 2 jets, and $ T at level 2 (see Sec. Ill for details). 
As discussed in Appendix C, for this analysis Main- Ring 
events were corrected and not rejected. Over the com- 
plete Run 1 data set, these triggers provided a total inte- 
grated luminosity of 130.2 ± 5.6 pb _1 . The event sample 
passing these triggers consists primarily of misidentified 
multijet and heavy flavor events. 

The backgrounds to this signature arise from Drell-Yan 
(Z/^f*) production that results in a dielectron final state 
(Z — ► ee, Z — > tt — ► ee, and 7* — > ee), WW — > ee, 
and multijet events containing one or more misidenti- 
fied electrons. The latter background consists primarily 
of W(-^> e^)+3 jet events in which one of the jets is 
misidentified as an electron. 

The offline selection cuts and their cumulative effect 
are summarized in Table 10. After passing the trig- 
ger requirement, events are required to have 2 electrons 
(loose electron identification, see Sec. IV A) with Et > 20 
GeV and \t]\ < 2.5. This initial selection has an accep- 
tance (e x B) of (0.26 ± 0.03)% (for m t = 170 GcV/c 2 ), 
and essentially eliminates any background from heavy 
flavor production and reduces the QCD multijet back- 
ground to a small fraction of the remaining dominant 
background from Z — > ee. The number of Z + n jet 
events is proportional to a™, and a similar steep falloff 
in jet multiplicity is observed for the other backgrounds 
present at this stage. Requiring 2 jets with Et > 15 GeV 
and 1 77 1 < 2.5 significantly reduces backgrounds from Z 
boson, Drell-Yan and WW production, and QCD mul- 
tijet events. Most of these (Z, Drell-Yan, and QCD 
multijet) do not contain high-p^ neutrinos. Therefore, 
a hard cut on the rfx brings these events to an even 
more manageable level. At this point the background is 
still dominated by Z — > ee events, so the next step re- 
quires that the dielectron invariant mass not be within 
the mass window of the Z boson (see Table 10). How- 
ever, since Z — » ee events have no real $ T , this cut is 
only made for events with $ T < 40 GeV, thereby re- 
claiming a considerable amount of ti efficiency. The final 
two cuts, H T > 120 GeV and N icts > 2 with E^ 1 > 20 
GeV and \rj' iet \ < 2.5, are obtained through the optimiza- 
tion procedure discussed in Sec. VI, and provide rejection 
against the remaining background from Z — > tt, WW, 
and Drell-Yan production, and QCD multijet events. Ta- 
ble 10 shows the number of data events, expected signal 
(m t = 172.1 GeV/c 2 ), and expected background surviv- 
ing at each stage of the selection. It is clear from this 
table that the Ifr requirement greatly reduces the back- 
ground. This is shown in Fig. 10, where Ip T is plotted 
vs M ee for all the major backgrounds (a-d), for ti Monte 
Carlo (e), and for data (f). Because of the presence of 
two neutrinos, the WW background is not reduced much 
by the selection on $ T . It is, however, reduced signif- 
icantly by the jet and H T requirements. The effect of 
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TABLE 10. Number of observed and expected ee events passing at each cut level of the offline analysis. Expected number of 
ti events are for m t = 172.1GeV/c 2 . Uncertainties correspond to statistical and systematic contributions added in quadrature. 



Number of ee events at each cut level 







Total 


mis- id 


physics 






Data 


sig + bkg 


bkg 


bkg 


ti 


2e, E C T > 20 GeV, + e id + trig 








4168 ± 1243 


1.9 ±0.3 


+ 2 jets, > 15 GeV 


112 


125 ± 36 


9.0 ±0.08 


114 ±35 


1.8 ±0.3 


+ ^ > 25 GeV 


3 


3.2 ±1.9 


0.23 ±0.06 


1.5 ±1.9 


1.5 ±0.3 


+ $r al > 40 GeV or 












M ee < 79 GeV/c 2 or M ee > 103 GeV/c 2 


2 


2.3 ±0.5 


0.22 ±0.06 


0.62 ±0.21 


1.5 ±0.3 


+ 2 jets, > 20 GeV 


2 


1.9 ±0.4 


0.20 ±0.05 


0.39 ±0.12 


1.4 ±0.3 


+ Wt > 120 GeV 
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1.7 ±0.2 


0.20 ±0.05 


0.28 ±0.09 
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FIG. 10. Scatter plots of ~$ T vs M ee for the ee channel: 
(a) Z — > ee events, (b) Z — > tt — > ee MC events, (c) QCD 
multijet events, (d) WW -» ee MC events, (e) ti -» ee MC 
signal (m t = 172.1 GeV/c 2 ), and (f) data. The signal region 
is defined as being above the solid line in each plot. 



the cut on WW events can be seen in Fig. 11(b), 
which gives the iff. distribution for ti — ► efi events, but 
is very similar to that for ti — ► ee events. After the above 
selection, only one ee candidate remains. 

The Z — > ee background is determined entirely from 
data. As noted above, Z(— > ee)±jets events have no 
real IpT, and due to the excellent electron momentum 
resolution, any I$t observed in the detector will arise 
from mismeasurement of jet Et and other noise in the 
calorimeter. Due to the extremely high rejection power 
of the IjjT requirement on Z — > ee±jet events, a rfr mis- 



TABLE 11. Expected number of ee signal and background 
events after all cuts in 130.2 pb _1 . Uncertainties are statis- 
tical and systematic contributions added in quadrature. The 
systematic uncertainty on the total background includes cor- 
relations among the different background sources. 



Expected number of ee events in 130.2 pb 



top MC m t (GeV/c 2 ) 




140 


2.34 ±0.34 


150 


1.96 ±0.29 


160 


1.62 ±0.23 


170 


1.25 ±0.18 


180 


1.02 ±0.15 


190 


0.79 ±0.12 


200 


0.62 ± 0.09 


Z — > tt — > ee 


0.08 ± 0.06 


multijet (mis-id e) 


0.20 ± 0.05 


Z -> ee 


0.06 ±0.01 


WW -> ee 


0.09 ± 0.03 


DY -> ee 


0.06 ± 0.03 


Total background 


0.48 ±0.10 
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measurement rate is determined from a sample of QCD 
multijet data selected to closely match the jet require- 
ments in this analysis: > 2 jets, E T > 20 GeV, H T > 70 
GeV (where the remaining 50 GeV contribution to the 
> 120 GeV is assumed to originate from the highcst- 
Et electron). The fraction of events in this sample that 
passes the Ifr > 25 GeV requirement is taken as the I$t 
mismeasurement rate (i.e., the fraction of the time that 
the detector resolution will result in a false I(it signal). 
Due to a slight dependence on jet multiplicity, the Ifr 
mismeasurement rate is determined as a function of the 
TpT cut and number of jets n in the event and is found 
to be (1.02 ±0.09)% for n = 2, (0.86 ±0.02)% for n = 3, 
and (1.12 ± 0.02)% for n = 4 for # T > 25 GeV; and 
(0.20 ± 0.04)% for n = 2, (0.14 ± 0.01)% for n = 3, and 
(0.17 ± 0.01)% for n = 4 for # T > 40 GeV. These fac- 
tors are then applied to the number of dielectron events 
that pass all selection requirements (including the Z bo- 
son mass window cut), except for that on tf,T, to obtain 
the total expected Z — > ee background of 0.058 ± 0.013 
events. The systematic uncertainty on this determina- 
tion is discussed in Sec. IX. 

The background from multijet events is also obtained 
entirely from data. The probability for an extra-loose 
electron to pass the loose electron identification cri- 
teria (see electron misidentification rate discussion in 
Sec. IV A) is applied to both the full Run 1 sample (not 
including Main-Ring, MR, events) of dielectron events in 
which one electron candidate passes the loose identifica- 
tion and the other fails the loose identification but passes 
the extra-loose identification, and to that where both 
electron candidates fail the loose identification but pass 
the extra-loose identification. The resultant misidentifi- 
cation background is then scaled up by the (nonMR + 
MR) /nonMR luminosity ratio to account for the misiden- 
tification background expected in the MR data. 

Backgrounds from Z — > tt — > ee, WW — > ee, and 
7* — > ee are obtained from PYTHIA and isajet Monte 
Carlo samples via Eq. 6.5, and are normalized either to 
experimental or theoretical values. 

The Z — ► tt — > ec Monte Carlo samples are normal- 
ized to DCs Z boson cross section measurement and its 
measurement of (to obtain more Z+jets events and 
thus enhance the final statistics, generator-level cuts are 
placed on pf) [83,84] and corrected for the Z — > tt and 
t — > ev e v T branching fractions [85]. The 7* — > ee Monte 
Carlo sample is likewise normalized to DCs measure- 
ment of the Drell-Yan (7* — > ee) cross section in the 
dielectron mass range 30 GeV/c 2 < M ee < 60 GeV/c 2 
[86]. The WW — > ee Monte Carlo samples are normal- 
ized to theory [87], and a 10% uncertainty is assigned [88]. 

For the Z — > tt — > ee background, the associated jet 
spectrum in PYTHIA, herwig, and isajet does not agree 
with that found in the Z — ► ee data. This is corrected 
by incorporating the jet cut survival probabilities from 
the Z{— ► ee)±jet data (where the Ht cut is taken as 70 
GeV, as in the mismeasured I^t calculation) rather than 
from Monte Carlo. 



As described in the previous section, the tt acceptances 
are computed via Eq. 6.4 using Monte Carlo events gen- 
erated with herwig and passed through the D0 detector 
simulation (see Sec. V). The trigger efficiency is obtained 
from Z — > ee data but cross checked with the trigger sim- 
ulator (see Sec. V). Both approaches result in a trigger 
efficiency of 99 ± 1% [68]. 

The acceptance values after all cuts for seven top quark 
masses (for all channels) are given in Sec. X. The ex- 
pected numbers of tt events, determined via Eq. 6.3, are 
given in Table 11 for each of these seven masses. Fi- 
nally, a cross section of 2.4 ± 4.6 pb is obtained for the 
ee channel. 

To test the robustness of the background predictions, 
comparison is made of data and expectations in regions 
dominated by background (i.e., at earlier steps along the 
selection chain). Making use of Eqs. 6.3 - 6.5 for the 
different stages of the selection, Table 10 shows that the 
expectation from background and tt compare well with 
what is observed in the data at the various stages of the 
selection procedure. 

B. The e/i channel 

The signature for an event in the e/i channel consists 
of one high-i?T isolated electron, one high-py isolated 
muon, two or more jets (from the b quarks and initial 
and final state radiation), and significant $ T (from the 
neutrinos). The trigger for this channel required one of 
the following level 2 terms to be satisfied: 

• ELE- JET ( la), ELE- JET-HIGH (lb), Or ELE- JET-HIGHA 

(lc), which required an electron, 2 jets, and $ T . 

• MU-ELE(la and b) or MU-ELE-HIGH(lc), which re- 
quired an electron and a muon. 

• MU-JET-HiGH(la and b) or mu-JET-CENt(Ic), 
which required a muon and a jet. 

Details of these triggers are discussed in Sec. III. Main- 
Ring events are not included in this analysis. Over the 
complete Run 1 data set, these triggers provided a total 
integrated luminosity of 112.6 ± 4.8 pb -1 . 

The backgrounds to this signature arise from Drell- 
Yan production of tt which can lead to e/i final states 
(Z — ► tt — > ep and 7* — ► tt — > e/i), WW — > e/i, 
and multijet events containing an isolated muon and a 
misidentified electron. The latter background consists 
primarily of W(— > ^v)+3 jet events, where one of the 
jets is misidentified as an electron. Backgrounds con- 
taining a real electron and a misidentified isolated muon, 
and those containing both a misidentified electron and 
a misidentified isolated muon were discussed in Ref. [58] 
and found to be negligible. 

The offline selection cuts and their cumulative effect 
are summarized in Table 12. After passing the trigger 
requirement, events are required to have > 1 electron 
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TABLE 12. Number of observed and expected e/t events passing at each cut level of the conventional analysis. Expected 
number of tt events are for writ = 172.1GeV/c 2 . Uncertainties correspond to statistical and systematic contributions added in 
quadrature. 



Number of e/i events passing cuts 







Total 


Mis-id 


Physics 






Data 


sig + bkg 


bkg 


bkg 


tt 


Et > 15 GeV, p£ > 15 GeV 












+ e id ± fi id + trig 


130 


93 ±7 


50 ±2 


39 ±6 


4.3 ±0.9 


+ Ai?(/i,jet) > 0.5 


60 


59 ±6 


17.8 ±0.9 


38 ±6 


3.4 ±0.7 


+ $ T > 10 GeV 


41 


38 ±3 


13.5 ±0.7 


21.4 ±3.3 


3.4 ±0.7 


+ ^ > 20 GeV 


22 


21.8 ±2.2 


4.5 ±0.4 


14.0 ±2.1 


3.2 ±0.6 


+ AR(e,fi) > 0.25 


20 


19.5 ±2.2 


2.3 ±0.3 


14.0 ±2.0 


3.2 ±0.6 


+ 2 jets, > 15 GeV 


4 


3.4 ±0.6 


0.32 ±0.14 


0.34 ±0.09 


2.7 ±0.6 


+ tff. > 100 GeV 


4 


2.8 ±0.5 


0.11 ±0.12 


0.24 ±0.08 


2.5 ±0.5 


+ H e T > 120 GeV 


3 


2.6 ±0.5 


0.08 ±0.12 


0.20 ±0.08 


2.3 ±0.5 


+ 2 jets, E^ > 20 GeV 


3 


2.5 ±0.5 


0.08 ±0.12 


0.19 ±0.10 


2.2 ±0.5 



(loose electron identification, see Sec. IV A) with Et > 15 
GeV, \r)\ < 2.5 and > 1 muon (loose muon identifica- 
tion, see Sec. IV B) with pt > 15 GeV/c. This initial 
selection has an acceptance (e x B) of 0.68 ± 0.15% for 
m t = 170 GcV/c 2 . At this stage, the background is dom- 
inated by QCD multijet events containing a jet misiden- 
tified as an electron and a non-isolated muon from the 
scmi-leptonic decay of a b or c quark. This background 
is reduced significantly by requiring the muon to be iso- 
lated, Ai?(/i,jet) > 0.5. To further reduce the misidenti- 
fication background, the next two steps require $ T > 10 
GeV and iff 1 > 20 GeV. The cut on iff 1 is particu- 
larly effective against background from W(-^ /i^)±jets 
events (where one of the jets is misidentified as an elec- 
tron) due to the fact that iff provides a measure of 
the transverse momentum of the W boson since both 
of its decay products deposit little or no energy in the 
calorimeter. Studies also show that QCD multijet events 
that contain a highly electromagnetic jet (misidentified 
as an electron) which gives rise to an isolated muon from 
the semi-leptonic decay of a b or c quark, can easily enter 
this analysis (as can W(— > /xi/)+jets events where there 
is significant bremsstrahlung from the muon as it passes 
through the EM calorimeter) . Such events typically have 
the e and /i very close in (77, <j>) space, and a require- 
ment of AR(e, fi) > 0.25 effectively eliminates this class 
of misidentification background. 

After the above requirements, the background is pri- 
marily from Z — > tt — ► e/i events and, to a lesser extent, 
from WW — ► e/i events. The jets associated with these 
processes arise from initial state radiation (recoil) and 
are therefore softer in Et than the b jets in a tt event. In 
addition, as noted above (see Sec. VIA), the number of 
Z + n jet events is proportional to a™, and a similar steep 
falloff in jet multiplicity is observed for the Drell-Yan 
(and presumably WW) backgrounds. Requiring two jets 
with > 15 GeV and \i? et \ < 2.5 significantly reduces 



these backgrounds and that from QCD multijet produc- 
tion. The final cuts on H T > 120 GeV and Ar jets > 2 for 
E^ > 20 GeV and \^ ct \ < 2.0 are obtained through the 
optimization procedure discussed in Sec. VI and provide 
further rejection against the remaining backgrounds. Af- 
ter the above selection, three e/i candidates remain in the 
data. 

Table 12 shows the number of data events, expected 
signal (m t = 172.1 GeV/c 2 ), and expected background 
surviving at each stage of the selection. It is clear from 
this table that the H T cut is the most effective cut during 
the final stages of the analysis. This is also shown in 
Fig. 11, where the H T distributions are given for the 
three major backgrounds (a-c), for tt Monte Carlo (d), 
and for data superimposed on the total background and 
expected tt signal (e). 

As in the case of the ee channel, the background from 
multijet events is obtained entirely from data. The prob- 
ability for an extra-loose electron to pass the loose elec- 
tron identification criteria (see misidentification rate dis- 
cussion in Sec. IV A) is applied to the full Run 1 sample of 
e/i events, where the electron candidate passes the extra- 
loose electron identification but fails the loose electron 
identification, with all the other kinematic cuts applied. 
As shown in Table 13, the QCD multijet (misidentified 
e) background is determined to be 0.08 ± 0.12 events. 

Background estimates for Z — > tt — > e/i, WW — > e/i, 
and 7* — ► e/i events are obtained via Eq. 6.5 using nor- 
malized pythia and ISAJET Monte Carlo samples. The 
Z — > tt — > e/i Monte Carlo samples are normalized to 
DCs measurement of a(pp —fZ + X)B(Z — > ee) and 
the associated measurement of [83,84], and incorpo- 
rate the Z — > tt, t — > ev e v Tl and r — > [iVy,y T branching 
fractions [85]. The 7* — ► tt Monte Carlo sample is like- 
wise normalized to D0's measurement of the Drell-Yan 
(7* — ► ee) cross section in the dielectron mass range 30 
GeV/c 2 < M ee < 60 GeV/c 2 [86] also incorporating the 
t — ► ev e v r and t — > iiv^v-r branching fractions [85]. The 
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TABLE 13. Expected number of efx signal and background 
events in 112.6 pb -1 after all cuts in the conventional analy- 
sis. Uncertainties are statistical and systematic contributions 
added in quadrature. The systematic uncertainty on the total 
background includes correlations among the different back- 
ground sources. 
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FIG. 11. 

pected background (hatched), expected signal (open), and 
data (solid) after all cuts except > 120 GeV (shown by 
solid vertical line) and 2 jets with E 1 ^ * > 20 GeV (correspond- 
ing to line 6 of Table 12). Plots (a) - (c) show the individual 
contributions of the three leading backgrounds and give the 
expected number of events. Plot (d) gives the expected tt 
contribution (m t = 170 GeV/c 2 ), and plot (e) overlays the 
total expected background, expected signal, and data (xO.l). 



tt MC m t (GeV/c 2 ) 




140 


4.07 ±0.88 


150 


3.32 ±0.72 


160 


2.77 ±0.60 


170 


2.29 ±0.49 


180 


1.84 ±0.40 


190 


1.48 ±0.32 


200 


1.12 ±0.24 


Z — > tt — + efi 


0.10 ±0.09 


QCD multijet (mis-id e) 


0.08 ±0.12 


WW-> en 


0.08 ±0.02 


DY — > TT — ► 6/i 


0.006 ± 0.004 


Total background 


0.26 ±0.16 



WW — > efi Monte Carlo samples are normalized to the- 
ory [87], and a 10% uncertainty assigned [88]. 

As for the ee channel, the Z — > tt — ► eji Monte Carlo 
samples are not used to model the jet and ti^ require- 
ments. Instead, survival probabilities for these cuts are 
obtained from Z(—> ee)±jet data. 

The tt acceptances are computed via Eq. 6.4 using 
Monte Carlo events that are generated with herwig and 
passed through the D0 detector simulation (see Sec. V). 
The trigger efficiency is obtained from the trigger sim- 
ulator and is dependent on the detector region of the 
electron and muon, giving (95 ± 5)% for CC(e)CF(/i), 
(93±5)%forEC(e)CF(/i), (90±4)% for CC(e)EF(/i), and 
(93 ± 5)% for EC(e)EF(/x). The acceptance values after 
all cuts for seven top quark masses (and for all channels) 
are given in Sec. X. The expected number of tt events 
passing this selection is determined via Eq. 6.3 and are 
given in Table 13 for these same seven masses. Finally, a 
cross section of 6.8±4.6 pb is obtained for the ey, channel. 



C. The channel 

The signature for an event in the /x/x channel consists 
of two isolated high-py muons, two or more jets (from the 
b quarks and initial and final state radiation) , and signif- 
icant $ T (from the neutrinos). The trigger for this chan- 
nel required one of the following level 2 terms to be sat- 
isfied: MU-JET-HIGH(la and lb), MU-JET-CAL(lb), MU- 
JET-CENT(lb and lc), or MU-JET-CENCAL(lb and lc). 
Each of these required a muon and one jet at level 2 
(see Sec. Ill for details). Main- Ring events arc not in- 
cluded in this analysis. Over the complete Run 1 data 
set, these triggers provided a total integrated luminosity 
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of 108.5 ±4.7 ptT 1 . 

The backgrounds to this signature arise from Drcll- 
Yan production with dimuon final states (Z — > /lit, 
Z — > tt — > /x/it, and 7* — » /j/j), ITW — ► iz/z, and multijct 
events containing misidcntified isolated muons. The lat- 
ter background consists primarily of four-jet events where 
the semi-leptonic decay of b and/or c quarks result in 
two muons that pass the isolation requirement, and of 
W{-^ fiv)+3 jet events where one of the jets gives rise 
(through the semi-leptonic decay of a b or c quark) to a 
muon that passes the isolation requirement. 

The offline selection cuts and their cumulative effect 
are summarized in Table 14. After passing the trig- 
ger requirement, events are required to have two muons 
(loose muon identification, see Sec. IV B) with pr > 15 
GeV/c and |^| < 1.0 (\t]\ < 1.7 in Run lbc postclean) 
and one jet with > 20 GeV and \rf\ < 2.5. This initial 
selection has an acceptance (e x B) of 0.35% (m t = 170 
GeV/c 2 ). At this stage, the dominant background is from 
cosmic rays. This is minimized by rejecting tracks that 
are back-to-back in both -q and (f>: 

AtOiufo) < 165° for \v(pi) + V(P2)\ < 0.3. (6.6) 

It is necessary to exclude background from J/ip — > fj.fi. 
As discussed below, the muon momentum resolution pro- 
hibits an efficient cut on M MM at the Z boson mass peak. 
However, at lower muon pt, it is an effective quantity 
and is used to reject low-mass pairs resulting from high- 
Pt J/ip production with recoil jets: M MM > 10 GeV/c 2 
is required. At this stage, the background is dominated 
by QCD multijet events rich in heavy flavor with muons 
originating from semi-leptonic decays of b or c quarks. By 
requiring both muons to be isolated (Ai?(/i,jet) > 0.5), 
this background is reduced to a negligible level. The 
remaining background is mainly from events containing 
isolated dimuons from Z/-f* and WW production. The 
jets associated with these processes arise from recoil and 
are thus softer in Et than the b jets in a tt event. Also, 
as noted in Sec. VIA, the number of Z+ > n jet events 
is proportional to a™, and a similar steep falloff in jet 
multiplicity is observed for the Drell-Yan and WW back- 
grounds. The next step in the analysis therefore requires 
a second jet with Et > 20 GeV and \rj\ < 2.5, reducing 
the dimuon background from these sources. The require- 
ment of Ht > 100 GeV is obtained through the optimiza- 
tion procedure, as discussed in Sec. VI, and provides fur- 
ther rejection against the remaining background, leaving 
only the contribution from Z — > /z/z at a non-negligible 
level. 

As noted above, because of limitations on the momen- 
tum resolution of the D0 muon system, the invariant 
mass peak of the Z boson is smeared and a simple cut on 
M M/J is ineffective in reducing this background. Instead, 
rejection is achieved using the result of a % 2 minimization 
procedure that involves a refitting of the muon momenta 
with a constraint that the transverse momentum of the 
dimuon system balance the remaining transverse energy 



in the event: 

(_>__ 1 ) 2 1 ) 2 

2 _ W yr) 

V Pfil / V Pfj,2 / 

+ ^W) (6J) 

with the constraint that = Mz- 

Ml = 2(p>» 2 -^ 1 .^ 2 ) (6.8) 

where p^i is the measured momentum for the z-th muon, 
p ^ is the fitted value ofp^i, c(— ) is the measured muon 

momentum resolution (see Eq. 2.1), E^ 1 and E^ 1 are the 
x and y components of E T al , and c(^ al ) and cr(^ al ) 
are their measured resolutions (see Eq. 4.9). This \ 2 
is minimized as a function of p ^ and p°^ 2 - An event 
is considered to be a Z — > pp candidate, and is thus 
rejected, if Prob(x 2 ) > 0.01. This procedure is also used 
to remove Z — > fip background from the tt — > p+jets + p 
tag channel (see Sec. VII B). 

Table 14 shows the number of observed events, ex- 
pected signal (for m t = 172.1 GeV/c 2 ), and expected 
background surviving at each stage of the selection. It is 
clear from this table that the Ht and Prob(x 2 ) cuts pro- 
vide significant background rejection in the final stages 
of the analysis. This is shown in Fig. 12, where Ht vs 
Prob(x 2 ) is plotted for Z — > /z/z and Z — > tt — > fifi MC 
events (a,b), for tt MC events (c), and for data (d). 

One tt — ► fifi candidate survives the above selection. 
Both muons in the event are central, and each track has 
the maximum of ten hits in the muon chambers, the case 
where the momentum resolution is best modeled and un- 
derstood. An interesting feature of this event is that all 
the muons and jets are in one hemisphere in cj> in the 
detector, leaving only $ T m the other half; this topology 
is highly unlikely to come from the main background of 
Z — > fifi production. 

The background from multijet events is determined en- 
tirely from data. The probability for a jet to give rise to 
an isolated muon is determined separately for the CF 
and EF regions of the muon system using a sample of 
of multijet events. These probabilities are then applied 
to the jets in a sample of muon (loose identification, sec 
Sec. IV B) + jet events to obtain the background ex- 
pected from W(-^> fiv)+jets, QCD multijet production, 
and Z — > tt — > fi + hadrons where the second muon 
originates from the semi-leptonic decay of a b or c quark 
from initial or final state radiation. 

In a manner analogous to the background calcula- 
tions used for the ee and efi channels, backgrounds from 
Z — > fiji, Z — ► tt — ► fifi, WW — > pfi 1 and 7* — > fifi 
are obtained via Eq. 6.5 from pythia and ISA jet Monte 
Carlo samples which are normalized to experimental or 
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TABLE 14. Number of observed and expected fifi events passing at each cut level of the offline analysis. Shown are results 
for Run lb+lc (CF-CF) only. Expected number of ti events are for nit = 172.1GeV/c 2 . Uncertainties correspond to statistical 
and systematic contributions added in quadrature. 



number of /j/j events passing cuts 







Total 


Mis-id 


Physics 






Data 


sig + bkg 


bkg 


bkg 


tt 


2tt, ptjn > 15 GeV/c, + it id 












+ trig + 1 jet, > 20 GeV 


606 






174 ± 50 


1.6 ±0.2 


+ A0(/2i,At2) < 165° for |J7 M1 ±r, M2 | < 0.3 


207 






146 ± 42 


1.5 ±0.2 


+ M m > 10 GeV/c 2 (J/V rej) 


165 


187 ±43 


40 ±9 


146 ± 42 


1.5 ±0.2 


+ Ai?(/i,jet) > 0.5 


105 


136 ± 39 


0.70 ± 0.33 


134 ± 39 


0.9 ±0.1 


+ 2nd jet, > 20 GeV 


19 


13.6 ±8.0 


0.22 ±0.10 


12.7 ±8.0 


0.72 ±0.09 


± H T > 100 GeV 


6 


5.1 ±3.3 


0.03 ±0.02 


4.5 ±3.3 


0.53 ±0.07 


+Z fit prob( X 2 ) < 1 % 
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0.9 ±0.3 


0.03 ±0.02 


0.42 ±0.16 


0.48 ±0.06 
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FIG. 12. Scatter plots of Ht vs Prob(x 2 ) for the fifi chan- 
nel: (a) Z — > fifi background, (b) Z —> tt — > ii/i background, 
(c) ti — > signal, and (d) data. The signal region is shown in 
the upper left corner of each plot (Prob(x 2 ) < 1%, Ht > 100 
GeV). 



theoretical values. In particular, the Z — > tt — ► fj.fi MC 
samples are normalized to the D0 Z boson cross section 
measurement but incorporate Z — > tt and r — ► iiv ll v T 
branching fractions from elsewhere [85]. Similarly, the 
7* — > fifi Monte Carlo sample is normalized to D0's mea- 
surement of the Drell-Yan (7* — > ee) cross section in the 
dielectron mass range 30 GeV/c 2 < M ee < 60 GeV/c 2 
[86] . The WW — ► Monte Carlo sample is normalized 
to theory [87] and a 10% uncertainty assigned [88]. 

As for the ee and e/i channels, the Z — > tt — > /j/j 
Monte Carlo samples are not used to model the jet and 
Ht requirmcnts. Instead, survival probabilities for these 



cuts are obtained from Z{— > ee)±jet data. 

As described in Sec. VI, the ti acceptances are com- 
puted via Eq. 6.4 using Monte Carlo events that are gen- 
erated with herwig and passed through the D0 detector 
simulation (Sec. V). The trigger efficiency is computed 
using data-derived trigger turn-on curves applied to ti 
Monte Carlo and is determined to be (95 ±5)%. The ac- 
ceptance values after all cuts for seven top quark masses 
(and for all channels) are given in Sec. X. The expected 
numbers of ti events passing this selection are determined 
via Eq. 6.3 and are given in Table 15 for these same seven 
masses. Finally, a cross section of 2.1 ±8.8 pb is obtained 
for the /i/i channel. 

To test the robustness of the background predictions, 
comparisons are made between the data and expecta- 
tions in regions dominated by background (i.e., at earlier 
steps along the selection chain). Equations. 6.3 - 6.5 
give, for the different stages of the selection, the results 
in Table 14, which show that that the expectation from 
background and ti compare well with what is observed in 
the data at the various stages of the selection procedure. 



D. The ev channel 

The ev channel is based on the assumption that one 
of the W bosons decays to ev and that the remaining 
ti decay products conspire to give rise to significant Jjf^ 
(> 50 GeV). As can be inferred from Figs. 3 and 4, this is 
most probable for ee and e/i events but will also occur in 
some fraction of the e±jets events. To eliminate overlap 
with the dilepton channels, it is further assumed that for 
e/it(ee) events, the muon (second electron) is either too 
low in pt(Et) to pass the selection or escapes detection. 
The signature for an event in the ev channel is therefore 
one, and only one, high-i/p electron, two or more jets 
(from the b quarks and initial and final-state radiation), 
and very large lfi T (from the neutrinos and possibly a 
lost lepton). The virtue of this channel is that it can 
recover some of the ti cross section not seen by the other 
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TABLE 15. Expected number of fifi signal and background 
events after all cuts in 108.5 pb _1 . Uncertainties are statis- 
tical and systematic contributions added in quadrature. The 
systematic uncertainty on the total background includes cor- 
relations among the different background sources. 



Expected number of /i/i events in 108.5 pb 



ti MC rm (GeV/c 2 ) 




140 


1.02 ±0.15 


150 


0.88 ±0.13 


160 


0.78 ±0.11 


170 


0.67 ±0.09 


180 


0.54 ±0.08 


190 


0.44 ± 0.06 


200 


0.33 ± 0.05 


Z — *■ TT — ► flfj, 


0.03 ± 0.03 


QCD multijet (mis-id fj>) 


0.07 ±0.01 


Z — > 


0.58 ±0.22 


WW -> lili 


0.007 ±0.004 


DY -> nix 


0.07 ± 0.04 


Total background 


0.75 ± 0.24 



channels. Indeed, investigating herwig ti Monte Carlo 
events (at m t = 170 GeV/c 2 ), the final ev sample is 
found to consist of one-half dilcpton (ee and e/i) events, 
one-third e±jets events, and one-sixth e + hadronic-tau 
events. 

The trigger for the ev channel required one of the fol- 
lowing level 2 terms to be satisfied (see Sec. Ill): 

• ELE-JET(la), ELE- JET-HIGH (lb) , EMl-EISTRKCC- 

MS(lb) or ele-jet-higha(Ic), all of which required 
an electron, 2 jets, and $ T - 

• MISSING-ET(lab) or missing-et-high(Ic), both of 
which required only very large ■ 

Note that Main-Ring events were not included in this 
analysis. Over the complete Run 1 data set, these trig- 
gers provided a total integrated luminosity of 112.3 ± 
4.8 pb" 1 . 

The primary backgrounds to this signature arise from 
W(-^ ev)+2 jet events and QCD production of three-jet 
events where one jet is misidentified as an electron and 
the $x is an artifact of jet Et mismeasurement. An ad- 
ditional source of background is WW + n jets production 
where one of the W bosons decays to ev and, in the case 
of n = or 1, the other W decays hadronically. Similarly, 
backgrounds from WZ ± n jets also contribute, but to a 
lesser extent. 

The offline selection cuts and their cumulative effects 
are summarized in Table 16. After passing the trigger re- 
quirement, events are required to have one electron (min- 
imal electron identification, see Sec. IV A) with Et > 20 
GeV and |r/| < 1.2. This channel differs from the other 
ti channels both in choosing its initial electron identifi- 
cation to be minimal (loose electron identification is re- 
quired at a later stage) and in the restriction of electrons 



to the CC region of the calorimeter (to suppress QCD 
multijet background, which increases in the forward re- 
gion). This initial selection has an acceptance (e x B) 
of (11.1 ± 3.2)% (for m t = 170 GcV/c 2 ). The next step 
requires $ T > 50 GeV to select high-$ T ti events, reject 
QCD multijet background, and decrease the number of 
W(-^ ev) and WW events. To further decrease these 
backgrounds, two jets with Et > 30 GeV and \-q\ < 2.0 
are required. At this stage the background is dominated 
by W{— * ev)+2 jet events and a cut on the e, $ T trans- 
verse mass, Mjf > 115 GeV, brings it down to approxi- 
mately one event. The transverse mass is defined by 

(e, $r) = \j (|4| ± \$ T \) 2 - (4 + t T f. (6.9) 

This cut is also effective against QCD multijet back- 
ground, being similar to the Ej, (= E T + $ T ) cut which 
will be described in Sec. VII A, and tends to reject events 
where the electron is parallel to the $ T in <j>. The back- 
ground that remains is dominated by 3-jet events, where 
one of the jets is misidentified as an electron and the $ T 
is an artifact of jet Et mismeasurement. A topological 
cut, A0(^t, 2 nd Et object) > 0.5 rad, rejects two-jet- 
like events where the $ T is aligned with one of the jets 
due to an upward fluctuation of the highest Et jet or a 
downward fluctuation of the second-highest Et jet. Note 
that the electron is treated as a jet in this Et ordering. 

The next step requires that the loose electron identifi- 
cation criteria be applied to all electron candidates and 
brings the remaining QCD multijet background down to 
an acceptable level. The final step in the selection re- 
quires, for the purpose of obtaining a combined cross 
section, that this channel be orthogonal with the other 
top channels with which it overlaps: ee, e/z, and e+jets. 
This is accomplished by vetoing any event that passes 
the selection requirements of any one of these channels. 
As shown in Table 16, four events pass all ev selection re- 
quirements. One of the events has four jets with Et > 15 
GeV, as would be expected for a £±jets event, and the 
remaining three events have only two jets, which is more 
characteristic of dilepton events. 

The background from VK±jets is modeled with VECBOS 
Monte Carlo distributions that are scaled to match the 
jet E T , $ti an d spectra found in data. The Monte 
Carlo sample is normalized to the number of W{— * ev)+2 
jet events found in data and Eq. 6.5 is used to compute 
the expected background of 0.5 ± 0.3 events, as shown in 
Table 17. 

The QCD multijet background estimate is obtained 
from data and is defined as the mean of the results 
from two independent methods. In the first method, 
the probability for a jet to be misidentified as a loose 
electron is determined from a sample of multijet data to 
be (0.0091 ± 0.0012)% in the CC region of the calorime- 
ter. This probability is then applied to the number of 
jets with Et > 20 GeV in a sample of three or more 
jet events where all requirements except that of elec- 
tron identification have been applied. This method re- 
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TABLE 16. Number of observed and expected ev events passing at each cut level of the offline analysis. Expected number of 
tt events are for m t = 172.1GeV/c 2 . Uncertainties correspond to statistical and systematic contributions added in quadrature. 



Number of ev events passing cuts 







Total 


Mis-id 


Physics 






Data 


sig + bkg 


bkg 


bkg 


tt 


le, E C T > 20 GeV, + min e id + trig 


119,263 








71.5 ±20.2 


+ ^j? 1 > 50 GeV 


3941 




434 ± 74 




36.0 ± 10.2 


+ 1 jet, > 30 GeV 


1422 




357 ±61 




35.5 ±10.1 


+ 2 nd jet, E^ > 30 GeV 


192 


244.4 ± 39.0 


92.9 ± 16.0 


121.2 ±35.6 


30.3 ±8.6 


+ M^y > 115 GeV/c 2 


25 


29.3 ±4.8 


24.4 ±4.7 


1.0 ±0.4 


3.9 ± 1.1 


+ A(j)(# T ,2 nd E T object) > 0.5 


12 


18.1 ±3.0 


13.7 ±2.9 


0.9 ±0.4 


3.6 ± 1.0 


+ loose e id 


5 


4.1 ±0.8 


0.69 ±0.12 


0.75 ±0.35 


2.7 ±0.8 


+ orthogonality to other channels 


4 


2.9 ±0.7 


0.47 ±0.15 


0.72 ± 0.34 


1.7 ±0.5 



suits in an estimate of the QCD multijet background 
of 0.576 ± 0.077 (stat) ± 0.076 (sys) events. In the sec- 
ond method, the standard rate for an extra-loose can- 
didate to be misidentified as a loose candidate (see Ta- 
ble 6) is applied to a sample of electron ± jet events 
(extra-loose electron identification) to which all other 
other kinematic cuts have been applied. This method 
results in an estimate of the QCD multijet background 
of 0.367 ±0.129 (stat) ±0.005 (sys) events. The mean of 
these two approaches yields an expected QCD multijet 
background of 0.47 ± 0.15 events, as shown in Table 17. 

The backgrounds from WW and WZ events are ob- 
tained via Eq. 6.5 from PYTHIA Monte Carlo normalized 
to the theoretical cross section [87] , and are given in Ta- 
ble 17. 

As shown in Table 16, the cuts on Jff^ and Mjf are 
most effective in reducing the background. This is shown 
in Fig. 13, where $f vs Mjf is plotted for the VF±jets 
and QCD multijet backgrounds (a,b), for tt Monte Carlo 
events (c), and for data (d). It can be seen that the four 
candidate events are well inside the signal region and far 
from the cut boundaries. 

As described in Sec. VI, tt acceptances are computed 
via Eq. 6.4 using Monte Carlo events generated with her- 
wig and passed through the D0 detector simulation (see 
Sec. V). The trigger efficiency is obtained from the Trig- 
ger Simulator (see Sec. V) and found to be 99.4+g®%. 
The final acceptances for seven top quark masses (and 
for all channels) are given in Sec. X. The expected num- 
bers of tt events passing this selection are determined via 
Eq. 6.3 and are given in Table 17 for these same seven 
masses. Finally, a cross section of 9.1 ±7.2 pb is obtained 
for the ev channel. 

To test the robustness of the background predictions, 
a comparison is made between the data and expectations 
in regions dominated by background (i.e., at earlier steps 
along the selection chain). Making use of Eqs. 6.3 - 6.5 
for the different stages of the selection, Table 16 shows 
that the expectation from background and tt compare 
well with what is observed in the data at the various 
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FIG. 13. Scatter plots of $ T vs for the ev channel: 

(a) H / ±jets background, (b) background from multijet events 
with a misidentified electron, (c) tt signal (mt = 170 GeV/c 2 ), 
and (d) data. The signal region is shown in the upper right 
corner of each plot (Mjf > 115 GeV/c 2 , $ T > 50 GeV). 
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TABLE 17. Expected number of ev signal and background 
events after all cuts in 112.3 pb _1 . Uncertainties are statis- 
tical and systematic contributions added in quadrature. The 
systematic uncertainty on the total background includes cor- 
relations among the background sources. 



Expected number of ev events in 112.3 pb 



« 10 



it MC m (GeV/c 2 ) 




140 


2.96 ±0.88 


150 


2.64 ± 0.77 


160 


2.06 ±0.60 


170 


1.72 ±0.50 


180 


1.49 ±0.43 


190 


1.15 ±0.33 


200 


0.91 ±0.27 


WW 


0.16 ±0.05 


WZ 


0.017 ±0.005 


W+jets 


0.54 ±0.32 


QCD multijet 


0.47 ±0.15 


Total background 


1.19 ±0.38 



stages of the selection procedure. 



VII. ANALYSIS OF LEPTON + JETS EVENTS 

As discussed in Sec. I, the lepton±jets signatures are 
characterized by one isolated, high-pT charged lepton, 
and four or more jets. This signature is similar 
to that of W+jets production. Figures 4 and 6 include 
Monte Carlo distributions for the lepton and jet Et/pt 
and \r)\, and Ifix expected in it lepton±jets events. As 
shown in Table 18, requirements based on these charac- 
teristics form the initial selection for all four channels. 

TABLE 18. Initial selection for ^±jets analyses. The 
^(VF)! cut is introduced and described in Sec. VII A. 



Selection cut 


Topological 
e±jets n+jets 


Muon tag 
e±jets p±jets 


1 isol e, E e T > 20 GeV, 






|?f | < 2.0 + tight e 


id 


yes no 


yes no 


1 isol fi, p 1 ^ > 20 GeV/c, 






W l \ < 1.7(1.0) + ti; 


;ht /j, id 


no yes 


no yes 


H tag 




veto 


yes 


^ al (GeV) 




> 25 > 20 


> 20 


# T (GeV) 




> 20 


> 20 


\r,(W)\ 




< 2.0 




N jcts 




> 4 


> 3 


(GeV) 




> 15 


> 20 


h(jet)| 




< 2.0 


< 2.0 



The triggers used to select the candidate events require 
at least one high-pT lepton and some combination of $ T 
and jet requirements (see Sec. Ill for details). The run 
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FIG. 14. Jet multiplicity distribution for e + $ T +jets data 
(triangle points) and ti Monte Carlo (hatched histogram) af- 
ter initial selection. Trigger inefficiency is not included in the 
Monte Carlo samples. 



ranges and luminosities for the four channels are given in 
Table 19. 

TABLE 19. ^±jets run ranges and luminosities. Channel 
names are as defined in the text. 





e±jets/topo 


/i±jets/topo 


e+jets/ju 


/^±jets//i 


Run range 


la,lb 


la,lb 


la,lb,lc 


la,lb 


Lum. (pb~ ) 


119.5 


107.7 


112.6 


108.0 



The primary background sources are VF±multijet pro- 
duction and QCD multijet events with a misidentihed 
isolated lepton and mismeasured lfi T . As indicated in 
Table 18, the initial selection requires a high-p^ tight lep- 
ton (which dramatically reduces the QCD multijet back- 
ground), large $ T) and several jets. 

Figure 14 shows the number of events as a function of 
the number of jets in the event for e±jets inclusive data 
and for ti MC events after the initial selection. As can 
be seen, the signal to background ratio is still very low. 
It is, therefore, necessary to further exploit the differ- 
ences between signal and background. The most obvious 
differences are in the event topology and the presence 
or absence of a b quark jet. The b quark is inferred 
in the D0 detector by the presence of a non-isolated 
muon (muon tag). Therefore, two orthogonal analyses 
are employed beyond this point: (1) a purely topolog- 
ical analysis, which by construction does not contain a 
muon tag, and (2) an analysis that relies primarily on the 
presence of a muon tag, but also makes use of some topo- 
logical cuts. These channels are denoted respectively as 
^±jets/topo and £+jets//x. The initial selection for these 
channels is given in Table 18. 
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for the /z/u channel (see Eq. 6.2), Ht is defined for the 
lcpton+jets channels as 



2.5 5 7.5 10 2.5 5 7.5 

Expected Background (events) 



10 



FIG. 15. Results of the random grid search in terms of 
expected signal vs expected background for the e+jets topo- 
logical analysis for four possible variable sets: (a) A(W +)ets) 
and h, (b) A(W+ jets), h, and $ T , (c) A(W + jets), H T , and 
tfr, (d) A(W +jets), Ht, andpriW). See text for definitions 
of these variables. 



In order to obtain the most precise measurement of the 
ti production cross section possible, an optimization was 
performed to find those topological variables that pro- 
vide the best separation between signal and background. 
This was accomplished through the use of a random grid 
search [89] in which many possible cut points were tested 
on the signal and background models. Many variables 

jets; 



were investigated in this way: pt(W) = \pr\, Et, ^Yj 



i?T (lepto 



'H£T , two types of aplanarity (A) , and two 



^ — H T Qeta)+pr(W) ' 

types of Ht- Aplanarity is essentially a measure of the 
"flatness" of an event and is defined to be | of the small- 
est eigenvalue of the normalized laboratory momentum 
tensor (M), where this tensor is defined by [90] 



where p is the three momentum of object o, i, j corre- 
spond to the x,y, and z coordinates, and the objects in- 
cluded in the sum depend on the type of aplanarity under 
consideration: (1) only the jets, „4(jets), and (2) the jets 
and the reconstructed leptonic W, A(W + jets). Large 
values of A are indicative of spherical events, whereas 
small values correspond to more planar events. Events 
due to ti production arc quite symmetric as is typical for 
the decay of a heavy object. W+jet and QCD multijet 
events are more planar, owing primarily to the fact that 
the jets in these events arise from gluon radiation. 

Analogous to the transverse-energy variable defined for 
the ee and en. channels, and identical in form to that used 



H 7 



5> 

jots 



(7.2) 



The sum is over all jets with E T > 15 GeV and \rj\ < 2.0 
(recall that the p^i channel uses \r}\ < 2.5). The sec- 
ond transverse-energy variable is simply the sum of the 
standard Ht and the magnitude of the W boson trans- 
verse momentum vector, f/-r(all) = Ht +Pt(W). Events 
due to ti production tend to have much higher values of 
Ht than background. This is due to the fact that the 
jet Et is typically much harder for jets originating from 
the decay of a heavy object than are those from gluon 
radiation. 

The ti sample used in the optimization of all four 
channels is generated using herwig with m t = 180 
GeV/c 2 . The appropriate combination of VF+jets and 
QCD multijet events is used for background. The 
Z — > \m background to the ^+jets/^ channel is not in- 
cluded in the optimization. For the £+jets/(i channels, 
both the VF+jets and QCD multijet background esti- 
mates are based entirely on data. For the topological 
channels, the QCD multijet background is based on data 
and the VF+jets contribution is modeled using the VEC- 
BOS Monte Carlo. These background samples are used to 
investigate the region of phase space remaining after the 
initial selection (see Table 18), and thus differ somewhat 
from the samples used in the full background determina- 
tion to be discussed in Sees. VII A and VII B. 

All of these variables are studied in pairs and 
in different combinations, and for each set of 
cut points a corresponding point in the expected 
(5(signal),B(background)) plane is found. When all such 
points are plotted, they define a boundary that maxi- 
mizes the expected signal for a given background level, 
which is termed the "optimal boundary" (see, for exam- 
ple, Fig. 15). Comparison of the optimal boundaries for 
the various combinations of variables shows that the pair 
A(W + jets) and Ht provides the best signal to back- 
ground ratio for a given signal efficiency. 

After determining that A(W + jets) and Ht are the 
best variables, it is necessary to select which cut point (on 
the optimal boundary) results in the most precise cross 
section measurement. Contours of constant uncertainty 
on the measured cross section (da /a) can be derived from 
the relation 



a = 



e-C 



N-B 
e-C 



(7.3) 



where iV, S, and B are the number of observed, expected 
signal, and expected background events, respectively, e is 
the signal efficiency, C is the integrated luminosity, and 
a is the measured cross section [91]. The cut points on 
the optimal boundary with the smallest da /a and best 
significance (s/b) are (see Fig. 16): 

£+jets/topo: 
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FIG. 16. Expected signal vs expected background plots 
for A(W + jets) and Ht optimization variables for (a) 
/i+jets/topo and (b) e+jets//x. The solid curves are contours 
of constant uncertainty on the cross section (da /a). Arrows 
indicate chosen cut points. 

• H T > 180 GeV 

• ^(W+jets) > 0.065 
£+jets/yu: 

• H T > 110 GeV 

• A(W+jets) > 0.040. 

Following the initial selection and optimization it is 
necessary to make several additional channel-specific re- 
quirements. These requirements, along with the results 
and expectations from signal and background, are dis- 
cussed in the next two sections (VII A and VII B). 

Acceptances for all four £+jets channels are computed 
from Monte Carlo events generated by the herwig [54] 
program for 24 top quark mass values (m t — 90-230 
GeV/c 2 ) and then passed through the full D0 detec- 
tor simulation (see Sec. V). The expected number of 
ti events passing the selection for a given channel is 



N = a t t(m t ) Mh3^ t ) ■ &J (7.4) 

i— runs j—dct 

where a t i(mt) is the theoretical ti cross section at a top 
quark mass of m t [45]; £jj is the integrated luminosity 
for run i and detector region j (CC and EC for electrons, 
CF and EF for muons); and the acceptance is 

A(i, j, m t ) = £ t rig • £ P id • £ sc \ G- B, (7.5) 

where etn g (i,j,mt) is the trigger efficiency, e p id{i,j) is 
the efficiency for lepton identification (isolated leptons 
and muon tag), e so i(i, j, mt) is the efficiency of the selec- 
tion cuts, G(i,j) is the geometrical acceptance, and B is 
the branching fraction for the sample in question. Trig- 
ger efficiencies are obtained from data or Monte Carlo 
events, depending on the channel, and are discussed in 
more detail below. Particle identification efficiencies are 
obtained from data for the case of electrons (as discussed 
in Sec. IV A) and from a combination of data and Monte 
Carlo in the case of muons (as discussed in Sec. IV B). 
The selection efficiencies e so i and the geometrical accep- 
tances G are obtained from Monte Carlo events. As dis- 
cussed in Sec. X, the acceptance, rather than the ex- 
pected number of ti events, is used in the calculation of 
the ti cross section. Typical values for the acceptance, 
often denoted as the "efficiency times branching fraction" 
(ex B), for all eight leptonic channels, are given in Sec X 
for seven top quark masses. The numbers of ti events ex- 
pected in the four £+jets channels are given in Tables 22 
and 25 Sees VII A and VII B for the same set of top quark 
masses. The systematic uncertainties on the acceptances 
and backgrounds are discussed in Sec. IX. 

A. Topological tag 

As described in the previous section, the first two 
stages of the £+jets/topo selection require the cuts de- 
scribed in Table 18 followed by the cuts on A(W + jets) 
and Ht- There is, however, one cut in Table 18 which has 
not yet been discussed. This cut on rj(W), the pseudora- 
pidity of the lepton and ~$ T fit to a W boson hypothesis, 
is designed to remove from consideration those regions of 
phase space where the ly+jets vecbos Monte Carlo does 
not model the W+jets data very well. As can be seen 
in Fig. 17, the VECBOS prediction is considerably below 
the data in the forward region [92] . Therefore, the initial 
selection requires that |»j(W)| < 2.0. It should be noted 
that only a few percent of ti events have |?j(W)| > 2.0, 
so this cut does not represent a serious reduction in ac- 
ceptance. It should further be noted that these analyses 
determine the W+jets backgrounds primarily from the 
data. The VECBOS Monte Carlo is only used to determine 
the survival probability for the cuts on «4(W+jets), Ht, 
and Ej, which is the scalar sum of the lepton Et and 
$ T - As can be seen in Fig. 18, a requirement of E^, > 60 
GeV provides significant rejection against QCD multijet 
background while having little effect on the ti signal. 
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FIG. 17. \il(W)\ distribution for ^+jets/topo data (his- 
togram) for the sum of predicted signal and background (filled 
circles), and background alone (open triangles), after appli- 
cation of all selection criteria except the ij(W) cut. 



As noted above, the primary backgrounds to the 
£+jets/topo channels are from W(— ► £f)+jets and QCD 
multijet events which contain a misidcntihcd electron or 
isolated muon and mismeasured $ T . The mismeasured 
ifjrp arises primarily from mismeasurement of jet Et or 
vertex z position. 

The background calculation proceeds in four steps: 



1. The QCD multijet background is determined as a 
function of the inclusive jet multiplicity from data 
samples in which the A, Ht,t](W), and E^ cuts 
have not been applied. Due to the different pro- 
cesses that give rise to a misidcntified electron or 
isolated muon, these backgrounds are handled dif- 
ferently: 

• Jets that have a large electromagnetic frac- 
tion can sometimes pass the electron identi- 
fication criteria and be misidcntified as elec- 
trons. To determine the background from 
multijet events containing such misidcntified 
electrons and $ T , one begins with the $ T 
spectrum from n + 1 jet (n > 0) events with 
|f?(W0l < 2.0 which pass an electron trigger 
but fail the full electron identification cuts 
(mis- id e + $ T sample). This sample cor- 
rectly describes (with sufficient statistics) the 
$ T distribution for the QCD multijet back- 
ground, but the normalization is not correct 
since the electron identification requirement 
has not been made. The correct normalization 
is obtained by matching the number of events 
at low $ T ($ T < 10 GeV) to that found in 
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FIG. 18. Elf distributions for tt Monte Carlo (m t = 170 
GeV/c 2 ) (dashed histogram), and for QCD multijet data 
(solid histogram), after application of all selection criteria ex- 
cept those on E^, A, and Ht- The distribution for IT+jets 
is similar to that for tt. The solid vertical line at Eff = 60 
GeV indicates the cutoff value. 



a complementary sample that passes the nor- 
mal electron identification criteria. Requir- 
ing $ T > 25 GeV then provides the expected 
number of QCD multijet background events to 
the e + n jet selection. Uncertainties on this 
procedure are dominated by the statistics of 
the samples used and range from 9.5% (13%) 
for the Run la (Run lb) e + 1 jet selection to 
27% (54%) for the Run la (Run lb) e + 4 jet 
selection. 

Muons from the semi-leptonic decay of a b 
or c quark are normally accompanied by an 
associated jet (non-isolated). However, occa- 
sionally the decay kinematics are such that 
there is insufficient hadronic energy to pro- 
duce a jet. In these cases the muons from 
semi-leptonic b and c decays will appear to be 
isolated. The probability that a muon origi- 
nating from the decay of a heavy quark will 
appear isolated varies with jet multiplicity, 
run period, and detector region, and is de- 
noted by 7 m i s _id (run, det). Typical CF values 
are 11% for /i+ > 1 jet events and 6% for 
fi+ > 2, > 3, > 4 jet events (the correspond- 
ing EF values are 22% and 15% respectively). 
For a given jet multiplicity, n, these probabili- 
ties are measured using samples of QCD mul- 
tijet events with $ T < 20 GeV as the ratio 
of the number of isolated-^ + > n jet events 
to the number of non-isolated- fi + (> n + 1) 
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jet events. The QCD multijet background is 
defined by the product of this probability and 
the number of non-isolated- \i + (> n + 1) jet 
events with $ T > 20 GeV. The primary un- 
certainty in this method stems from the de- 
termination of the above misidentified muon 
isolation probabilities. The value of 30% as- 
signed to this uncertainty is dominated by the 
statistical precision of the control sample used 
to derive the false isolation fraction for four-jet 
events. 

These procedures are carried out for each inclu- 
sive jet multiplicity, thereby providing the expected 
QCD multijet contribution to the £+ > n jet selec- 
tions (n = 1,2,3,4), as defined in Table 18. For the 
£+ > 4 jet selection, the expectation is 4.4 ± 2.2 
events in the e+jets/topo channel and 6.44 ± 2.08 
in the ^+jets/?rra channel. 

2. The background from W(— ► £v)+jets is computed 
by performing a fit to the jet-multiplicity spectrum 
that remains following the subtraction of the QCD 
multijet background. Inherent in the fit is the as- 
sumption of "Berends (iVjets) scaling" [93,94] which 
suggests that there is a simple exponential relation- 
ship between the number of events and the jet mul- 
tiplicity: 



<j(W + n jets) 
<r{W + (n- 1) jets] 



= a, 



(7.6) 



where a is a constant (for any given jet Et and r\ 
requirements) and n is the inclusive jet multiplic- 
ity. For any given inclusive jet multiplicity i, the 
number of events which are observed following the 
QCD multijet subtraction is given by 



N° hs = NY ■ a 1 ' 1 + /f p • N top , 



(7.7) 



where is the number of W + 1 jet events, N top 
is the number of ti events in the sample, and /* op 
is the fraction of ti events with jet multiplicity i 
(obtained from Monte Carlo). The values of N° hs 
are plotted in Fig. 19. Fits to Eq. 7.7 determine 
the values of a given in column 2 of Table 21 (N^ 
and 7V top are also obtained from this fit). Once 
a is known, the number of W + 4 jet events that 
pass the initial selection can be determined from 
the equation 



N w = N w 



a 



(7.8) 



The resulting W+jets background after the £+4 jet 
selection is 37.2 ± 4.5 events for the e+jets channel 
and 18.8 ± 3.2 events for the ju+jets, as indicated 
in Table 21. This method, solely based on data, 
is independent of theoretical calculations of W + n 
jet cross sections which have large uncertainties at 
high jet multiplicities. 
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FIG. 19. Number of events as a function of inclusive jet 
multiplicity for the e+jets/topo and /i+jets/topo analyses. 
All cuts have been applied except A and Ht- The linear 
nature of the distributions is known as Berends scaling. Note 
that since the cut has been applied, the values here differ 
from those in Table 20. 

3. For the e+jets channel only, a correction factor 
of 1.09 ± 0.39 (1.71 ± 0.12) is applied to the Run 
la (Run lb) QCD multijet background results to 
account for trigger differences between the back- 
ground method and the actual data selection and 
for the increased luminosity from the inclusion of 
the Main-Ring data (see Sec. Ill and Appendix C) 
in the Run la and Run lb data sets. A similar cor- 
rection factor of 1.09±0.17 (1.22±0.06) is applied to 
the Run la (Run lb) IT+jets background. Follow- 
ing these corrections, the backgrounds to the e + 4 
jets selection are found to be 7.2 ± 2.2 events from 
QCD multijet and 44.8 ± 8.6 events from VF+jets. 

4. To determine the expected background following 
the final three cuts on E^., A, and Ht (see Ta- 
ble 21), a cut survival probability f is computed for 
each background. This probability factor is applied 
to the results obtained after the £+ > 4 jet selec- 
tions, thus giving the final expected QCD multijet 
and VF+jet backgrounds: 

A r (total bkg) = N?™ ■ /qcd + NJl 4j ■ f w (7.9) 

where N^P and A^ 4j are the QCD multijet and 
W+jet background estimates following the £+ > 
4 jet selections, and /qcd and fw are the sur- 
vival probability factors for the QCD multijet and 
VF+jets backgrounds respectively. 

• For the e+jets channel, /qcd is determined 
from the combined E^, A, and Ht pass rate 
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TABLE 21. Steps in e+jets/topo and p+jets/topo background calculation: column 2, row 1 gives the expected number of 
QCD multijet background events (£+4 jets); column 1, row 2 gives the value of a determined from the fit to Eq. 7.7; column 
2, row 2 gives the expected number of VT±4 jet events; column 3 gives the trigger and Main-Ring (MR) correction factors; 
column 4 gives the result of multiplying column 2 by column 3 (step 3 in the text); column 5 gives the E^f, A, Ht cut survival 
probabilities; and column 6 gives the final expected background obtained by multiplying column 4 by column 5. Note that 
Runs la and lb are treated separately for the e+jets channel whereas they are treated as a single run for the /^+jets channel. 





a 


Exp # of evts 
Steps 1-2 


Trigger & 
MR. corr 


Exp # of evts 
Step 3 


By , A, H T cut 
survival prob. (/) 


Exp # evts 
Step 4 


e+jcts QCD multijet la 
lb 
Total 




0.7 ± 0.8 
3.7 ± 2.0 
4.4 ± 2.2 


1.09 ± 0.39 
1.71 ± 0.12 


0.76 ± 0.91 

6.4 ± 2.0 
7.16 ± 2.20 


0.071 ± 0.040 
0.051 ± 0.010 


0.054 ± 0.072 
0.325 ± 0.119 
0.379 ± 0.139 


W+jcts la 
lb 
Total 


0.17 ± 0.02 
0.18 ± 0.01 


5.45 ± 1.53 
31.77 ± 4.24 
37.21 ± 4.50 


1.09 ± 0.17 
1.22 ± 0.06 


5.9 ± 1.9 
38.9 ± 8.3 
44.8 ± 8.6 


0.092 ± 0.061 
0.092 ± 0.061 


0.544 ± 0.185 
3.590 ± 0.799 
4.135 ± 0.899 


/^+jcts QCD multijet 




6.44 ± 2.08 




13.9 ± 4.4 




0.993 ± 0.498 


W+jcts 


0.19 ± 0.02 


18.8 ± 3.2 


1.37 ± 0.07 


25.8 ± 4.6 


0.129 ± 0.027 


3.324 ± 0.911 



TABLE 20. Number of £+jets/topo data events passing at 
each cut level. Note that the e+jets luminosity of 90.9 pb -1 
does not include recovered Main- Ring data (see Appendix C) 
- the Main-Ring contribution is given in parentheses. Simi- 
larly, the luminosity for the /i+jets channel does not include 
Run la or recovered Main-Ring data. The Main-Ring contri- 
bution plus that from Run la is given in parentheses. 



TABLE 22. Observed and expected number of £+jets/topo 
signal and background events after all cuts. Uncertainties 
shown are statistical and systematic contributions added in 
quadrature. The total background systematic uncertainty in- 
cludes correlations among the different background sources. 





e+jets 


/i+jets 


Lum (pb L ) 


90.9 


76.6 


N jets > 1 


6604 


2127 


iVjets > 2 


1225 


537 


iVjcts > 3 


223 


124 


iVjets > 4 


39 


28 


Eff > 60 GeV, 


39 


22 


A > 0.065, 


18 


10 


H T > 180 GeV 


7(2) 


4(6) 



on a sample of misidentified electron+4 jet 
events that satisfy the $ T and rj(W) require- 
ments. 

• For the /x+jets channel, the prescription is 
simply an extension of the QCD multijet back- 
ground computation described above for the 
/i + n jet selection. Specifically, the selection 
criteria are applied to five-jet events, where 
the jet associated with the non-isolated muon 
is not included in the A and Ht calculations. 

For both channels, fw is determined using the vec- 
BOS Monte Carlo program to measure the final ef- 
ficiency (including the 1+ > 4 jet, Ej,, A, and 
Ht cuts) relative to that for the £+ > 4 jet se- 
lection. To investigate the systematic uncertain- 
ties associated with this Monte Carlo based pro- 
cedure, samples are generated with two different 
Q 2 scales, M w and (p^(jet)), and with two differ- 
ent hadronic fragmentation prescriptions, ISAJET 
and HERWIG. Comparison with the background- 





e+jets 


/i+jets 


Lum (pb 1 ) 


119.5 


107.7 


Observed 


9 


10 


top MC m t (GeV/c 2 ) 






140 


12.06 ±5.20 


8.22 ±3.56 


150 


11.20 ±3.72 


7.83 ±2.98 


160 


10.11 ±2.35 


7.12 ± 2.40 


170 


8.97 ± 1.61 


5.72 ± 1.72 


180 


7.44 ± 1.04 


4.80 ± 1.27 


190 


5.70 ±0.68 


3.84 ±0.92 


200 


4.60 ± 0.47 


3.14 ±0.69 


W+]ets 


4.14 ±0.90 


3.32 ±0.91 


QCD multijet 


0.38 ±0.14 


0.99 ±0.50 


Total background 


4.51 ±0.91 


4.32 ± 1.04 



enriched sample of data indicates that vecbos gen- 
erated at Q 2 = (p^(jet)) and fragmented through 
herwig provides the best match. This choice is 
therefore used to compute the values of fw- 

These four steps are summarized in Table 21. 

Figure 20 shows the distribution of A vs Ht for £±jets 
(combined e±jets and /x±jets) events for data, herwig 
ti Monte Carlo (m t = 170 GeV/c 2 ), QCD multijet, and 
VECBOS I4 y ±jets Monte Carlo events. From this figure it 
is clear that A and Ht provide significant discrimination 
between signal and background. 

As described in Sec. VII, ti acceptances are computed 
via Eq. 7.5 using Monte Carlo events generated with 
herwig and passed through the D0 detector simulation. 
The trigger efficiency for the e±jets channel is obtained 
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FIG. 20. Scatter plots of A vs i/ T for £+jets data (d) 
compared to expectations from higher-luminosity samples of 
it MC (m t = 170 GeV/c 2 ) (c), and QCD multijet (b) and 
VK+4 jet MC (a) backgrounds. The dashed lines represent 
the threshold values used for selection. The effective lumi- 
nosity given for plot (b) is determined as the product of the 
luminosity of the selected multi-jet sample and the inverse of 
the appropriate misidentification rate. 



from W+jets data and determined to be 98.2^4;4%. For 
the zi+jets channel, the trigger efficiency is computed 
using data-derived trigger turn-on curves applied to it 
Monte Carlo and is determined to be 89 ± 5%. The ac- 
ceptance values after all cuts for seven different top quark 
masses (and for all channels) are given in Sec. X. 

Following Eq. 7.4, the expected numbers of it events in 
the f+jets/topo channels are given in Table 22 for these 
same seven masses. Also shown are the final numbers of 
events observed in the data, 9 in the e+jets channel and 
10 in the /z+jets channel. Table 20 shows the observed 
number of data events passing at the different stages of 
the selection procedure. Note that for this table, the 
e+jets luminosity does not include Main-Ring data and 
the zt+jets luminosity does not include Run la or Main- 
Ring data. Finally, the cross section obtained from the 
e+jets/topo and /z+jets/topo channels are 2.8 ± 2.1 pb 
and 5.6 ± 3.7 pb, respectively. 



B. n tag 

The initial selection for ij+jets//x events is described 
in Sec. VII and summarized in Table 18. All events are 
required to have a /x tag as defined in Sec. IV B. 

The dominant backgrounds that remain after the ini- 
tial selection arise from W{— * £^)+jets production, QCD 
multijet events that contain a misidentified electron or 
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FIG. 21. Scatter plots of A0(/z, tf T ) vs # T for (a) e+jets//z 
QCD multijet background, (b) it — > e + jets//i, (c) /i+jets//z 
QCD multijet background, and (d) it — > /z+jets//z. The solid 
lines define the cut boundaries. 

isolated muon and mismeasured $ T , and also Z(— > 
/x/x)+jets for the /i+jets//j channel. 

For events that have no genuine source of#r, the pres- 
ence of a muon, as a consequence of the muon system's 
modest momentum resolution, may lead to mismeasured 
ifjj* which is aligned or anti-aligned with the muon px- 
Indeed, in multijet data, the distribution of the angle (f> 
between the muon momentum and the direction of the 
$ T , A0( / u, $ T ), peaks at 0° and 180°, whereas for it 
events this distribution rises monotonically from 0° to 
180° as indicated in Fig. 21. In order to reduce back- 
ground from QCD multijet events, both /U-tag channels 
make a cut on the allowed region in the $ T , A0(/i, $ T ) 
plane: 

JpT > 35 GeV, if |A</>(/j, $ T )\ < 25°, for e+jets, 

(7.10) 



and, 



A0(/i, $ T ) 
|A<K M ,# T )-90°| 



< 170° 



and 



90° 



< 



45 GeV 



for /j, + jets. (7.11) 



The effectiveness of these cuts is displayed in Fig. 21, 
which shows the distributions in the $! T ,A<j)(jj,, $ T ) plane 
for QCD multijet events and it Monte Carlo events for 
both /z-tag channels. 

In addition to the QCD multijet and H^+jets back- 
grounds noted above, the /x+jets//x channel, by virtue of 
the fact that it requires two muons, has a non-negligible 
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FIG. 22. \ 2 probability distribution for the fi+jets/fj, chan- 
nel after all cuts except P{x 2 )'- ( a ) Z(— * wO+jets MC and 
(b) it MC. 

background from Z(— * /x/z)+jets production. Although 
the muons from Z boson decay are, in principle, isolated, 
there is a small probability that one of them will overlap 
with one of the jets in the event and thus appear to be 
non-isolated. The /z+jets//z channel relies therefore on 
a kinematic fitting procedure to reduce this background. 
As described in Sec. VI C, a kinematic fit to the Z — > fifi 
hypothesis is performed and a x 2 is obtained (see Eqs. 6.7 
and 6.8). Events with a x 2 probability greater than 1%, 
P(x 2 ) > 0.01, are considered likely Z boson candidates 
and are therefore rejected. As can be seen in Fig. 22, 
this procedure provides very good rejection against the 
Z(— ► fi/i)+jets background and has essentially no effect 
on the it signal. 

The general scheme for background calculation pro- 
ceeds in three steps which are first outlined and then 
discussed in detail: 

1. Compute the QCD multijet background: 

• For the e+jets//x channel, the QCD multi- 
jet background is computed by applying an 
electron misidentification rate to a ^-tagged 
multijet control sample passing all cuts except 
tight electron identification. 

• For the /z+jets//z channel, the QCD multijet 
background is computed by applying isolated- 
muon and muon-tag misidentification rates 
to an untagged QCD multijet control sample 
passing all other cuts except the isolated muon 
requirement. 

2. Compute the M^+jets background: 

For both channels, the background from VF+jcts 
events is computed by applying a muon tag rate to 
the number of untagged multijet-subtracted I + > 
3 jet data events and then subtracting the expec- 
tation from it: 

N w {bkg) = N(data - QCD) • P tag - JVf (7.12) 
where 

"data" is the number of events passing all cuts ex- 
cept /i-tag; 



"QCD" is 

(the number of (extra-loose-e+> 3 jet)/(> 4 jet) 
events passing rfr,A, and Ht cuts) ■ (e/n mis-id 
rate); 

"-ftag" is the probability (as a function of jet Et 
and T], and run period) for a jet to contain a tagged 
muon, and is determined from QCD multijet data; 
"iVg*" is the expected top quark contribution after 
all cuts and is computed differently for the e+jets//x 
and /U+jets/jU channels: 

• for the e+jets/n channel, the expected top 
quark contribution is determined from data 
by fitting the jet spectra of the multijet- 
subtracted untagged e + n jet data under the 
assumption of jet scaling and measuring the 
excess for n > 3. A tag rate derived from it 
MC is applied to this excess to obtain TV". 

• for the n+jets/f/, channel, the expected top 
quark contribution (JV| ) is determined from 
herwig MC normalized to the theoretical 
cross section [95]. 

3. For the /z+jets//z channel only, determine the back- 
ground from Z — > /!/! using VECBOS MC events. 

The key elements of this procedure, namely the QCD 
multijet background calculations and the parameteriza- 
tion of the muon-tagging probability, are motivated and 
developed below. 

The estimation of the multijet background differs 
somewhat in the e+jets//i and /u+jets//u channels. The 
calculation for the e+jets/^ channel is similar to that 
used for the ^+jets/topo channels. Namely, the QCD 
multijet background is determined by relaxing the elec- 
tron identification criteria and observing the number of 
additional events that pass the selection. It is assumed 
that the number of events in the extra- loose electron sam- 
ple, Ni, consists of both real, N e , and misidentified (often 
referred to as "fake"), Nf, electrons 

Ni = N e + N f . (7.13) 

The probability for a real electron to pass from the loose 

sample into the tight sample, e\, is determined from 

Z — > ee data. Similarly, the probability for a misiden- 

f 

tified electron to make this transition, e t , is defined as 
the ratio of tight to loose electron events in a sample 
of "loose electron+1 jet" events without I^t [69]. These 
probabilities are determined separately for the CC and 
EC regions of the calorimeter and are given in Table 23. 
Applying these probabilities to the number of real and 
misidentified electrons in the loose sample gives the ex- 
pected number of events in the tight sample: 

N t = e e t N e + e{N f . (7.14) 

Eqs. 7.13 and 7.14 can be solved for the number of 
misidentified electron events in the loose sample: 
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N 



f 



s e t N t 



£* - si 



(7.15) 



The expected number of misidcntificd electron events in 
the final sample is the product of the number in the loose 
sample and the probability for a misidentified electron to 
pass the tight requirement, s\Nf. Values for the CC and 
EC regions of the calorimeter are given in Table 23. The 
combined (CC+EC) QCD multijet background for the 
e+jets//i channel, including additional systematic uncer- 
tainties (see Sec. IX) not given in Table 23, is tabulated 
later in this section. 

TABLE 23. e+jets//u QCD multijet background calcula- 
tion parameters 





CC 


EC 


N t 


4 


1 




8 


6 


et 


0.828 ±0.010 


0.453 ±0.015 




0.027 ±0.009 


0.053 ±0.012 


N f 


3.28 ±0.11 


4.30 ±0.31 


4-Nf 


0.088 ± 0.030 


0.228 ±0.054 



0.02 
0.015- 
0.01- 



bo 0.005 
bp 



a 

o 



(a) 




20 



40 



60 80 100 

JetE T (GeV) 



120 



(b) 




0.01 
0.0075 

0.005 - 
0.0025 


1 2 

JetT) 

FIG. 23. Parameterization of the muon tag rate, for muons 
in the CF region from Run lb, as a function of (a) jet Et and 
(b) jet r). 




The calculation of the QCD multijet background 
for the /U+jets//« channel is an extension of that 
used for the ^t±jets/topo channel. As described in 
Sec. VII A, the QCD multijet background calculation for 
the ^±jets/topo analysis applied the probability for a 
muon from a 6 or c quark decay to appear isolated to the 
number of non-isolated-/Lt±jet events to determine the ex- 
pected number of misidentified isolated muon events in 
the signal sample. The /z+jets//i analysis extends this by 
applying an additional tag rate function. This tag rate 
function is based on a Monte Carlo sample containing 
a high fraction of 6-quark jets, and is parameterized in 
terms of the jet E T as: 

, n r-w , n , ( e t - 15.0 GeV\ 
h(E T , run, dct) = £>(run, det) • tanh ( — 49 q Q c y — j 



(7.16) 

where D(run, det) is a scale factor that depends on the 
run period and detector region under consideration. The 
QCD multijet background to the [i+jets/n channel is 
then determined from the product 

Nqcd = ^2 ^2 N ° ' ^mis-id (run, dct) • h(E T , run, dct), 

run, dct jets 

(7.17) 

where N is the number of events which pass all selection 
criteria except for the isolation requirement on the high- 
Pt \i and the /x-tag requirement, and / m j s _id(run, det) 
is the misidcntificd-isolated-^. probability discussed in 
Sec. VII A. The final value, including systematic uncer- 
tainties, is tabulated at the end of this section. 



The jets produced in association with W boson pro- 
duction originate primarily from final state gluon radi- 
ation. Therefore, except for a small contribution from 
gluon splitting (g — ► bb), W+jets events are expected 
to contain very few b quarks and thus very few muon 
tags. In order to estimate this background, it is assumed 
that the heavy flavor (b and c quark) content in FV±jets 
events is the same as in QCD multijet events [58]. The 
expected number of VF±jets±/x tag events is therefore 
computed from the product of the number of untagged 
W+jet events and a muon-tag probability (Ptag) 



N w (tagged) = N w (not tagged) • P t 



tag- 



(7.18) 



This probability is defined in a control sample of mul- 
tijet events by the fraction of jets that contain a muon 
within a cone of AR — 0.5 around a jet axis. The con- 
trol sample consists of events collected with a multijet 
trigger (jet-MULTI, see Table 5) that have four or more 
jets reconstructed offline (Et > 15 GeV, \rf\ < 2). These 
events were collected under essentially the same detec- 
tor and accelerator conditions as the signal sample. The 
multijet and untagged W / ±jets samples have similar jet 
Et and r\ distributions, and, since both samples owe their 
high jet multiplicity to gluon radiation, they should also 
have similar quark- flavor content. 

This fraction, also known as the tag rate, is param- 
eterized explicitly as a function of jet Et and 77, and 
is handled separately for the CF and EF regions of the 
muon system. The r\ dependence is fit independently for 
the different run intervals used in the two analyses (see 
Table 19). The tag rate as a function of jet Et and X] for 
muons in the CF region for Run lb is shown in Fig. 23. 
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The tag rate increases with jet Et because higher-energy 
jets have, on average, higher energy muons that are more 
likely to penetrate the calorimeter and magnet and be de- 
tected. The shape of the 77 distribution is primarily due 
to the geometrical acceptance of the muon system, but 
varies somewhat over the different run intervals. As a 
function of jet Et, the data are fit to the functional form 

f ( Tp \ J A 1 +A 2 E T + A 3 E T , for E T < p 
f{ET) -\A l+ A 2 p + A,p\ for£ T >p (7 ' 19) 

where p = — \A 2 jA 3 , and the parameters A\, A 2 and A3 
are free. The resulting curves for muons in the CF and 
EF regions are denoted / CF and / EF respectively. As a 
function of 77, the data for muons in the CF region are fit 
to the functional form 

g CF { V , r) = B hr (1 + B^if) [ erf ( V B 2 , r + B 3 , r ) (7.20) 

- eri(r]B 2 . r - B 3 , r )], 

where r labels the three periods of the run as specified in 
Sec. IV B, erf(a;) = 2/y/n J* exp(— t 2 ) dt, and the param- 
eters Bi ri B 2 r , B 3 ^ r , and B^ r are free to vary. Similarly, 
for muons in the EF region, the data are fit to 

g FF (rj) = & { erf [(M - C 4 )C 2 + C 3 ] (7.21) 
-crf[(M-C 4 )C2-C 3 ]}, 

with free parameters C\, C 2 , C 3 , and C 4 . There is no 
run dependence in Eq. 7.21, since, as noted in Sec. IV B, 
the EF region of the muon system was only used during 
the final run period (Run lb+c postclean). The complete 
tag rate function is 

Pt, s (E T ,V,r) = D? F f CF (E T )g CF (v,r) (7.22) 
+ £ ,EF /E F ( £ T)5 EF (r?); 

where D^ F and Df F are constants that normalize the 
predicted number of tagged jets in the control sample 
to the actual number. The values of the parameters in 
Eqs. 7.19-7.22 are given in Table 24. 

The accuracy of this procedure has been studied by 
comparing the predicted to observed number of events 
having a tagged jet for a variety of data samples 
representing different trigger conditions, physics pro- 
cesses, and jet multiplicities. These studies are sum- 
marized in Fig. 24, which shows the (Observed - Pre- 
dicted) /Predicted values for data samples that originate 
from nine different triggers (see Table 5 for the definitions 
of these triggers): 

• The inclusive multijet samples with minimum jet 
multiplicity of two, three, four, and five were taken 
with the triggers jet-min, jet-3-mon, jet-4-mon, 
and JET-MULTI, respectively. The last sample, with 
five jets selected offline, is a complete subset of the 
four jet sample used in the actual tag rate calcula- 
tion, comprising about one-third of the jets in the 
control sample. 



>2jets 

■ >3 jets 

>4jets 
■ >5 jets 
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■ e+>2jets 
— ■ — y+>l jet 
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FIG. 24. Tests of the muon tag rate. Shown are (Ob- 
served - Predicted)/Predicted values for data sets that orig- 
inate from nine different triggers. Some of the scatter is due 
to statistics, as indicated by the horizontal error bars; the 
remainder is ascribed to systematic effects as described in 
Sec. IX A 16. The solid vertical line is the overall mean value 
and the dashed vertical lines are the uncertainty on the overall 
mean. 

• The electron samples consist of events with a 
tight electron candidate, taken with the ele-1-mon 
(gis-DUEt) trigger for the case of one (two) or more 
additional jets. Almost all of the "electrons" are 
false. The purpose of examining these events is to 
check for an excess of tags due to bb or cc produc- 
tion, where one heavy quark decays to an electron 
and the other to a muon. There is no evidence of 
such an excess, and none is expected because of the 
isolation and high Et requirements imposed on the 
electron. 

• The photon samples consist of events with a tight 
photon candidate (see Sec. IV A), taken with the 
same triggers as the electron samples. About 30% 
of the 7+ > 1 jet events arc from direct-photon pro- 
duction and the rest are from multijet background 
[96]. The purity is less in the 7+ > 2 jet data. 

• The Z + jet data were obtained with the EMl- 
eistrkcc-esc trigger, by requiring two loose elec- 
tron candidates including at least one tight candi- 
date. The invariant mass of the electron pair is 
required to be between 80 and 100 GeV/c 2 . The 
background in this sample is low (10%); but unfor- 
tunately only four events with a tagged jet survive, 
so the statistical uncertainty is quite large. 
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TABLE 24. e+jets+^ tag parameters from Eqs. 7.19 - 7.22. 



/ Parameters 




CF 

9 


Parameters 




g EF Param. 


Normalization Param. 




CF value 


EF value 




r = l 


r = 2 


r = 3 


r = 3 


r = l r = 2 r = 3 


Ai 


-0.243E-2 


-0.902E-3 


Bl }T 


0.386E-2 


0.363E-2 


0.395E-2 


Ci 


0.349E-2 




249.6 248.7 223.4 


A 2 


0.170E-3 


0.84 7E-4 


B2,r 


11.5 


2.26 


4.78 


c 2 


3.92 


D r EF 


528.8 


A 3 


-0.397E-6 


-0.368E-6 




12.4 


2.17 


4.85 


c 3 


1.54 












B4 >r 


-0.483 


-0.477 


-0.198 


c 4 


1.43 







The horizontal error bars shown in Fig. 24 reflect the sta- 
tistical uncertainty on each comparison. As discussed in 
Sec. IX, that portion of the scatter that cannot be at- 
tributed to the statistical uncertainty is taken as a mea- 
sure of the systematic uncertainty of the tag rate proce- 
dure. 

The functional dependence of the tag rate is important 
only to the extent that the target sample differs from the 
control sample. It should therefore be noted that the 
test samples with low jet multiplicity have significantly 
steeper jet Et spectra than either the control sample or 
the W^+jets data after application of the A and Ht cuts. 

Because these analyses are concerned with the number 
of tagged events that remain in a data sample following 
selection cuts on Ht and A, it is important to confirm 
that the tag rate does not depend on these variables in 
an unexpected way. Figure 25 shows a comparison of the 
predicted and observed numbers of tagged events as a 
function of Ht and A for the > 3 jet and > 4 jet test 
samples. The aplanarity distributions are in good agree- 
ment. Differences in the Ht distributions suggest that 
a cut could result in a discrepancy of a few percent be- 
tween the predicted and observed number of events. This 
is among the contributors to the tag-rate uncertainty that 
are discussed in Sec. IX A 16. 

As noted in the outline at the beginning of this sec- 
tion, contamination from QCD multijet and ti events 
requires that the background from VF+jets be com- 
puted via Eq. 7.12. The QCD multijet contribution to 
the untagged sample is estimated by applying the lep- 
ton (e//i) misidcntification rate to a sample of (loose- 
e+ > 3 jet)/(> 4 jet) events that have passed the $t, A 1 
and Ht requirements. The ti contribution (iV") for the 
e+jets//i channel is determined from data by fitting the 
jet spectra of the QCD-multijet-subtracted e + n jet data 
under the assumption of jet scaling and measuring the 
excess for n > 3. Following the hypothesis of jet multi- 
plicity scaling, the number of W+jet events can be de- 
scribed by a function of the form 



n, a 



if M fa 



(7.23) 



where rii is the number of events with i or more jets, ri^ 
is the number of W boson events with three or more jets, 
ft is the number of events in the ti MC sample with i or 
more jets, and a is a free parameter. A fit to Eq. 7.23 



finds 



to be 19.2 ± 9.5 events. 
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FIG. 25. Predicted (histogram) and observed (filled circles) 
Ht and A distributions in multijet data: (a) Ht distributions 
for > 3 jet data, (b) A distributions for > 3 jet data, (c) Ht 
distributions for > 4 jet data, and (d) A distributions for > 4 
jet data. 



(Nf) for the \x- 



applying the tt tag rate (P t " g ) to n\ . The tt contribution 



Hjets/^i channel is determined from her- 
WIG MC normalized to the theoretical cross section [95]. 

As given in Table 25, the W+jets backgrounds for 
the e+jets/ii and /x+jets//x channels are determined via 
the multi-step procedure above to be 0.74 ± 0.30 and 
0.73 ± 0.14 events respectively. Systematic uncertainties 
on the VF+jets background arise primarily from uncer- 
tainties in Berends scaling and ti MC tag rate (e+jets//x 
channel only) and the tag-rate parameterization. These 
are discussed in Sec. IX. 

The background from Z — > /ifi to the /U+jets//x channel 
is determined from vecbos Z+jets Monte Carlo events 
in a fashion similar to the Monte Carlo background cal- 
culations used for the dilepton channels (see Eq. 6.5) and 
is given in Table 25. 

Backgrounds from single top, WW, and WZ produc- 
tion were also studied and found to have a negligible con- 
tribution to the total combined background, and there- 
fore are not included in this discussion. 
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TABLE 25. Total observed and expected number of 
£+jets/fi events after all cuts. 





e±jets//i 


/i±jets//i 


Lum (pb 1 ) 


112.6 


108.0 


Observed 


5 


6 


it MC m t (GeV/c x ) 






140 


6.93 ± 1.35 


4.65 ± 1.19 


150 


6.18 ±1.06 


3.31 ±0.83 


160 


4.51 ±0.73 


2.60 ±0.63 


1 70 


Q 70 _1_ A C7 
J.IO I U.JI 




180 


3.11 ±0.46 


1.84 ±0.43 


190 


2.44 ± 0.36 


1.40 ±0.32 


200 


1.83 ±0.27 


1.08 ±0.25 


IF+jets 


0.74 ±0.30 


0.73 ±0.14 


QCD multijet 


0.32 ± 0.26 


0.50 ±0.17 


Z — > nn 




0.17 ±0.08 


Total background 


1.05 ±0.40 


1.40 ± 0.23 



The inclusive jet multiplicity spectrum of the £±jets//i 
data obtained prior to enforcing the A and Ht require- 
ments is compared with that for the expected background 
in Fig. 26. Good agreement is seen in the background- 
dominated 1 and 2 jet bins, but for 3 or more jets, the 
excess due to it production is evident in both ^-tagged 
channels. 

Figures 27 and 28 show the distributions of A vs Ht 
for e+jets//i and /z+jets//i events for data, herwig it 
Monte Carlo (m t = 170 GcV/c 2 ), QCD multijet data, 
and VECBOS !F±jets Monte Carlo events. From these 
figures it is clear that the cuts on A and Ht provide 
a significant improvement in the discrimination between 
signal and background for these channels. 

As described in Sec. VII, the it acceptances are com- 
puted via Eq. 7.5 using Monte Carlo events generated 
with herwig and passed through the D0 detector simu- 
lation. The trigger efficiency for the e±jets/^( channel is 
obtained from the Trigger Simulator (see Sec. V) and has 
been compared with that found for I / F±jets data to esti- 
mate its systematic error, resulting in a value of 99l" 5 %. 
For the /i+jets/fi channel, the trigger efficiency is com- 
puted in the same fashion as for the /z±jets/topo chan- 
nel using data-derived trigger turn-on curves applied to 



it Monte Carlo events and is determined to be 96^5 %. 
The acceptance values after all cuts for seven different top 
quark masses (and for all channels) are given in Sec. X 
(Table 32). Following Eq. 7.4, the expected number of 
it events in the ^±jets//i channels are given in Table 25 
for these same seven masses. Also shown in Table 25 are 
the final numbers of events observed in the data, 5 in 
the e±jets//^ channel and 6 in the /x±jets//x channel. Fi- 
nally, the cross section obtained from the e±jets//x and 
/i.+jets/ /x channels are 6.0 ± 3.6 pb and 11.3 ± 6.6 pb, 
respectively. 
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FIG. 26. Inclusive jet multiplicity spectra for £+jets/n data 
(circles) and expected background (triangles) obtained prior 
to applying the A and Ht requirements. Note that good 
agreement is seen for the > 1 and > 2 jet bins, but the > 3 
jet bin shows a clear excess in the data. 
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FIG. 27. Scatter plots of A vs Ht for the e+jets//i channel 
for (a) VECBOS W+jets MC background, (b) QCD multijet 
background, (c) herwig ti MC events (m t = 170 GeV/c 2 ), 
and (d) data. 
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FIG. 28. Scatter plots of A vs Ht for the /z+jets//n chan- 
nel for (a) VECBOS W+jets MC background, (b) QCD multijet 
background, (c) herwig tt MC (m t = 170 GeV/c 2 ), and (d) 
data. The effective luminosity given for plot (b) is determined 
as the product of the luminosity of the selected multi-jet sam- 
ple and the inverse of the muon misidentification rate. 



VIII. ANALYSIS OF ALL-JETS EVENTS 

As noted in Sec. I, the all-jets channel is discussed in 
detail in Ref. [56] and is only summarized here. 

The signature for the all-jets channel is characterized 
by the presence of six or more high transverse momentum 
jets. Given the overwhelming nature of the background 
to this channel, primarily from QCD multijet produc- 
tion, the challenge of this analysis is to develop selec- 
tion criteria that provide maximum discrimination be- 
tween signal and background, together with an estimate 
of the residual background in the signal region. Several 
kinematic and topological properties of the events were 
investigated, and neural networks employed to properly 
combine all possible sources of discrimination between 
signal and background. In order to improve the signal 
to background ratio, the analysis requires the presence 
of at least one muon-tagged jet in every event. Because 
the data provide an almost pure sample of background 
events, the background model is determined entirely from 
data. The modeling uses untagged events that are made 
to represent tagged events by adding muon tags to one 
of the jets in the event. The cross section is determined 
using fits to the neural network output, and checked us- 
ing a conventional counting method. The cross section 
obtained for m t = 172.1 GeV/c 2 is: 



a ti = 7.3 ± 2.8(stat) ± 1.5(sys) pb. 



(8.1) 



This cross section differs slightly from the value reported 
in Ref. [56] due to an update of the luminosity normal- 
ization. The significance of the excess of tt signal over 
background is estimated by defining the probability P 
of having the expected background fluctuate up to the 
observed number of events. This corresponds to a 3.2 
standard deviation effect, sufficient to establish the exis- 
tence of a ti signal in multijet final states [56,97]. 



IX. SYSTEMATIC UNCERTAINTIES 

The individual uncertainties which affect the accep- 
tance and background are discussed below. A discussion 
of the treatment of the correlations between the uncer- 
tainties can be found in Appendix E. 



A. Sources 

1. Luminosity 

As noted in Sec. Ill, the luminosity is determined with 
the level hodoscopes and is normalized to a world 
average total pp inelastic cross section from CDF [61], 
E710 [62], and E811 [63] Collaborations. The systematic 
uncertainty on the luminosity stems from both the level 
measurement and the world average total pp inelastic 
cross section and is found to be 4.3%. 
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2. Energy scale 

Uncertainty in the jet energy scale affects the cross 
section determination only via the uncertainty in the rel- 
ative scale between data and MC. This uncertainty is de- 
termined by comparing Z(— > ee)+jet events in data and 
MC [98]. Events are selected by requiring two electrons 
with E T > 15 GeV, 82 GeV/c 2 < m ee < 102 GeV/c 2 , 
and at least one jet with Et > 15 GeV. The azimuthal 
bisector of the two electrons is determined and the trans- 
verse momentum of the Z boson is projected along this 
bisector using the electron momentum vectors. The jet 
transverse momenta are also projected along this bisector 
with the contribution from each jet in the event summed 
to form the jet projection. The jet energy projection 
versus the Z — > ee projection is plotted for MC (her- 
wig and vecbos) and data from Run lb, and a lin- 
ear regression fit performed to determine the slope and 
offset of each sample. Comparison of the ratios of the 
slopes (MC/Data) and the differences in the offsets (MC 
- Data) indicate an uncertainty in the jet energy scale 
slope of 4% and an uncertainty in the jet energy scale 
offset of 1 GeV. 



3. Electron identification 

The procedure for determining the electron identifica- 
tion efficiencies is discussed in Sec. IV A 9. The primary 
source of uncertainty in this technique stems from the 
method used to subtract the background under the Z 
boson mass peak. Comparison of several different back- 
ground subtraction schemes [68] is used to determine the 
systematic uncertainties given in Table 6. 

4- High pr and tag muon identification 

As described in Sec. IV B 6, the muon identification 
efficiencies arc determined from a modified version of 
D0GEANT which has additional corrections to account 
for time dependent detector inefficiencies and incorrect 
modeling of the muon track finding efficiency. The time 
dependent correction is applied only to Run la and Run 
lb(preclean) with an uncertainty of 5%, arising primar- 
ily from statistical considerations. The track finding ef- 
ficiency correction varied with detector region with an 
uncertainty of 1.5% in the CF and 2.2% in the EF, also 
arising primarily from statistical considerations. The 
uncertainty arising from the detector simulation is de- 
termined by comparing Z — > fi/j, MC events which are 
passed through the modified version of D0GEANT with 
Z — > fifi data, the difference being a measure of the un- 
certainty This uncertainty varies with run period, detec- 
tor region, and muon identification choice, and includes 
uncertainties from the muon trigger efficiency. The un- 
certainties noted above are added in quadrature to deter- 



mine the systematic uncertainty on the efficiencies given 
in Tables 7-9. 



5. e+jets trigger 

This uncertainty accounts for systematic variations in 
the trigger efficiency for those signal and background MC 
samples that rely primarily on electron triggers (see Ta- 
ble 1). The determinations of the trigger efficiencies for 
each channel are discussed in the subsections of Sees. VI 
and VII. For electron trigger efficiencies determined 
via the Trigger Simulator (e/x: signal and all MC back- 
grounds; ev: signal and all MC backgrounds; e+jets//x: 
signal), the systematic uncertainty is determined by com- 
paring the trigger efficiency of e+jet data events (ob- 
tained from an unbiased trigger) with that found passing 
W(-^ ef)+jet MC events through the Trigger Simulator. 
For electron trigger efficiencies determined directly from 
data: for the ee channel, comparison of the Z(— ► ee)+jets 
trigger rate obtained from unbiased data with that ob- 
tained from passing Z(— > ee)+jet MC through the Trig- 
ger Simulator found a difference of 1% which was taken 
as a measure of the uncertainty; for the e+jets channel, 
studies of the efficiency variation using different samples 
and cuts led to the assignment of an uncertainty of 3%. 

6. l^r+jets trigger 

This uncertainty accounts for systematic variations in 
the efficiency of the $ T triggers (see Table 4). Trigger 
efficiencies from the $ T triggers were obtained from mea- 
sured turn-on curves convoluted with kinematics from 
MC events. The systematic uncertainty is determined 
from the differences in efficiency due to variations in top 
quark mass (for signal) and variations in the A and Ht 
of the events (background). Note that efficiencies for the 
muon triggers were determined from a parameterization 
of the turn-on curves of the muon+jet triggers and the 
systematics have been folded into the uncertainty on the 
muon identification efficiency. 

7. Multiple interactions 

As discussed in Sec. Ill, there were, on average, 1.3 
pp interactions per bunch crossing during Run 1, giving 
rise to additional minimum bias events produced along 
with the high-pr interactions of interest to the present 
analyses. These additional minimum bias events were 
not included in the MC models although they can con- 
tribute to mismeasurement of the primary interaction 
vertex and thus to mismeasurement of lepton and jet 
transverse energies/momenta. For £+jet events, such ef- 
fects were found to be negligible since the presence of 
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three or more hard jets from a single interaction ver- 
tex minimized any potential confusion in determining 
the correct vertex. For the dilepton channels the ef- 
fect is more pronounced, and a systematic uncertainty 
is estimated for all signal and MC-based backgrounds. 
To make this estimate, additional signal and background 
MC samples were produced with one and two minimum 
bias events added. The efficiencies and background pre- 
dictions from these samples are then weighted accord- 
ing to the luminosity distribution of the Run 1 data set 
and compared to the samples for which no minimum bias 
events had been added. The deviations, which vary sig- 
nificantly from channel to channel and between signal 
and background, are taken as an estimate of the uncer- 
tainty. 



8. tt Monte Carlo generator (kinematics) 

The uncertainty on the modeling of kinematic quanti- 
ties (high-pT leptons, jets, and Et) due to imperfections 
in the MC generator is based on efficiency differences 
between the herwig and isajet generators. This uncer- 
tainty is calculated separately for each channel. The pro- 
cedure, which is the same for each channel, is to generate 
a smooth curve summarizing the observed generator dif- 
ference (isajet-herwig/herwig) for top quark masses 
from 140 GeV/c 2 to 200 GcV/c 2 , ignoring any 6-tag or b- 
tag-veto cuts. As seen in Table 26, the dilepton channels 
are parameterized using a constant relative uncertainty 
and the lepton+jets channels are parameterized using an 
exponential function of the top quark mass. The aspect 
of the generator to which the kinematic acceptance is 
most sensitive is the parton showering, herwig has been 
shown to reproduce jet properties well at both the Teva- 
tron [99] and LEP [100]. Reference [99] describes a study 
of the topological properties (spectra of angles and en- 
ergy distribution among jets) in inclusive three and four 
jet events and the authors find that "[a]part from the 
cos(t9*) distributions, the herwig event generator pro- 
vides a reasonably good description of the data while the 
differences between the data and the predictions of [the] 
isajet and pythia event generators are large in many 
distributions." 

TABLE 26. Smoothed kinematic generator uncertainties 
for the eight leptonic channels. 



Channel 


Relative Uncertainty 


Fit 




Applied 


ee 


5.5% 




5% 


e/i 


-4.9% 


5% 


nn 


3.3% 




5% 


ev 


-11.1% 


12% 


e + jets 




exp(4.59 - 


0.0407m t ) 


fj, + jets 




exp(0.546 - 


- 0.0120m t ) 


e + jets/n 




cxp(-0.279 


- O.OlSOmt) 


/j, + jets//i 




cxp(-0.293 


- 0.0124m* ) 



9. tt Monte Carlo generator (b-tagging) 

In addition to kinematic quantities (high-pT leptons, 
jets, and $r), generator imperfections can contribute to 
the uncertainty in the probability that a soft muon will be 
produced and subsequently pass the identification and pr 
cuts (see Sec IV B). Potential sources of uncertainty in- 
clude the branching fraction of b — > ji + X, the branching 
fraction of c — > /j + X for cascade decays, b quark frag- 
mentation, B hadron decay form factors, and uncertain- 
tics associated with misidentified tags. Only the effect of 
the branching fraction of b — > /i + X has been considered. 
In herwig, all b hadrons decay via a spectator model 
with a branching fraction to muons B(b — > fi) = 0.11. 
The particle data book [85] lists the following inclusive 
measurements of B hadron scmilcptonic branching frac- 
tion: 



T(4S I ) inclusive B — ► fj, 
T(4S) inclusive B -> I 
High energy inclusive B 
High energy inclusive B 



10.3 ±0.5% 
10.43 ±0.24% 
H 10.7 ±0.7% 
t 11.13 ±0.29% 



The errors on the inclusive B — > I branching fraction are 
quite small, although the T(4S I ) and high energy mea- 
surements are inconsistent at two standard deviations. 
The uncertainty due to this variation has been increased 
to account for the remaining sources of uncertainty, re- 
sulting in the assignment of a fractional uncertainty of 
10%. 



10. VECBOS 

As discussed in Sec. VII A, the €±jets/topo channels 
use VECBOS to determine the A(W + jets), H T , and Ej, 
cut survival probability for W+jets backgrounds. The 
systematic uncertainty for this procedure is estimated by 
comparing the ,4(W+jets), Ht, and Ej, distributions of 

> 2 and > 3 jet events in data and VECBOS (after adding 
contributions from tt and QCD multijet production to 
the VECBOS sample in the appropriate proportions). For 

> 2 jet events, a 6% difference is seen and for > 3 jet 
events, a 10% difference is seen. Extrapolated to > 4 jet 
events, a 15% uncertainty is estimated. 



11. Background cross section 

As described in Sees. VI and VII B, backgrounds deter- 
mined from MC have their initial cross sections normal- 
ized to either measured or theoretical values and the un- 
certainties are therefore taken from the cited references. 



12. Other simulation 

This uncertainty accounts for additional, channel spe- 
cific, systematic effects due to the simulation and is only 
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TABLE 27. Expected Run 1 dilepton backgrounds and the corresponding statistical and systematic uncertainties (number 
of events) . 









ee 










e/i 
















Zee 


Ztt 


WW 


DYee 


multijet 


Ztt 


WW 


DYrr 


multijet 


Zfifi 


Ztt 


WW 




multijet 


# of evts 


0.058 


0.081 


0.086 


0.056 


0.197 


0.103 


0.077 


0.006 


0.077 


0.579 


0.030 


0.007 


0.068 


0.068 


Statistical 


0.009 


0.008 


0.008 


0.011 


0.044 


0.051 


0.006 


0.004 


0.121 


0.141 


0.015 


0.003 


0.030 


0.010 


Luminosity 




n nn/i 
U.UU4 


U.UU4 


u.uuz 




n nri/i 
U.UU4 


U.UUo 


n nnn 
U.UUU 




n no^ 
U.Uzo 


U.UU1 


n nnn 
U.UUU 


n nnQ 
U.UUo 




Energy Scale 




0.020 


0.022 


0.014 




0.026 


0.010 


0.000 




0.133 


0.007 


0.002 


0.016 




e id 




0.004 


0.004 


0.003 




0.005 


0.002 


0.000 














High pt /i id 












0.012 


0.007 


0.001 




0.040 


0.002 


0.001 


0.005 




e+jets trig 




0.000 


0.000 


0.000 




0.005 


0.004 


0.000 














Mult. Int. 




0.008 


0.009 


0.006 




0.017 


0.012 


0.001 




0.018 


0.001 


0.000 


0.002 




Bkg crsec 




0.010 


0.009 


0.028 




0.010 


0.008 


0.001 




0.059 


0.005 


0.001 


0.010 




Other Sim 




0.050 








0.064 










0.019 








Mis-id e 










0.015 








0.003 












Mis-meas Ifr 


0.009 




























Z fitter 




















0.060 


0.001 


0.000 


0.002 




Total 


0.013 


0.056 


0.027 


0.034 


0.046 


0.089 


0.021 


0.004 


0.121 


0.218 


0.026 


0.004 


0.036 


0.010 



TABLE 28. Expected Run 1 i?+jets backgrounds and the corresponding statistical and systematic uncertainties (number of 
events) . 





e+jets/topo 


/i+jets/topo 


e + jets/ fj, 




H+jets/n 






W+jets 


multijet 


IV+jets 


multijet 


W^+jets 


multijet 


W+jets 


Zfifi 


multijet 


# of evts 


4.135 


0.379 


3.324 


0.993 


0.738 


0.316 


0.726 


0.170 


0.500 


Statistical 


0.464 


0.139 


0.437 


0.347 


0.044 


0.246 


0.118 


0.036 


0.052 


Luminosity 
















0.007 




Energy Scale 


0.207 




0.179 










0.017 




High-pT A* id 
















0.022 




Tag n id 
















0.010 




^T+jets trig 






0.166 










0.002 




Mult. Int. 
















0.000 




VECBOS 


0.616 




0.665 














Bkg crsec 
















0.051 




Berends scaling 


0.413 




0.369 




0.292 










Mis-id e 












0.066 








Tag rate 










0.074 




0.073 






Mis-id n 








0.298 










0.100 


fj, multijet 








0.199 










0.100 


Tag probability 


















0.075 


Z fitter 
















0.042 




Total 


0.899 


0.139 


0.911 


0.498 


0.304 


0.255 


0.139 


0.081 


0.168 
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TABLE 29. Expected Run 1 ev and all-jets expected back- 
grounds (number of events) and the corresponding statistical 
and systematic uncertainties. Uncertainties labelled Tag rate 
norm, Tag rate fn, and ti corr are for the all-jets channel only 
and correspond respectively to uncertainties associated with 
the normalization of the muon tag rate, the functional form 
of the muon tag rate, and corrections to the background for 
ti signal. The systematic uncertainties on the all-jets channel 
are discussed in detail in Ref. [56] . 









ev 




all-jets 




WW 


wz 


W+jets 


multijet 


multijet 


# of evts 


0.161 


0.017 


0.543 


0.471 


24.8 


Statistical 


0.028 


0.002 


0.272 


0.103 


0.7 


Luminosity 


0.007 


0.001 


0.023 






Energy scale 


0.040 


0.004 


0.136 




1.0 


e id 


0.004 


0.000 


0.013 






e+jets trig 


0.005 


0.001 


0.016 


0.014 




Mult. Int. 


0.009 


0.001 


0.030 






VECBOS 






0.086 






Bkg crsec 


0.016 


0.002 








Other Sim 








0.104 




Mis-id e 








0.034 




Tag rate norm 










1.2 


Tag rate fn 










1.2 


tt corr 










1.0 


Total 


0.053 


0.005 


0.319 


0.151 


2.4 



included for the Z — > tt background to the ee, e/i, and 
H/j, channels and for the QCD multijet background to the 
ev channel. As described in Sees. VI A- VIC, the jet cut 
survival probabilities for the Z — > tt — > it backgrounds 
are obtained from Z(— ► ee)+jet data. The primary lim- 
itation of this technique is the limited statistics of the 
Z(— ► ee)+jet data set, which is taken as the dominant 
uncertainty. As described in Sec. VI D, the QCD multijet 
background is obtained as the mean of two independent 
procedures. The difference between the two procedures 
is taken as a systematic uncertainty. 

13. Berends scaling 

As noted in Sec. VII A, the assumption of A?j ets or 
Berends scaling (see Eq. 7.6) is used by the £+jets/topo 
channels to compute the background from VF+jets. In or- 
der to investigate the validity of this assumption, a num- 
ber of data sets were examined: VF+jcts, QCD multijet, 
Z+jets, photon+jets, and VECBOS WK+jets production. 
For each sample the number of events with a minimum 
jet multiplicity of n — 1 and n— 2 were used to predict the 
number of events with a minimum jet multiplicity > n. 
These predictions were compared with observations and 
the maximum differences are given in Table 30. Based 
on these values an uncertainty of 10% is assigned for the 
uncertainty due to Berends scaling. 

As described in Sec. VII B, the calculation of the 
W+jets background for the e+jets//z channel is deter- 



TABLE 30. Maximum deviation between predictions from 
Berends scaling and observation for several data sets. 



Data Set 


Maximum Deviation (%) 


W+jets 


3.1 


QCD multijet 


< 10 


Z+jets 


< 4 


Photon + jets 


< 5 


VECBOS W+jets 


< 1 



mined via Eq. 7.12 where is obtained by applying 
the ti tag rate to the measured excess for e+ 3 or more 
jets as determined from Berends scaling (Eq. 7.23). In 
addition to the uncertainty from Berends scaling of 10%, 
there is a significant uncertainty in the ti tag rate de- 
termined from MC, leading to a total uncertainty of 40% 
which has been included under the Berends scaling head- 
ing for the e+jets/ \x channel. Note that Berends scaling 
is not used for the /z+jets//z channel. 

14- Electron misidentification rate (mis-id e) 

As described in Sees. IV A, VI, and VII B, determina- 
tion of the background from multijet events in which a jet 
is misidentified as an electron is based on an independent 
measurement of the electron "misidentification rate." For 
the ee, e/i, and ev channels, these misidentification rates 
were determined by counting the number of loose electron 
candidates found in a sample of QCD multijet events 
containing one electromagnetic cluster that passed the 
extra-loose electron identification requirements. The un- 
certainties on this procedure are dominated by the statis- 
tics of the extra-loose electron sample. For the e+jets/^ 
sample, the misidentification rate described in Sec. VII B 
depends on the jet multiplicity from which an uncertainty 
of 21% was estimated. Note that for the e+jets/topo 
channel, the background from QCD multijet events is 
handled differently and did not make use of an electron 
"misidentification rate." 



15. Mismeasured B~t 

As noted in Sec. VIA, for the ee channel the back- 
ground from Z{— > ee)+jets is determined directly from 
data, but since Z(— > ee)+jet events have no real Ifr, a 
l^T mis-measurement rate, computed from QCD multi- 
jet data as a function of jet multiplicity, is applied. The 
uncertainty on this procedure is obtained by varying the 
triggers and selection criteria used to collect the initial 
multijet sample, and is assigned a value of 15%. 
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TABLE 31. Efficiency times branching fraction (e x B) and statistical and systematic uncertainties (all in percent) for 
rn t = 170 GeV/c 2 . 









fifi 








P -4- lf>t ^ / 1 1 


/ 1 -4- ir>t ^ / 1 1 


sll-jcts 


£ X B 


0.165 


0.349 


0.106 


0.263 


1.288 


0.911 


0.568 


0.371 


1.963 


Statistical 


0.002 


0.004 


0.002 


0.008 


0.020 


0.046 


0.017 


0.037 


0.151 


Energy scale 


0.011 


0.020 


0.008 


0.066 


0.169 


0.137 


0.026 


0.008 


0.112 


Electron id 


0.008 


0.009 




0.006 


0.044 




0.022 






High-^TM id 




0.033 


0.007 






0.098 




0.048 




Tag id 










0.005 




0.022 


0.022 


0.137 


e+jets trigger 


0.001 


0.018 




0.008 


0.058 




0.028 






^T+jets trigger 












0.046 




0.019 


0.098 


Mult. Int. 


0.016 


0.057 


0.005 


0.014 


0.000 


0.000 


0.000 


0.000 




Generator (kin) 


0.008 


0.017 


0.005 


0.032 


0.126 


0.203 


0.034 


0.034 




Generator (6 tag) 










0.021 


0.017 


0.057 


0.037 




Z fitter 






0.003 










0.019 




Total error 


0.023 


0.074 


0.013 


0.076 


0.225 


0.272 


0.084 


0.086 


0.253 



16. Tag rate 

The VF+jets background to the £+)ets//j, channels is 
obtained, as a function of jet Et and 77, by multiply- 
ing the number of (QCD multijet and tt subtracted) un- 
tagged £+jets events by a tag rate determined from mul- 
tijet data. As described in Sec. VII B, the accuracy of the 
tag rate was studied by applying it to a number of differ- 
ent data sets and comparing the predicted and observed 
values (see Fig. 24). Variation not due to statistics is 
calculated to be 8.2% [69] and rounded upward to 10%. 

17. Muon misidentification rate (mis-id fj,) 

The /z+jets/topo and p+jets/p channels both employ 
the use of an "isolated muon misidentification rate" to de- 
termine the background from QCD multijet events. As 
described in Sec. VII A, this misidentification rate is de- 
pendent on the jet multiplicity and is computed from 
samples of QCD multijet events with $ T < 20 GeV as 
the ratio of the number of isolated-p + n jet events to 
the number of non-isolated- p + (n + 1) jet events. The 
primary source of uncertainty in this measurement is the 
statistical precision of the control samples, leading to an 
uncertainty of 30% for the four-jet samples used for the 
/z+jets/topo channel and 20% for the three-jet samples 
used for the p+jets/p channel. 

18. [i multijet 

Both the /x+jets/topo and p+jets/p channels have 
background from QCD multijet events which contain a 
muon from b or c quark decay that is misidentified as 
an isolated muon. Both channels rely on multijet con- 
trol samples to model this background. Differences in 
key kinematic distributions between the multijet control 
samples and the true background are accounted for in the 
uncertainty discussed here. As discussed in Sec. VII A, 



the QCD multijet background to the /x+jets/topo chan- 
nel is obtained by applying a survival probability to pass 
the E^, A, and Ht cuts (determined from n+1 jet data) 
to an n jet control sample. Comparisons of the A and Ht 
distributions for the n and n+ 1 jet sample lead to an es- 
timated uncertainty of 20%. Similarly, for the p+jets/p 
channel, the QCD multijet background is determined by 
applying a tag probability to the jets in a multijet con- 
trol sample of non-isolated p + 3 jet events on which all 
kinematic cuts (including A and Ht) have been applied. 
Differences in the A and Ht distributions between the 
multijet control sample and the true background sample 
lead to the assignment of an uncertainty of 20%. 

19. p tag probability 

As described in Sec. VII B, for the p+jets/p channel 
the QCD multijet background is determined by applying 
a tag probability, derived from ti MC events, to a multi- 
jet control sample. An uncertainty of 15% is assigned to 
this tag probability to account for the fact that the prob- 
ability is averaged over the CF and EF detector regions 
and that the MC sample has not been subjected to the 
corrections described in Sec. IV B 6. 



20. Z boson mass fitter (Z fitter) 

As described in Sees. VIC and VII B, the mi and 
/x+jets/jti channels reduce their background from Z — > fifi 
events by cutting on a minimized \ 2 fit for the muon 
pair mass to give Mz and for $^ al to equal the pt of the 
Z boson, in effect "fitting for the Z." Consideration of 
the muon momentum resolution and variation of the JfiT 
resolution parameterizations used for both data and MC, 
lead to the estimate of a systematic uncertainty of 10% 
for this procedure. 
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TABLE 32. Efficiency x branching fraction (in percent) for all eight leptonic channels for m t = 140-200 GeV/c 2 . Uncer- 
tainties correspond to statistical and systematic contributions added in quadrature. 



m t (GoV/c 2 ) 


140 


150 


160 


170 


180 


190 


200 


ee 
efj, 

MM 

ev 

e+jcts/topo 
/^+jcts/topo 
e+jots//j 
^i+jcts/^ 


0.106 ± 0.015 
0.214 ± 0.045 
0.055 ± 0.008 
0.156 ± 0.046 
0.597 ± 0.256 
0.451 ± 0.194 
0.364 ± 0.069 
0.255 ± 0.064 


0.129 ± 0.018 
0.252 ± 0.053 
0.069 ± 0.010 
0.201 ± 0.058 
0.801 ± 0.264 
0.621 ± 0.235 
0.469 ± 0.078 
0.262 ± 0.065 


0.152 ± 0.021 
0.302 ± 0.064 
0.088 ± 0.012 
0.224 ± 0.065 
1.036 ± 0.236 
0.810 ± 0.271 
0.491 ± 0.077 
0.295 ± 0.071 


0.165 ± 0.023 
0.349 ± 0.074 
0.106 ± 0.013 
0.263 ± 0.076 
1.288 ± 0.225 
0.911 ± 0.272 
0.568 ± 0.084 
0.371 ± 0.086 


0.187 ± 0.026 
0.389 ± 0.082 
0.119 ± 0.016 
0.315 ± 0.091 
1.479 ± 0.197 
1.058 ± 0.276 
0.656 ± 0.095 
0.405 ± 0.093 


0.199 ± 0.028 
0.429 ± 0.092 
0.133 ± 0.018 
0.335 ± 0.096 
1.558 ± 0.174 
1.164 ± 0.276 
0.709 ± 0.099 
0.423 ± 0.096 


0.209 ± 0.029 
0.440 ± 0.093 
0.137 ± 0.018 
0.360 ± 0.104 
1.703 ± 0.158 
1.291 ± 0.278 
0.721 ± 0.102 
0.443 ± 0.100 



X. CROSS SECTION RESULTS 

The preceding sections describe nine analyses that ex- 
tract data samples rich in ti events. For an individual 
channel i, the cross section is determined from the rela- 
tion 



A(m t )i ■ Li 



(10.1) 



where A(nit) is the acceptance (efficiency times branch- 
ing fraction) for a top quark mass of m t , Li is the inte- 
grated luminosity, N t is the number of observed events, 
and Bj is the number of expected background events 
from source j. The efficiency times branching fraction 
values for all eight leptonic channels for m t = 140-200 
GeV/c 2 are given in Table 32. The numbers of observed 
events, along with those expected from signal and back- 
ground, the integrated luminosity, and the final measured 
cross sections (for m t = 172.1 GeV/c 2 ) for each chan- 
nel are summarized in Table 33. The value of m t — 
172.1 GeV/c 2 is D0's combined dilcpton and lepton+jets 
mass measurement [92,101]. The cross section results 
for the various channels (and several combinations) are 
compared in Fig. 29, and are seen to be in good agree- 
ment with one another and with theoretical expecta- 
tions [45-47]. Complete details of the 39 observed lep- 
tonic events are given in Ref. [102]. 

The combined ti production cross section is determined 
from the analog of Eq. 10.1: 



o-(m t ) t 



J~] A(m t )j ■ A 



(10.2) 



All Channels Combined 
Leptonic Combined 

ee 
ell 

mi 

ev 

Dilepton Combined 

e+jets/topo 
il+jets/topo 
e+jets/il 
(l+jets/(J, 

L+jets combined 

All jets 

Theory 



10 



~~ 1 
15 



% (pb) 



5.7 ± 1.6 pb 

5.3 ± 1.7 pb 

2.4 + 4.6 pb 

6.8 ± 4.6 pb 
2.1 ±8.8 pb 
9.1 ±7.2 pb 

6.0±3.2pb 

2.8 ±2.1 pb 
5.6±3.7pb 

6.0 ± 3.6 pb 
11.3 ± 6.6 pb 

5.1 ± 1.9 pb 
7.3±3.2pb 
4.7 - 5.9 pb 
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FIG. 29. D0 measured ti production cross section values 
for all channels, assuming a top quark mass of 172.1 GeV/c 2 . 
The vertical line corresponds to the cross section for all chan- 
nels combined and the shaded band shows the range of theo- 
retical predictions [45-47,50]. 



where the sum i is over all nine channels and the sum j is 
over all background sources in all nine channels. Recall 
(see Sec. I) that all channels are, by construction, orthog- 
onal. As discussed in Appendix E, the determination of 
the cross section takes into account the correlated un- 
certainties between the inputs to Eq. 10.2. Plotting the 
cross section values and their uncertainties for a range of 
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TABLE 33. Summary of observed number of events, expected signal and background, integrated luminosity, and cross 
section for all nine channels at ra t = 172.1 GeV/c 2 . Uncertainties correspond to statistical and systematic contributions added 
in quadrature. 







Total 


Signal 


Total 








AobB 


sig+bkg 


{rat = 172.1) 


background 


J Cdt (pb" 1 ) 


a(pb) 


ee 


1 


1.68 ±0.23 


1.20 ±0.17 


0.48 ±0.10 


130.2 ±5.6 


2.37 ±4.58 


e/i 


3 


2.45 ±0.53 


2.19 ± 0.47 


0.26 ±0.16 


112.6 ±4.8 


6.81 ±4.59 




1 


1.39 ± 0.30 


0.64 ± 0.09 


0.75 ± 0.24 


108.5 ±4.7 


2.11 ± 8.79 


ev 


4 


2.87 ± 0.71 


1.68 ± 0.49 


1.19 ± 0.38 


112.3 ± 4.8 


9.12 ± 7.23 


Tji 1 £i7~i^ i~\y\ r'Avn nincin 
-L/llCjJ LUI1 UUlIlUIlicU. 


Q 


O.Oa ZC _L .*iO 


\ 71 4- 1 f!7 


9 RQ 4- fl (SIS 




(S 09 4 1 91 


e+jets/topo 


9 


13.16 ± 1.67 


8.64 ± 1.47 


4.51 ±0.91 


119.5 ±5.1 


2.83 ±2.05 


/it+jets/topo 


10 


9.84 ± 1.62 


5.52 ± 1.62 


4.32 ± 1.04 


107.7 ±4.6 


5.60 ±3.71 


e+jets//i 


5 


4.65 ±0.54 


3.59 ±0.55 


1.05±0.40 


112.6 ±4.8 


5.98 ±3.56 


/i±jets//i 


6 


3.62 ±0.52 


2.22 ±0.52 


1.40 ±0.23 


108.0 ±4.6 


11.27 ±6.60 


£±jets combined 


30 


31.27 ±3.52 


19.98 ±3.52 


11.28 ± 1.97 




5.10 ± 1.85 


Leptonic combined 


39 


39.66 ± 4.65 


25.69 ± 4.41 


13.97 ± 2.22 




5.31 ± 1.72 


all-jets 


41 


37.40 ± 2.92 


12.60 ±2.12 


24.80 ± 2.37 


117.9 ±5.1 


7.33 4 3.20 


All channels total 


80 


77.06 ± 6.19 


38.29 ± 5.34 


38.77 ± 3.32 




5.69 ± 1.60 



top quark masses gives the band shown in Fig. 30. Also 
shown are the theoretical expectations for the ti cross 
section as a function of m t [45-47,50]. Combining this 
cross section result with the combined D0 dilepton and 
lepton±jets mass measurement [92] gives the point with 
error bars shown in Fig. 30. 

In addition to the final cross section and mass result, 
it is also instructive to compare the properties of the 
ti candidate events with expectations. These is exam- 
ined in Figs. 31-34 which show the distributions of 
the ti candidates (shaded histograms), tt Monte Carlo 
(unshaded histogram), expected background (open tri- 
angles), and expected signal plus background (solid cir- 
cles) for various quantities. Overall, these plots show 
better agreement between the candidate and ti ± back- 
ground distributions than between the candidate and the 
background-only distributions. 



XI. CONCLUSIONS 

Nine analyses have been described which select event 
samples dominated by tt production. A total of 39 events 
are found in the leptonic channels with an expected back- 
ground of 14.0 ± 2.2. Combining these results with the 
integrated luminosity and signal efficiency (at m t = 172.1 
GeV/c 2 ), the ti production cross section for the leptonic 
channels is determined to be 



5.31 ± 1.34(stat) ± 1.08(sys) pb. 



(11.1) 



This cross section differs slightly from the value reported 
in Ref. [103] due primarily to an updated luminosity nor- 
malization, and to a lesser extent to minor changes in the 
background estimation for some channels and to the use 
of a slightly different top mass. 



15 



3 
a 

o 
in 

8 
u 



10- 



5 - 








150 160 170 180 190 200 
Top Quark Mass (GeV/c 2 ) 

FIG. 30. D0 measured ti production cross section as a 
function of the top quark mass (mt) for the leptonic channels 
(shaded band) and at the D0 measured top quark mass (point 
with error bars). Also shown are the upper (fi = rrit/2) and 
lower (fi — 2m t ) bounds for four theoretical predictions of the 
ti cross section as a function of m t : Laenen et al. [LL] (dotted 
lines) [45], Berger et al. [LL] (solid lines) [46], Bonciani et al. 
[NLL] (dashed lines) [47], and Kidonakis [NNLL] (dot-dash 
lines) [50], where the labels LL, NLL, and NNLL indicate lead- 
ing-log, next-to-leading-log, and next-to-next-to-leading-log 
resummation calculations respectively. 
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FIG. 31. Et(pt) and r\ distributions for electrons (muons) for leptonic tt candidates (shaded histogram), ti MC [herwig, 
nit = 170 GeV/c 2 ] (unshaded histogram), expected background (open triangles), and expected signal plus background (solid 
circles). The measured muon pr for e\x candidate 58796-7338(417) is 280.0 GeV/c and is therefore off scale in plot (b). As 
given in Ref. [102], the event label corresponds to run number and event number (in an early event numbering scheme, this 
event became well know as "event 417" and thus retains this label parenthetically). 
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FIG. 32. Jet Et distributions for dilepton (a)-(b) and i!+jets (c)-(f) ti candidates (shaded histogram), ti MC [herwig, mt = 
170 GeV/c 2 ] (unshaded histogram), expected background (open triangles), and expected signal plus background (solid circles). 
The dilepton candidate histograms [(a)-(b) shaded] have been multiplied by a factor of 0.25 for presentational clarity. 
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FIG. 33. Jet r\ distributions for dilepton (a)-(b) and £+jets (c)-(f) tt candidates (shaded histogram), ti MC [herwig, m t = 
170 GeV/c 2 ] (unshaded histogram), expected background (open triangles), and expected signal plus background (solid circles). 
The dilepton candidate histograms [(a)-(b) shaded] have been multiplied by a factor of 0.25 for presentational clarity. 
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FIG. 34. distributions for the ee (a) and e+jets/topo (c) channels, and Ifr distributions for the e/i (b) and /i+jets/topo 
and /x+jets//i (d) channels: ti candidates (shaded histogram), ti MC [herwig, m t = 170 GeV/c 2 ] (unshaded histogram), 
expected background (open triangles), and expected signal plus background (solid circles). The ee and e/i candidate histograms 
[(a)-(b) shaded] have been multiplied by a factor of 0.25 for presentational clarity. The measured $ T for e/i candidate 
58796-7338(417) is 182.9 GeV and is therefore off scale in plot (b). 
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For the all-jets channel, summarized in Sec. VIII and 
described in detail in Ref. [56], a total of 41 events arc 
found with an expected background of 24.8 ± 2.4 events. 
Combining the lcptonic and all-jets channels gives a total 
of 80 candidates with an expected background of 38.8 ± 
3.3 events. This combination results in a tt production 
cross section of 



5.69 ± 1.21(stat) ± 1.04(sys) pb. 



(11.2) 



As can be seen in Fig. 29, the tt production cross sec- 
tions obtained for the individual channels are in good 
agreement with the combined cross section and with that 
from theory [45-47]. And as shown in Fig. 1(b), the com- 
bined cross section is in excellent agreement with D0's 
previously reported values. The current level of uncer- 
tainty on QCD predictions for the tt production cross 
section [46,47] is seen in Fig. 30 to be about ±0.3 pb, 
less than 20% of the current experimental uncertainty 
Run II of the Fermilab Tevatron is expected to provide 
an experimental uncertainty on the tt cross section of 
around ±9% (~ 0.6 pb) in 2 fb~ , limited by system- 
atic uncertainties [104]. This will begin to place restric- 
tions on the various QCD predictions and provide strin- 
gent tests for non-standard production and decay mech- 
anisms. In the longer term, the systematic limitations 
on the measurement of the tt production cross section at 
the CERN Large Hadron Collider is expected to be less 
than 10% [105]. 
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APPENDIX A: ENERGY SCALE CORRECTIONS 

Gluon radiation and fragmentation can alter a parton's 
original energy and direction before its remnants interact 
and are measured in the calorimeter (E^l as ). Also, ac- 
companying spectator interactions, not associated with 
the hard scattering, can deposit energy within a jet. In 
addition, fluctuations in interactions in the detector can 



provide changes to i^eas- For example, emitted par- 
ticles, especially hadrons, can produce very wide show- 
ers in the calorimeter that can affect the fraction of en- 
ergy (1 — 5*) contained within any fixed size cone. Also, 
most of the absorber is composed of uranium, the ra- 
dioactive decay of which can deposit significant energy 
in the calorimeter. Finally, the signal response (R) of 
the calorimeter to a jet is dominated by any difference 
of response to electrons (or photons) relative to charged 
hadrons [106,107], and by any energy deposited in unin- 
strumented or non-uniform parts of the detector. The en- 
ergy from spectator interactions and uranium noise pro- 
vides a total offset (O) that must be corrected. 

Other than correcting for spectator interactions, only 
detector effects are considered in the energy calibration 
of jets. A jet's particle level energy (-E^'j) is defined as 
the energy of a jet found from final state particles using 
a similar cone algorithm to that used at the calorimeter 
level. The calibration procedure [107] provides E 3 ^ cl from 
i^ as through the relationship: 



rnjet _ (^incas 

p' c1 " R(l-S) ' 



(Al) 



The calibration is performed separately but identically 
in data and in the Monte Carlo, with the O and 5 correc- 
tions applied to jet energies to extract the particle-level 
values E 3 ^ cl . 

The offset O is estimated as follows. The difference 
in E T density in (77, <f>) space between single and double- 
interaction events, which was obtained with a minimum 
bias trigger, is defined to be the contribution of the un- 
derlying event to single interactions. The contribution 
from noise is obtained from this same sample by sub- 
tracting the Et for the underlying event from the Et 
density in single interactions. The total systematic un- 
certainty for the offset in Et density varies from 100 MeV 
to 300 MeV, depending on the value of r/. 

The showering of a jet's fragments in the calorimeter 
causes energy to leak out of, or into, any jet cone. To 
quantify this effect, jets are generated using the herwig 
program [54] , and reconstructed from their original final- 
state particles. These are subsequently replaced with 
electron or hadron showers from test beam data, and 
reconstructed using our cone algorithm, thereby defining 
a jet shower. The total shower energy is normalized to 
that of the original final-state particles. The ratio of the 
contained shower energy to that of the original energy 
(=1 — 5) is calculated as a function of AR. For central 
jets with AR = 0.5, 5 lies between 0.01 and 0.03, de- 
pending on jet energy, with a systematic uncertainty of 
1% on 1 - 5. 

The IpT in direct-photon candidate events (composed 
of true direct photon events and background dijet events 
where a tt° is back to back with a hadronic jet) is used to 
determine the response of the calorimeter to jets. Differ- 
ences in response between the photon and the recoiling 
hadronic system produce an overall imbalance in trans- 
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verse energy in the calorimeter, giving rise to I^t- In 
these events, the absolute response R of the leading jet 
can be determined from other well-measured quantities 
in the event: 



R = 1 + 



T ■ nj, 
Et 1 ' 



(A2) 



where E T (> 15 GeV) is the transverse energy of the pho- 
ton and tit is the unit vector along the photon's trans- 
verse momentum. Since both the Et of the photon and 
the direction of the probe jet are well-measured, the en- 
ergy estimator E' can be defined 



E 1 = E T cosh(rj iet ). 



(A3) 



Measuring the correlation of R with E', and E^ cas with 
E' } determines the dependence of R on E^ cas . 

Backgrounds to direct photons are a source of uncer- 
tainty for this analysis, particularly in collider data. In- 
strumental background from highly electromagnetic jets 
is limited by tight isolation criteria. The residual bias 
to the measured response is 1.4%. The remaining back- 
ground consists mostly of W{— > ef)+jets production, 
and corresponds to about 0.5%. 

In the calibration, because of the rapidly falling pho- 
ton cross section, energies of central jets are limited to 
< 150 GeV. Exploiting the uniformity of the detector, 
events with EC jets are used to measure the response to 
higher energy jets. Sensitivity to the number of multiple 
interactions in an event results in a 2% systematic un- 
certainty. Because uncertainties in the measurement of 
the energy scale of low Et jets are quite large, a Monte 
Carlo direct-photon sample is used for this region, and 
provides a systematic uncertainty of about 3.5%. 

The total correction is shown in Fig. 35. It rises to a 
maximum of 1.18 at Et — 70 GeV, followed by a slow 
fall to 1.12 at E T ^ 500 GeV. The upper and lower 
dashed lines correspond to one standard deviation up- 
per and lower excursions on the total uncertainty, taken 
as the addition in quadrature of the independent effects 
discussed above. 



APPENDIX B: MAIN-RING VETO 

As noted in Sec. Ill, particles lost from the Main 
Ring can affect the measurements of the outer hadronic 
calorimeter and muon system. The primary losses occur 
every 2.4 seconds when the protons are injected into the 
Main Ring and 0.3 seconds later as the beam, which is 
being accelerated, passes through transition [108]. The 
injection from the Booster into the Main Ring causes the 
bunch to widen, and, consequently, a greater amount of 
beam leaks out of the beampipe. After a few full circuits 
of the beam in the Main Ring, the bunch coalesces and 
is mainly confined to the beampipe. Additional losses 
need to be accounted for in the case when the passage 
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FIG. 35. Total correction to the energy scale for central 
(H < 0.5) jets. 



of the proton beam coincides with the pp crossing in 
the Tevatron (which occurs every 3.5 /its). With each 
pass, errant particles from the bunch scatter outside the 
beampipe causing energy deposition in the outer layers 
of the calorimeter and multiple tracks in the muon sys- 
tem. Because the electromagnetic calorimeter and track- 
ing systems are shielded from these losses, the electron 
triggers are not significantly impacted. However, jet and 
especially muon triggers are affected, and it is necessary 
to veto events from jet, iff 1 , and muon triggers that oc- 
cur during periods of Main-Ring activity. During the 
course of the run, several schemes were used to eliminate 
such events without introducing unnecessary deadtimc: 

• MRBS-LOSS (mrbs): The trigger is disabled for 0.4 s 
after a proton bunch is injected into the Main Ring. 
This vetoes events during injection and transition 
and provides a brief recovery time for the muon 
and calorimeter systems. The typical deadtime for 
MRBS-LOSS veto is rts 17%. 

• micro-blank (mb): The trigger is disabled for 
events where Main- Ring bunches are present during 
the livetime of the muon system which is « ±800 
ns centered on the pp crossing time. The calorime- 
ter livetime is somewhat longer (rs 2 /is), so this is 
therefore not completely efficient for vetoing events 
with Main-Ring energy in the calorimeter. The typ- 
ical deadtime for micro-blank is w 7 %. 

• max-live (ml): The trigger is disabled during pe- 
riods of overlap between MRBS and MB. This cor- 
responds to the first few passes of newly injected 
beam through the detector. 
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• GOOD-CAL (gc): The trigger is disabled during pe- 
riods of overlap between MRBS and MB and during 
MB periods of highest intensity beam leakage. This 
leakage is measured by a set of scintillator arrays 
surrounding the Main-Ring beampipe upstream of 
the D0 detector. 

• GOOD-BEAM (gb): The trigger is disabled during 
periods of either mrbs or MB. GOOD-BEAM is the 
cleanest possible running condition. 

The Main-Ring veto used for each trigger is given in 
Tables 1-5. However, by default, all channels required 
GOOD-BEAM for the offline analyses. As will be noted 
in Sees. C, VIA, and VII A, for the ee and e+jets/topo 
channels it is possible to remove this offline requirement 
on GOOD-BEAM and recover a significant fraction of the 
data lost to it. 



APPENDIX C: MAIN-RING RECOVERY 

As noted in Sec. Ill, all triggers used in the present 
analyses, being combinations of electron, muon, jet, and 
$^ al triggers, suffer some loss from the vetoing of events 
that coincide with activity in the Main Ring. Due to the 
location of the Main- Ring beam pipe within the detector, 
the fine hadronic (FH) and electromagnetic sections of 
the calorimeter and the tracking systems are well shielded 
from this background, so electron and photon measure- 
ments are not significantly affected. However, hadronic 
jet (and thus If^) an d muon measurements are affected. 
The effect on the hadronic calorimeter gives rise to fake 
jet backgrounds and mismeasured $^ arising from ei- 
ther large positive signals, if the Main-Ring losses coin- 
cide with the Tevatron beam crossing (micro-blank), or 
from large negative signals for Tevatron beam crossings 
that were preceded by Main-Ring losses (mrbs). In the 
latter case, the output voltage of the calorimeter preamps 
slowly decreases toward zero, causing the difference be- 
tween a peak and the baseline to become negative. As 
discussed in the following paragraphs, these effects on the 
hadronic calorimetry can be minimized with the proper 
corrections. The effect on the muon system is to decrease 
the overall muon-finding efficiency by less than 10% dur- 
ing periods of Main-Ring activity, with most of the inef- 
ficiency localized to the regions near the Main Ring. 

The ee, e+jets/topo, and /Lt+jets/topo channels all re- 
trieve Main-Ring events and correct jets and Tj^ X in the 
same way. 

Normal jets (those from periods when the Main Ring 
is not active [good-beam]) typically have at most 10% 
of their energy in the outer, coarse hadronic (CH) region 
of the calorimeter. During periods of Main- Ring activity 
(micro-blank), a significant enhancement is seen in the 
number of jets with <f> values close to that of the Main 
Ring (<fi w 1.7), and the vast majority of these jets have 
CH energy fractions between 60 and 90%. Therefore, for 
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FIG. 36. Effect of Main-Ring Jff 1 corrections: Jff^ 1 vs <j> f° r 
mrbs events for (a) no correction, (b) CH correction, and (c) 
total correction. 



jets in the vicinity of the Main Ring (1.5 < <fi < 2.0) that 
have CH energy fractions greater than 20%, the CH en- 
ergy is simply removed [68]. This correction causes the 
jet Et to be biased low due to the fact that some "real" 
CH energy is also removed, but as this only affects a 
small fraction (< 2%) of jets in Main-Ring events, it is 
not a significant concern. Since jets in top quark events 
are very energetic, the removal of the CH energy typi- 
cally leaves the jet Et well above threshold. Therefore, 
the loss in efficiency is minimal, affecting only a small 
fraction of the 2% of jets in Main-Ring events that are 
corrected. For events with large negative signals (mrbs) 
there is also only a small reduction in efficiency, so jets 
in these events are not corrected. 

For If^ the situation is more complicated and re- 
quires corrections for both the large positive signals in 
MICRO-BLANK events and the large negative signals in 
MRBS-LOSS events. The vast majority of these events 
are corrected simply by removing the CH energy from 
the $^ al calculation. This can be seen in Fig. 36 which 

shows TjS^ vs <j) for mrbs events. Fig. 36(a) is without 
any correction and shows a large number of events with 
large $^ al pointing towards the Main Ring. As shown 
in Fig. 36(b), where the CH energy has been removed, 
most of the events with large iff pointing towards the 
Main Ring have been corrected. Although this procedure 
does remove some positive energy that would normally 
be included, it docs not degrade the iff 1 resolution ap- 
preciably due to the fact that normal (non-Main-Ring) 
events characteristically have a low (< 10%) CH energy 
fraction. Unfortunately, some events with large iff 1 in 
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the vicinity of the Main Ring persist after the removal 
of the CH energy. These events appear primarily in the 
region of the intercryostat detector (ICD) and massless 
gap (MG). To correct such events, a vector sum is cal- 
culated for all cells in the ICD and MG that have neg- 
ative energy below a given threshold, and this vector is 
then subtracted from the $^ vector. These thresholds 
were determined from comparisons of the negative en- 
ergy spectra of the ICD and MG cells for good-beam 
(non-Main- Ring) and MRBS-LOSS events [68]. In addi- 
tion to removing all unwanted negative energy, as seen 
in Fig. 36(c), this procedure brings the resolution to 
an approximately normal level. 

APPENDIX D: e/x NEURAL NETWORK 
ANALYSIS 

To further explore this channel, a more sophisticated, 
independent analysis is also performed. The basic scheme 
begins with a loose selection and then uses a neural net- 
work (NN) to maximize the significance. 

This newer analysis is based on the same data set and 
trigger requirements described in Sec. VI B, and the ini- 
tial selection is similar. After passing the trigger require- 
ment, events are required to have at least 1 loose elec- 
tron with Et > 15 GeV, |ry| < 2.5 and at least 1 loose 
muon with p T > 15 GeV/c. A cut of AR(fj,,]et) > 0.5 
is then applied to reduce background from QCD multi- 
jet events containing a misidentified electron and a non- 
isolated muon. To remove QCD multijet events in which 
a misidentified electron and an isolated muon arise from 
the same jet, a cut of AR(e, /i) > 0.25 is applied. As can 
be seen in Table 34, at this stage the backgrounds from 
QCD multijet events containing a misidentified electron 
and an isolated muon from the semi-leptonic decay of a b 
or c quark and W(— > /z^)+jets events in which one of the 
jets is misidentified as an electron is still non-negligible. 
A cut $p al > 15 GeV eliminates the multijet events with 
low .^^and rejects the majority of the W(— > fj,v)+jet 
events (as noted above, for W(— > fj,v)+]et events, IfS^ 1 
is a measure of the transverse momentum of the W bo- 
son). The background at this stage consists primarily of 
dijet events with a misidentified electron and an isolated 
muon from semileptonic b or c quark decay (note that 

the muon momentum contributes to the measured ]ff ). 
This background is effectively eliminated by requiring 
two jets with Et > 15 GeV. At this stage the background 
is a mixture of QCD multijet (including W(-^> (iv)+jct 
events), Z/j* — ► tt — > e/i, and WW — > e/j, events. For 
the remaining stages of event selection, neural network 
techniques are used. 

The optimal discrimination between signal and back- 
ground can be achieved using three separate net- 
works [109]. Each of these discriminates between the 
signal and one of the dominant backgrounds: 



• Network 1 (NN1): tt vs QCD multijet events 

• Network 2 (NN2): tt vs WW -> e/z events 

• Network 3 (NN3): tt vs Z —> tt ^ efi events 

Training is performed on large samples of data (QCD 
multijet) and MC (ti, WW, Z -> tt) events. To reduce 
bias, these samples are prepared with a minimal selection 
criteria of E T > 10 GeV, p T > 10 GeV/c, and A jcts > 1 
with E^ 1 > 10 GeV. From these, a small sub-sample of 
1000-2000 events is selected at random to provide the 
training sample. The number of nodes and the input 
parameters for each network are selected to maximize 
discrimination between signal and background. The best 
results are obtained using three identical networks, each 
with six input nodes, seven hidden nodes, and one out- 
put node. The input parameters, which consist of five 
energy and one topological variable for each of the three 
networks, are listed below: 

• Variables used in NN1 and NN2 

— E T , transverse energy of leading electron 

— _Ey t2 , transverse energy of next to leading jet 

— $r , missing transverse energy as measured 
by the calorimeter 

— H^ tB , scalar sum of jet transverse energies 

H£ ts = E t\ with lV ct | < 2.5 

all jets 

and E^ > 15 GeV 

— M(e/i), electron-muon invariant mass 

— A(j>(efi), azimuthal separation of the leading 
electron and muon 

• Variables used in NN3 

— same as NN1 and NN2 except that E^ n re- 
places Ej, t2 (transverse energy of leading jet) 

Each of the three networks is trained for 2000 training 
cycles. Training is started with a set of random weights 
and thresholds which are adjusted using back propaga- 
tion as the training proceeded. During training the tar- 
get outputs are set to unity for signal and zero for back- 
ground. For simplicity, the outputs of the three networks, 
Onni,Onn2,Onn3, are combined into an overall dis- 
criminant, 

<m b = i | t | — . (di) 

0JVJV1 0JVJV2 OnN3 

which gives the probability that a given event is signal. 
The output from such a combination is equivalent to that 
from a single network that was trained on each of the 
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TABLE 34. Number of observed and expected events passing at each cut level of the e\i neural network analysis. Expected 
number of tt events are for mt = 170GeV/c 2 . Uncertainties correspond to statistical and systematic contributions added in 
quadrature. 



Number of events passing e/j, NN selection 







Total 


Mis-id 


Physics 






Data 


sig + bkg 


bkg 


bkg 


tt 


E* T > 15 GeV, p£ > 15 GeV 












± e id + fi id ± trig 


130 


98 ± 12 


54 ±2 


40 ±9 


4.3 ±0.9 


+ AR(fi,jet) > 0.5, AR(e,fi) > 0.25 


58 


54 ±9 


12 ± 1 


39 ±8 


3.4 ±0.7 


+ $r al > 15 GeV 


44 


42 ±8 


5.8 ±0.5 


32 ±7 


3.3 ±0.7 


+ 2 jets, > 15 GeV 


6 


4.4 ±0.9 


0.68 ±0.17 


0.85 ±0.21 


2.8 ±0.7 


+ O c ^ h > 0.88 


4 


2.5 ±0.7 


0.04 ±0.12 


0.19 ±0.07 


2.3 ±0.5 



TABLE 35. Expected number of signal and background 
events after all cuts in 112.6 pb -1 for the e/i neural network 
analysis. Uncertainties are statistical and systematic contri- 
butions added in quadrature. The systematic uncertainty on 
the total background includes the correlations among the dif- 
ferent background sources. 



Expected number of e/i NN evts in 112.6 pb 



tt MC m (GeV/c 2 ) 




150 


3.51 ±0.86 


160 


2.84 ±0.68 


170 


2.30 ±0.53 


180 


1.81 ±0.41 


Z — ► TT — > 6/i 


0.10 ±0.07 


multijet (mis-id e) 


0.04 ±0.12 


WW -> e/i 


0.08 ±0.02 


DY — » tt — > e/i 


0.01 ±0.01 


Total background 


0.23 ±0.14 



three different backgrounds and the signal [109]. Test- 
ing on independent samples found that a requirement of 
Oj$m > 0.88 maximized the relative significance (which 
is defined to be the ratio of expected number of signal 
events to the measured uncertainty on the number of 
background events). 

After this selection four candidate events remain, three 
of which are also selected by the conventional analysis. 

Backgrounds and acceptances are estimated in much 
the same way as is done for the conventional analysis. 
The only real difference is that an additional correction is 
made for the effect of multiple interactions. This correc- 
tion is obtained by comparing special MC samples with 
one and two minimum bias events added with the stan- 
dard MC samples which have no minimum bias events 
added. The acceptance variation is parameterized as a 
linear function of the number of interactions and a cor- 
rection factor is obtained by applying this function to 
the distribution of the number of interactions throughout 
Run 1. A correction factor of 9% was found for events; 
since the Z/j* and WW backgrounds are kinematically 



Number of Events 






L 


— Signal + Background 

— Total Background 
• Data 








* ' i 1 







Output of the Neural Net (O NN ) 

FIG. 37. Distribution of signal±background (vertical 
hatching), background (diagonal hatching), and data (circles) 
as a function of neural network output. All initial cuts have 
been made except the requirement of 2 jets with Et > 15 
GeV. 

and topologically similar, they receive the same correc- 
tion. The QCD multijet background, being derived from 
data, does not require such a correction. The expected 
numbers of signal and background events passing the full 
selection are given in Table 35. Figure 37 shows a com- 
parison of data and the expected signal and background 
as a function of neural network output after all initial cuts 
except the requirement of 2 jets with Et > 15 GeV. A 
cross section of 9.75±5.18(stat)±1.95(sys) pb is obtained 
for the NN-based analysis which is in agreement with the 
value of 7.1 ±4.8 pb obtained for the conventional analy- 
sis. Comparison of acceptances and background expecta- 
tions between the two analyses finds the NN analysis with 
an increase in acceptance of 10% (for m t = 172 GeV/c 2 ) 
for the same background. Table 34 shows the number 
of data events, expected signal (m t = 170 GeV/c 2 ), and 
expected background surviving at each stage of the selec- 
tion. As with the conventional analysis, good agreement 
is seen between what is observed and what is expected. 

Systematic uncertainties are handled the same way as 
in the conventional analysis and are summarized in Ta- 
ble 36 and, with the exception of the uncertainty on the 
efficiency times branching fraction due to the top quark 
mass, are discussed in Sec. IX. The value of nit measured 
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TABLE 36. Efficiency times branching fraction (e x B) and 
statistical and systematic uncertainties (in percent), and ex- 
pected background and corresponding statistical and system- 
atic uncertainties (in number of events), for the e/i neural 
network analysis. 



only to the Z — > tt background. 



APPENDIX E: TREATMENT OF 
UNCERTAINTIES 







Expected # of Bkg 


events 




£ X B{%) 


Ztt 


WW 


DYtt 


Mis-id e 




0.351 


0.095 


0.077 


0.006 


0.044 


Statistical 


0.004 


0.055 


0.006 


0.004 


0.117 


i_i Limmobi ly 




0.005 


0.004 


0.000 




Energy Scale 


0.035 


0.026 


0.011 


0.004 




e id 


0.010 


0.004 


0.002 


0.000 




High-pT /i id 


0.037 


0.011 


0.008 


0.001 




e+jets trig 


0.020 


0.005 


0.004 


0.000 




MC generator 


0.025 


0.005 


0.004 


0.000 




Top quark mass 


0.020 










Bkg crsec 




0.010 


0.008 


0.001 




Other sim 




0.022 








Mis-id e 










0.003 


Total 


0.065 


0.067 


0.018 


0.005 


0.117 



by D0 is 172.1 ± 7.1GcV/c 2 [92,101] and the central 
value is used in the calculation of the efficiency times 
branching fraction. This uncertainty of ±7 GeV/c 2 is 
composed of an uncertainty of 4.0 GeV/c 2 due to jet en- 
ergy scale, 1.9 GeV/c 2 due to the tt MC generator, and 
6.1 GeV/c 2 due to statistics and other sources. The ef- 
fect of these uncertainties is determined by parameter- 
izing the efficiency times branching fraction as a linear 
function of top quark mass in the region between 165 
and 180 GeV/c 2 . This parameterization is used to con- 
vert the above uncertainties on m t into uncertainties on 
efficiency times branching fraction. The uncertainties on 
m t due to the jet energy scale and tt generator translate 
into uncertainties on e x B of 3.3% and 1.6%, respec- 
tively. These uncertainties are combined with the other 
jet energy scale and tt generator uncertainty contribu- 
tions (described in Sec. IX). The uncertainties on mt 
due to statistics and from other sources translate into an 
uncertainty of 5.8% on e x B, and are included as a sep- 
arate source of uncertainty in Table 36. As discussed for 
the conventional dilepton analyses, there is a significant 
discrepancy between data and MC for the jet Et spec- 
tra in Z+jet events. The conventional analyses correct 
for this by taking the jet cut survival probabilities from 
data and applying them to the MC. Such a procedure 
is not possible with a NN analysis. Fortunately the pri- 
mary disagreement between data and MC is in E^ t2 . It 
is for this reason that the variables used for NN3 differ 
from those for NN1 and NN2 in that NN3 employs E^ n 

,jet2 



instead of E } T . To account for the uncertainty due to 
the initial jet cuts of Nj cts > 2 with Ej? > 15 GeV, a 
data versus MC comparison was made and a difference 
of 21% was found. This uncertainty is listed under the 
category "Other simulation" in Table 36 and is applied 



As shown in Eq. 10.2, calculation of the tt production 
cross section requires as input the number of observed 
events found in all channels, the total expected back- 
ground, the individual channel acceptance for tt events, 
and the integrated luminosity for each channel. To sim- 
plify the discussion, Eq. 10.2 can be written in the form 



<j(m t ) t 



N — B 
AC 



where 



i 
3 

A = — ^ ' AiC{ 

i 



(El) 

(E2) 
(E3) 

(E4) 
(E5) 



with the sum i being over all nine channels and the sum 
j being over all backgrounds. 

It is assumed that the backgrounds and acceptances 
are subject to the same kinds of uncertainties (sec 
Sees. IX A 1 - IX A 20) and that no correlation exists 
among the different uncertainties. With these assump- 
tions, the background error matrix is given by 



s 2 Bij 



2 

Bixij i 



(E6) 



where i and j represent the various backgrounds in 
the different channels {e.g. IF+jets background in 
e+jets/topo channel), /U represents the source of uncer- 
tainty {e.g. electron identification), and the error ma- 
trix for a given uncertainty, Sg^j , is computed from first 
principles according to the equation 



(E7) 



where the symbol represents the average of the en- 
closed quantity when it is varied according to the uncer- 
tainty /i. Accordingly, the correlation matrix for a given 
uncertainty is given by 



C 



X2 



(E8) 



where Ss^i is the uncertainty on background i due to 
source /i. 
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With these definitions, the uncertainty on the total 
background is obtained from Eq. E3 by propagation of 
errors 



OBiJ \8BjJ ~ BlJ 



°Bi 3 



= E* 
= EE< 5 

M ij 



2 

B\xij 



where 



ij 

= J]] CBfj,ijSBfiidBfij- 



(E9) 
(E10) 
(Ell) 
(E12) 

(E13) 
(E14) 



In the case of uncorrelated errors (Cb^ = the unit 
matrix), Eq. E14 reduces to the usual quadratic sum for- 
mula, 



(*b)1 = E<V 



(E15) 



In the case of maximal positive correlation (Cb^ popu- 
lated entirely by l's), Eq. E14 reduces to a linear sum of 
errors, 



(E16) 



For the analyses in this paper, all uncertainties are 
handled according to one of these two limiting cases. Sta- 
tistical uncertainties are handled by the quadratic sum 
formula (Eq. E15) and systematic uncertainties are han- 
dled according to the linear sum formula (Eq. El 6). The 
total uncertainty on the background is therefore 



S 2 



b ■ E E '>■/••< 

fi— stat i 



'Bjii 



(E17) 



The importance of correlations for the background cal- 
culation as a whole depend on the extent to which dif- 
ferent backgrounds are affected by the same systematic 
uncertainties. The systematic uncertainties for all back- 
grounds to all channels are given in Tables 27 - 29. 

Applying the steps above to Eq. E4, the uncertainty 
on A is found to be 



0~A — £2 E E CA l xijO~A l x%b~A i ijC>iC>3 

= £2 E E ^Api^i 
(U=stat i 

+^E (e^Y- 

fi=sys \ i / 



(E18) 



(E19) 



Systematic uncertainties on the acceptance (e x B) are 
given for all channels in Table 31. Note that the accep- 
tance uncertainties are highly correlated due to the fact 
that the calculation for each channel is affected by essen- 
tially the same set of systematic uncertainties. The same 
relative uncertainty on the luminosity has been assumed 
for all channels (see Sec. IXA1). 

The total uncertainty on the top quark cross section is 
obtained by propagation of errors using Eq. El. The four 
inputs to the cross section can, in principle, give rise to 
six different correlation terms. However, the signal (N) 
has only a statistical uncertainty and the uncertainties 
on the acceptance (A) and the integrated luminosity (£) 
are uncorrelated. Therefore, the only correlation terms 
are those between the background (B) and the accep- 
tance (A) and between the background (B) and the inte- 
grated luminosity (£). The corresponding uncertainties 
are given by the equations 



s ba = E ( S ba)1 and 
5lc = J2^c)l- 



(E20) 
(E21) 



The error corresponding to a given uncertainty (/i) is 
calculated from first principles according to the equations 

(S BA )l = (BA), (B) tl (A) ft and (E22) 



(6bc)1=(B£),-(B),(£)„ 



(E23) 



where the symbol (...) represents the average of the en- 
closed quantities when they are varied according to the 
uncertainties \i. The correlation coefficients are given by 



Cba = 



°BA „i r _ °BC 

— — and Cbc — f f ■ 

ObOa ObOc 



(E24) 



In the linear/quadratic approximation, these correlation 
coefficients simplify to 



°BA = CbA ^ ( 5b \ ( 5a \ ' 



SbSa 



and 



Cbc = SbJI ^ Cbc "- {5b) » {5c) » ' 



K- E ^ Vc)> • 



SbSc 
1 



Sb$c 
ob 



(Sb) £ (Sc), 



(E25) 
(E26) 

(E27) 

(E28) 

(E29) 
(E30) 



GO 



The total uncertainty on the top quark cross section is 
therefore given by 



dN 

( ^JL 
V dC 

dftt 
dB 



5% 



Si 



don 
dB 

d<J t i 
dB 



5% 



dA 
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8A 



S 2 ba 



dC 



BC' 



(E31) 



6% 



(N-B) 
5b 



+ S\ + 5l 
2 A 2 L 2 



+ 



N-B 



Cba-j^ + Cbc~£ 



(E32) 
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